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PREFACE 


This  report  presents  the  results  and  final  analysis  of  the  local 
ice  load  measurement  conducted  on  the  four  deployments  aboard  the  USCGC 
POLAR  SEA  between  1982-84.  Data  were  collected  in  first  year  and 
multiyear  level  ice  in  McMurdo  Sound,  Antarctica.  The  first  and  second 
deployment  results  from  trips  to  the  Alaskan  Arctic  as  well  as  the 
instrumentation  and  data  analysis  techniques  were  presented  in  "Ice  Loads 
and  Ship  Response  to  Ice"  (SSC-329)  {reference  1) .  The  third  deployment 
results  from  the  Antarctic  were  presented  in  a  report  to  the  Maritime 
Administration  {Reference  2) .  Results  of  the  fourth  data  collection 
program  from  the  Beaufort  Sea  in  the  summer  of  1984  are  presented  in  "Ice 
Loads  and  Ship  response  to  Ice  -  A  Second  Season"  (SSC-339)  (Reference 
4) .  This  report  summarizes  the  previous  data  collection  programs  and 
provides  the  final  data  analysis  of  all  data  as  a  whole. 

A  statistical  analysis  of  extreme  pressures  and  forces  was 
performed  for  the  data  collected  on  all  four  deployments  and  is  present*  i 
in  this  report.  Pressures  over  one  subpanel,  four  subpanels,  and  forces 
on  frames,  stringers  (as  if  the  ship  were  longitudinally  framed) ,  and  the 
total  load  on  the  panel  were  fitted  to  3  parameter  extreme  value 
distributions.  The  results  of  the  extreme  value  statistics  performed 
were  then  used  to  suggest  ice  load  criteria  in  support  of  icebreaking 
ship  design  and  hull  design  regulations  for  icebreaking  ships. 
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1 . 


INTRODUCTION 


In  1982,  USCGC  POLAR  SEA  was  instrumented  with  an  array  of  strain  gages 
on  the  port  bow  for  the  purpose  of  measuring  ice  impact  pressures.  Two  trips 
to  the  Alaskan  Arctic  were  made  in  October  1982  and  March-April  1983  during 
which  time  about  1400  impact  events  were  collected.  The  research  was  carried 
out  on  behalf  of  the  Interagency  Ship  Structure  Committee,  the  U.S.  Maritime 
Administration,  and  Transport  Canada  (Transportation  Development  Centre) . 

Work  was  performed  in  conjunction  with  environmental  data  collection  programs 
sponsored  by  the  Alaskan  Oil  and  Gas  Association  and  the  U.S.  Maritime 
Administration . 

Ten  cant  frames  (CF  35  to  CF  44)  were  instrumented  at  8  vertical 
locations  by  strain  gaging  the  webs  of  the  frames  in  compression  perpendicular 
to  the  shell  plating  (Figure  1) .  A  total  of  sixty  active  channels  of  strain 
gages  allowed  contact  pressures  over  an  area  of  up  to  98  ft2  (9.1  m2)  to  be 
measured.  An  individual  strain  gage  channel  was  related  to  an  area  of  1.63 
ft2  (.15  m2)  for  which  a  uniform  pressure  was  computed  for  a  measured  strain. 

A  complete  description  of  the  data  acquisition  system  and  the  data  reduction 
procedures  as  well  as  the  results  of  the  two  deployments  can  be  found  in 
Reference  [ 1 ] * . 

The  POLAR  SEA's  trip  to  the  Antarctic  in  January  1984  offered  a  third 
opportunity  to  collect  ice  impact  data  in  thick  level  ice  in  conjunction  with 
resistance  tests  sponsored  by  the  Maritime  Administration  (MARAD) ,  Naval 
Engineering  Division  of  the  U.S.  Coast  Guard  and  Canadian  Transportation 
Development  Centre  (TDC) .  An  additional  310  ice  impact  events  were  collected 
by  this  effort  and  are  reported  under  contracts  to  MARAD  [2]  and  TDC  [3]. 

A  fourth  data  collection  program  was  conducted  in  October  and  November 
of  1984,  termed  the  1984  Summer  Deployment,  to  gather  additional  data  in 
summer  multiyear  ice  conditions  where  the  highest  loads  could  be  expected. 

This  deployment  recorded  337  impact  events  which  are  presented  and  analyzed  in 
SSC-339  [4] .  This  report  summarizes  data  from  all  four  deployments  and 
presents  further  analysis  of  the  complete  data  set. 


*  Numbers  in  brackets  refer  to  references  listed  in  Section  8. 
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Figure  1 

STRAIN  GAGE  LOCATIONS  ON  POLAR  SEA 


2.  SUMMARY  OF  THE  MEASUREMENT  PROGRAMS,  COLLECTED  DATA  AND  ICE  CONDITIONS 


The  local  ice  impact  loads  data  collection  program  has  made  use  of  four 
deployments  of  POLAR  SEA  between  the  fall  of  1982  and  the  fall  of  1984  to 
acquire  data  in  different  geographical  areas.  Seven  data  sets  are  identified 
by  geographical  area  and  date  of  data  collection.  The  data  sets,  representing 
2039  individual  impact  events,  are  summarized  in  Table  1.  A  listing 
summarizing  the  extremes  of  each  event  for  each  data  set  can  be  found  in 
Appendix  A,  sorted  by  the  highest  average  single  sub-panel  pressure.  Actual 
routes  of  the  ship  or  operating  areas  where  the  data  were  collected  are  shown 
on  the  maps  in  Figures  2  and  3. 

For  two  of  the  data  sets  involving  ice  conditions  that  included  both 
first-year  and  multi-year  ice,  it  was  possible  to  identify  subsets  of  known 
multi-year  impacts.  Sixty-seven  known  multi-year  known  multi-year  events  were 
identified  in  the  North  Chukchi  Winter  83  data  which  included  the  dedicated 
rams  of  multi-year  ridges  described  in  Reference  1.  An  additional  32  known 
multi-year  events  were  identified  in  the  Summer  3eaufort  84  data  set.  The 
multi-year  subsets  are  summarized  at  the  bottom  of  Table  1.  It  should  be 
noted  that  many  more  multi-year  events  occurred  and  were  recorded  in  the  South 
Chukchi  Winter  83,  North  Chukchi  Winter  83  and  Beaufort  Summer  84  data  sets, 
however  specific  multi-year  events  could  not  be  identified.  The  Beaufort 
Summer  82  data  were  collected  at  a  time  when  only  multi-year  ice,  with  the 
exception  of  light  refreeze,  existed  in  the  area. 

Table  2  presents  combined  data  sets  that  were  grouped  according  to  ice 
conditions  first  and,  secondly,  according  to  geoarapnic  area.  The  data  sets 
were  grouped  to  provide  the  largest  collection  of  data  of  similar  conditions 
for  the  extreme  value  analysis  of  Section  4.  In  this  type  of  analysis,  larger 
data  sets  provide  improved  extrapolation  to  the  longer  return  periods. 
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Figure  2 

LOCATIONS  OF  DATA  COLLECTION  EFFORTS  IN  THE  ALASKAN  ARCTIC 
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Figure  3 

LOCATION  OF  ANTARCTIC  ICE  LOADS  DATA  COLLECTION 


TABLE  1 


SUMMARY  OF 

DEPLOYMENTS,  DATA  SETS  AND  ICE  COND 

ITIONS 

TITLE  &  OATE 

LOCATION 

ICE 

TYPE 

NO  OF 
EVENTS 

Beaufort  Summer  82 

Sep  28  -  Oct  16 

100-150  nm  north  of  Prudhoe  3ay 
in  the  Alaskan  Beaufort  Sea 

MY 

167 

S  Bering  Winter  83 

Mar  24  -  Mar  26 

Transit  from  St. Paul  Is.  to  the 
west  end  of  St. Lawrence  Is.  in 
the  Bering  Sea 

FY 

173 

N  Bering  Winter  83 

Mar  27  -  Mar  28 

Transit  from  St.  Lawrence  Is.  to 
the  Bering  Strait  in  the  Bering  Sea 

FY 

241 

S  Chukchi  Winter  83 

Mar  29  -Apr  2 

Apr  28  -  May  2 

Transit  from  the  Bering  Strait  to 
Point  Hope  in  the  Chukchi  Sea  and 
return 

FY  ,MY 

299 

N  Chukchi  Winter  83 

Apr  3  -  Apr  27 

Round  trip  transit  Point  Hope  to 
Wainwright  in  the  Chukchi  Sea, 
operation  off  Wainwright 

FY  ,MY 

513 

Antarctic  Summer  84 

Jan  9  -  Jan  13 

McMurdo  Sound,  break-in  to  McMurdo 
Base 

FY 

309 

Beaufort  Summer  84 

Nov  18  -  Dec  1 

Operation  between  Barter  Is.  and 
Barrow  in  the  Beaufort  Sea,  transit 
through  the  Chukchi  Sea  to  the 
Bering  Strait 

FY  ,MY 

337 

SUBSETS  OF  KNOWN  MULTI-YEAR  EVENTS 

N  Chukchi  Winter  83  MY 
Apr  3  -  Apr  27 

North  Chukchi  Sea  off  Wainwright 

MY 

67 

Beaufort  Summer  84  MY 
Nov  12  -  Dec  1 

Beaufort  and  Chukchi  Seas 

MY 

32 

7/-f 


TABLE  2 


SUMMARY  OF  COMBINED  DATA  SETS 


TITLE 

COMBINED  FROM 

ICE  NO  OF 

TYPE  EVENTS 

Known  Multi-Year 

Beaufort  Summer  82 

N  Chukchi  Winter  83  MY 

Beaufort  Summer  84  MY 

MY 

266 

Heavy  Mixed  FY  &  MY 

Beaufort  Summer  82 

S  Chukchi  Sea  83  N  Chukchi  Sea  33 
Beaufort  Sutraier  84 

FY  ,MY 

1017 

Known  First-Year 

S  Bering  Winter  83 

N  Bering  Winter  83 

Antarctic  Summer  84 

FY 

723 

Sumner  Beaufort  Sea 

Beaufort  Summer  82 

Beaufort  Summer  84 

mostly  MY 

504 

Winter  Chukchi  Sea 

S  Chukchi  Winter  83 

N  Chukchi  Winter  33 

FY  ,MY 

398 

Winter  Bering  Sea 

S  Bering  Winter  83 

N  Bering  Winter  83 

FY 

812 
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3. 


ANALYSIS  OF  PRESSURE  VERSUS  CONTACT  AREA  RELATIONSHIPS 


The  intent  of  this  section  is  to  investigate  specific  relationships  that 
may  exist  among  the  data  collected.  Section  3.1  examines  an  actual  event  to 
illustrate  the  nature  of  the  measured  pressures.  Section  3.2  looks  at 
average  pressure  over  the  contact  area.  Three  of  the  data  sets  (North  Chukchi 
Winter  83,  South  Bering  Winter  83  and  Antarctic  Summer  84)  are  used  to 
illustrate  three  major  operating  scenarios;  high  Arctic  with  old  ice, 
first-year  ice  near  the  ice  edge  and  first-year  thick  level  ice.  Section  3.3 
shows  the  influence  (or  lack  there  of)  of  impact  force  on  velocity. 

3.1  Pressure  Imprint.  Descriptions 


Ice  pressure  is  calculated  at  each  time  step  within  an  event  by 
multiplying  the  sixty  measured  strains  by  a  60  x  60  matrix  to  produce  the 
sixty  average  ice  pressures  over  the  sub-panel  areas.  Results  from  each  event 
are  saved  in  the  form  of  an  impact  pressure  time-history.  Figure  4  shows 
these  calculated  pressures  for  five  sequential  time-steps  for  an  event  taken 
from  the  North  Chukchi  Winter  83  data  set.  Sampling  occurred  32  times  per 
second  so  the  time-step  shown  is  .031  seconds.  The  values  printed  are  in  psi 
(145  psi  =  1  MPa)  and  are  arranged  in  the  same  manner  as  the  sub-panel  areas 
on  the  hul 1 . 

This  event  took  place  during  April,  1983  in  the  North  Chukchi  Sea, 
resulting  in  a  peak  pressure  of  1141  psi  (7.9  MPa).  To  illustrate  the 
impact,  all  values  above  100  psi  (0.69  MPa)  were  highlighted.  The  sub-panels 
are  approximately  square  so  this  event  has  a  length  to  height  ratio  of  about 
4.  Part  of  the  event  may  be  below  the  panel  which  would  reduce  this  ratio. 

A  few  negative  values  can  be  seen  both  near  the  imprint  and  on  the 
"quiet"  portions  of  the  panel.  Two  factors  account  for  the  negative  values. 
One  is  a  shift  in  the  zeroes  of  all  channels  due  to  thermal  effects.  New 
zeroes  could  only  be  taken  when  there  was  no  load  on  the  panel  and  often  this 
was  not  possible.  Negative  zero  shifts  result  in  measured  values  below  the 
true  value.  Near  the  impact,  negative  values  could  also  result  from  an 
assumption  in  the  data  reduction  algorithm.  The  algorithm  assumes  that  the 
impact  pressure  is  uniform  over  the  sub-panel  area.  If  the  actual  impact  is 
concentrated  over  a  smaller  area  than  the  sub-panel,  the  uniform  pressure  for 
that  sub-panel  will  be  over-predicted  and  the  adjacent  sub-panels  will  be 
under-predicted  (negative  if  there  was  no  actual  load  on  them).  The  two 
effects  cancel  and  are  therefore  not  expected  to,  cause  any  significant  errors 
(i.e.  less  than  10  percent). 

Software  was  developed  during  the  1984  Antartic  deployment  to  correct  the 
thermal  drift  problem.  This  involves  viewing  and  zeroing  each  strain  gage 
time-history  prior  to  data  reduction.  All  data  collected  after  1983  has 
employed  this  method  as  part  of  the  data  reduction  process.  Additionally, 
revision  of  the  data  reduction  matrix  to  include  the  effects  of  non-uniform 
sub-panel  loading  has  been  studied.  While  reduction  of  all  the  data  sets  a 
second  time  with  an  improved  matrix  would  improve  the  accuracy  of  the 
predicted  pressures,  the  improvement  would  not  affect  the  final  results 
significantly  and  would  not  warrant  the  effort  required. 
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EVENT  RECORDED  IN  THE  NORTH  CHUKCHI  SEA 
ON  APRIL  24,  1983  AT  16:11:59 

Note:  Values  above  100  Psi  <.69MPa)  are  highlighted 


Figure  4 

EXAMPLE  OF  AN  ICE  PRESSURE  IMPRINT 
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3.2  Average.  Pressure  versus  Contact  Area 

There  are  many  ways  to  plot  pressure  versus  area.  In  this  section  only 
the  average  pressure  over  the  total  contact  area  will  oe  considered.  This 
should  not  be  confused  with  the  average  pressure  over  some  smaller  area  within 
the  contact  area.  In  this  section,  data  points  associated  with  a  single  sub¬ 
panel  area  imply  that  the  total  contact  area  was  only  one  sub-panel  area  for 
that  event  and  so  forth  for  larger  areas.  Higher  average  pressures  over  a 
portion  of  the  contact  area  are  not  plotted.  It  is  important  not  to  confuse  a 
plot  such  as  Figure  5  with  one  such  as  Figure  10  that  plots  highest  average 
pressures  within  the  contact  zone. 

Looking  at  average  pressure  and  contact  area  helps  to  understand  the 
mechanics  of  the  impact  event.  To  make  use  of  this  data  for  design,  it  is 
useful  to  know  the  contact  area  associated  with  a  given  force  (which  implies  a 
certain  average  pressure  over  the  contact  area).  In  this  section,  the  pressure 
versus  height  (line  loads  on  transverse  frames)  and  pressure  versus  length 
(line  loads  on  longitudinal  framing)  will  also  be  presented. 

Figure  5  shows  the  average  pressure  versus  contact  area  for  3  of  the  data 
sets.  Figure  5a  is  from  the  North  Chukchi  Sea  83  data  and  represents  the  high 
Arctic  with  considerable  multi-year  ice.  Clearly  the  large  imprints  have 
lower  average  pressures.  At  86  ft 2  (8  m2),  the  pressure  tends  to  cluster 
around  60  psi  (0.4  MPa).  Figures  5b  and  5c  show  similar  data  for  the  South 
Bering  Sea  83  and  Summer  Antarctic  84  data.  The  relatively  light  first-year 
ice  conditions  of  the  South  Bering  Sea  produced  low  average  pressures  over  the 
contact  area.  The  Antarctic  data  taken  in  thick  first-year  ice  fell  between 
the  other  two  data  sets. 

The  same  data  can  be  viewed  as  a  plot  of  force  versus  contact  area  shown 
in  Figure  6.  The  highest  force  in  the  North  Chukchi  Winter  83  data  set  (Figure 
6a)  occurred  at  only  43  ft 2  (4  m2).  One  would  expect  the  highest  forces  to  be 
associated  with  the  largest  contact  areas,  in  general,  as  the  other  data  sets 
show  (Figure  6b  and  6c).  Possibly  there  were  insufficient  high  energy  impacts 
in  hard  ice  to  generate  both  large  forces  and  contact  areas  at  the  same  time. 
The  randomness  of  the  ice  properties  could  also  be  responsible.  Figures  6b  and 
6c  show  a  much  clearer  trend  of  increasing  force  with  contact  area,  however. 

Figures  7  and  8  show  the  average  pressure  as  a  function  of  vertical  and 
horizontal  extent.  These  values  are  useful  in  predicting  the  loads  on  local 
framing  members.  Later  sections  will  further  examine  the  implications  to 
desi gn. 
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AVERAGE  PRESSURE  VERSUS  TOTAL  CONTACT  AREA 


FORCE  FORCE  FORCE 


LT 
500  r 


400 
300  f- 
200 

100  - 
oL 


MN 


20 


40  60 

AREA 


80 


100 


300 


LT  MN 


AREA 


LT  MN 


20 


40 

AREA 


60 


80 


M‘ 


Ft 


.  NORTH  CHUKCHI 
WINTER  '83 


.  SOUTH  BERING 
WINTER  '83 


.  ANTARCTIC 
SUMMER  '84 


Figure  6 

TOTAL  FORCE  VERSUS  TOTAL  CONTACT  AREA 
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Figure  7 

AVERAGE  PRESSURE  VERSUS  HEIGHT 
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AVERAGE  PRESSURE  VERSUS  LENGTH 
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3.3  Force  versus  Velocity 


Analytical  models  of  ice  impact  mechanics  predict  a  clear  relationship 
between  force  and  velocity  [5,  6].  Figure  9a  shows  the  force  versus  velocity 
data  for  the  POLAR  SEA  for  the  North  Chukchi  Winter  83  data  set.  The  maximum 
force  for  each  event  (one  data  point  per  event)  is  plotted  against  initial 
impact  speed.  No  clear  trend  with  velocity  is  evident  in  the  data.  Figure  9b 
shows  the  same  data  for  the  Antarctic  Summer  84  data  set.  In  the  latter  case, 
there  does  appear  to  be  an  increase  of  force  with  increasing  velocity,  however 
the  trend  could  be  masked  by  the  fact  that  the  data  were  collected  in  varying 
ice  thicknesses  from  3  to  6  feet  (1  to  2  m)  and  the  actual  ice  thickness  for 
an  individual  impact  is  not  known.  In  the  former  case,  the  ice  conditions 
included  a  range  of  first-year  and  multi-year  floes  witn  widely  varying 
thicknesses.  With  such  a  range  of  ice  conditions  and  the  tendency  for 
operators  to  be  more  cautious  in  heavier  ice,  it  is  not  surprising  that  no 
clear  trend  was  round.  It  is  evident  that  both  high  and  low  forces  occurred 
at  all  speeds.  Speed  control  (othe'*  than  that  which  was  already  imposed) 
would  not  have  lowered  the  force  levels. 

Two  general  conclusions  can  be  drawn  from  this  data.  One  is  the  need  to 
collect  ice  properties  data  in  as  much  detail  as  possible.  Secondly,  unless 
detailed  ice  geometry  and  properties  data  exist,  only  statistical  analysis  of 
the  data  is  valid.  Statistical  values  can  be  derived  from  the  data  that 
describe  the  ship  impact  process  in  an  overall  sense.  This  work  has  been  done 
and  is  presented  in  Section  4. 
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4.  STATISTICAL  ANALYSIS  OF  THE  DATA 


4.1  The  Shape,  of  the  Ice.  Impact.  Pressure.-Area  Curve,  and,  .its.  Effect  on  Ice 
Load  Development 


During  the  course  of  data  reduction  and  analysis  of  four  deployments  of 
collected  data  comprising  over  2000  impact  events,  a  multitude  of  pressure 
versus  area  plots  were  produced.  These  include  plots  for  each  time  step 
throughout  an  event,  plots  at  the  time  of  peak  pressure  and  peak  force  during 
each  event,  and  plots  of  the  highest  recorded  pressure  over  each  area  for  all 
events  in  a  data  set.  One  significant  fact  emerges;  the  pressure-area  curve 
has  a  consistent  and  character*; sti :  shape  whether  it  is  for  one  instant  of  one 
impact  or  the  extreme  envelope  of  many  impacts. 

The  implications  of  this  finding  have  a  profound  effect  on  simplifying  a 
statistically  based  ice  impact  load  algorithm.  The  pressure-area  curve  is 
typically  plotted  on  a  log-log  scale  as  shown  in  Figure  10.  At  small  impact 
areas,  pressures  tend  toward  a  line  which  decreases  slightly  with  increasing 
area.  At  large  impact  areas,  the  average  pressure  becomes  force  limited  and 
tends  toward  a  line  proportional  to  the  reciprocal  of  area.  The  upper  line  or 
the  pressure  asymptote  has  a  constant  slope  on  this  type  of  plot  that  is 
determined  by  area  to  a  power  in  the  range  of  -0.2  and  -0.3.  Events  occur 
randomly  and  the  effect  is  to  shift  the  asymptotes  up  or  down,  or  to  the  left 
or  right,  depending  on  the  severity  of  the  impact  and  the  type  of  ice 
encountered.  If  the  average  pressure  over  a  small  area  is  independent  of  the 
total  force  during  an  impact,  each  can  be  predicted  statistically  from 
measured  data  to  determine  the  asymptotes  of  the  pressure-area  curve.  The 
complete  design  curve  at  all  areas  can  therefore  be  generated  from  two  lines 
of  constant  slope,  one  associated  with  a  limiting  pressure  and  one  with  a 
1 imiti ng  force. 

The  slope  of  the  limiting  pressure  line  was  determined  from  a  number  of 
sources.  First,  envelope  pressure-area  curves  from  a  number  of  measured  data 
sets  are  shown  in  Figure  11  for  a  wide  range  of  areas.  As  one  can  see  a  line 
of  area  to  the  -0.2  power  fits  the  limits  of  the  data  well.  Secondly,  Figure 
12  presents  an  analysis  of  the  extreme  events  from  the  1982  Summer  Beaufort 
Sea  deployment.  Data  from  the  pressure-area  curves  for  these  events  up  to  6 
sub-panels  was  non-dimensional ized  and  plotted  in  the  figure;  that  is,  the 
ratio  of  the  highest  average  pressure  over  3  sub-panels  to  that  over  1  sub¬ 
panel  is  plotted  against  the  ratio  of  the  areas  (3),  for  instance.  Most 
combinations  of  sub-panel  ratios  are  considered  in  the  plot  up  to  6  sub¬ 
panels.  One  can  see  that  the  upper  bound  of  the  data  approaches  a  limiting 
slope  of  area  to  the  -0.2  power. 

There  is  an  apparent  correlation  with  theory  as  well.  If  it  is  assumed 
that  the  impact  pressure  (P)  is  proportional  to  compressive  ice  strength  [12] 
and  ice  strength  (ct)  is  proportional  to  strain  rate  (  k)  to  a  fractional  power 
b  [13],  then  it  can  be  shown  that  pressure  is  proportional  to  area  to  the  -b 
power  for  certain  shaped  indenters.  Unconfined  crushing  strength  shows  this 
behavior  at  low  strain  rates.  At  the  high  strain  rates  normally  encountered 
in  ship/ice  impacts,  unconfined  crushing  strength  becomes  constant,  but 
triaxially  crushing  strength  continues  to  exhibit  this  relationship  with  the 
same  fractional  power  as  the  unconfined  case  at  lower  strain  rates. 
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AREA  (Square  Meters) 


Figure  10 

A  TYPICAL  EVENT  SHOWING  THE  ASYMPTOTIC 
NATURE  OF  ICE  IMPACTS 
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Figure  12 

NON-DIMENSIONAL  ANALYSIS  OF  PRESSURE- AREA  DATA  FROM 
THE  BEAUFORT  SUMMER  82  DATA  SET 


Consider  a  spherical  shaped  indenter  where  the  strain  rate  at  the  failure 
zone  is  proportional  to: 

.  u 

e  = 

z 


where  u  is  the  indentation  velocity  and 
z  is  the  indentation  distance 

A 

z  =  - _ for  small  z 

2  7T  R 


where  A  is  the  contact  area  and 

R  is  the  radius  of  the  indenter. 


Therefore: 


P  ~  oq  =  e!5  «  u*3  A“b 


A  similar  result  occurs  if  it  is  assumed  that  the  strain  rate  is  proportional 
to  the  indenter  velocity  divided  by  the  contact  area. 

Cox,  et  al.  [13]  presents  uniaxial  compressive  strength  data  for  .ulti-year 
ice  from  the  Alaskan  Beaufort  Sea  tested  at  different  temperatures.  The  mean 
of  the  data  taken  at  23  degrees  Fahrenheit  (-5  degrees  Centigrade)  gives  a 
value  of  0.209  for  b.  This  compressive  strength  data  should  be  typical  of  the 
strengths  of  the  multi-year  ice  that  generated  the  ice  impact  loads.  The  slope 
of  the  pressure  limit  line  predicted  by  this  method  is  again  very  close  to  the 
-0.2  shown  above.  This  is  not  to  say  that  the  complex  interaction  of  impact 
of  a  ship's  side  with  ice  can  be  directly  compared  with  spherical  indentation 
or  triaxial  crushing  tests,  but  this  development  does  show  an  interesting 
correlation  to  the  measured  results. 
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4.2  Regression  o.f  Extreme.  Value.  Distributions 


Individual  impacts  from  all  deployments  consist,  in  reduced  form,  of 
time-histories  of  the  average  pressure  over  each  of  the  sixty  sub-panel  areas 
within  the  load  panel.  During  data  reduction,  each  sub-panel  time-history  is 
scanned  to  identify  the  time  of  highest  average  pressure  over  any  of  the  sub¬ 
panels  during  the  event  and  the  time  of  peak  force  on  the  entire  panel.  The 
corresponding  pressure  and  force  are  also  noted.  Pressure  versus  area  relation¬ 
ships  are  developed  for  both  of  these  times  during  the  event  by  identifying 
the  highest  pressure  on  any  sub-panel  at  that  time,  checking  all  contiguous 
areas  for  the  next  highest  pressure  and  so  forth  until  all  loaded  areas  have 
been  identified.  Since  each  sub-panel  area  is  the  same,  1.63  ft 2  (0.15  m2), 
the  result  is  a  tabular  listing  of  the  decay  in  average  pressure  as  a  function 
of  the  number  of  sub-panels  that  are  contiguously  loaded  at  the  time  of  peak 
pressure  and  peak  force  for  each  event.  The  number  of  sub-panels  can  easily 
be  multiplied  by  the  sub-panel  area  to  produce  plots  of  pressure  versus  area 
at  a  given  instant  of  time  as  shown  in  Figure  10. 

Additionally,  the  highest  load  along  a  frame  or  stringer  was  also  computed 
since  the  sub-p^nel  width  was  the  frame  spacing  (16  in  or  400  mm)  and  the 
sub-panel  height  was  almost  the  frame  spacing  (14.7  in  or  375  mm).  For  the 
load  on  the  frame,  the  highest  average  pressure  on  a  single  sub-panel  was 
located  first,  then  the  highest  average  pressure  over  two  adjacent  sub-panels 
arranged  vertically  one  above  the  other,  and  then  three  adjacent  sub-panels  in  a 
vertical  line.  The  process  continued  for  each  number  of  sub-panels  up  to  six 
in  a  vertical  line,  the  total  height  of  the  array  of  sub-panels  in  the  bow 
panel.  The  force  for  each  was  computed  as  the  average  pressure  times  the 
corresponding  measurement  area.  The  force  remained  relatively  constant 
regardless  of  the  length  of  the  measurement  area  (high  pressures  over  short 
lengths  and  low  average  pressures  over  longer  lengths)  but  the  maximum  typically 
occurred  at  a  length  of  about  half  the  panel  height.  The  fact  that  the  force  is 
relatively  constant  for  all  frame  lengths  allows  a  single  value  of  force  to  be 
used  to  characterize  each  event.  The  fact  that  the  maximum  force  on  a  frame 
occurred  at  a  span  roughly  half  the  height  of  the  panel  means  that  the  limited 
panel  height  was  sufficient  to  capture  the  maximum  load  on  the  frame  and  should 
not  effect  the  extreme  value  analysis.  A  similar  process  was  done  for  adjacent 
sub-panels  in  a  horizontal  line,  assuming  the  ship  was  longitudinally  framed. 

The  force  versus  stringer  length  for  up  to  10  sub-panels  arranged  in  a 
horizontal  line,  the  bow  panel  length,  was  computed  and  the  highest  was  saved 
for  the  extreme  value  analysis. 

The  above  described  procedures  have  been  performed  as  part  of  the  data 
reduction  process  for  each  event  on  each  deployment.  The  statistical  analysis 
conducted  for  and  described  in  this  report  starts  with  this  data  as  well  as  the 
peak  force  on  the  entire  panel  for  each  event  as  its  basis.  The  highest  average 
pressures  over  one  and  four  sub-panels  from  the  pressure-area  curve  and  the  peak 
forces  on  a  frame  and  stringer  as  well  as  the  force  over  the  total  panel  were 
analyzed  statistically.  Events  were  divided  into  data  sets  based  on  geographic 
area  of  operation  or  ice  conditions  as  described  in  Section  2.  Each  of  the  five 
variables  was  identified  for  the  events  in  each  set  of  data  and  ranked  from 
highest  to  lowest.  The  corresponding  probability  of  non-exceedance  was  computed 
based  on  the  formula: 


Probability  =  1-I/(N+1) 

where  I  is  the  rank  of  the  variable  in  the  data  set  and  N  is  the  number  of 
events  in  the  data  set. 
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A  three  parameter  extreme  value  distribution  was  then  fit  to  the  pressure 
or  force  versus  probability  data  [14].  The  curvature  parameter  in  the  extreme 
value  curve  fit  is  an  indication  of  the  type  of  asymptotic  distribution  as 
described  in  Appendix  B.  Three  types  of  distributions  are  possible  [15];  a 
Gumbel  (Type  I)  which  is  unbounded  and  linear  on  extreme  value  probability 
plots,  a  Frechet  (Type  II)  which  is  unbounded  and  linear  on  log-extreme  value 
plots,  and  a  Weibull  (Type  III)  which  has  an  upper  bound.  The  trends  in  the 
parameters  of  the  distribution  and  therefore  the  type  of  distribution  with 
respect  to  ice  conditions  have  been  studied  and  indicate  that  the  more  severe 
ice  conditions  have  a  Frechet  type  distribution.  As  ice  conditions  decrease 
in  severity,  a  Gumbel  type  extreme  value  distribution  appears  to  be  more 
appropriate.  First-year  ice  conditions  often  exhibit  an  upper  bound  in 
extreme  loads,  indicating  a  Weibull  distribution  gives  the  best  fit  to  the 
data.  Appendix  C  gives  plots  of  the  extreme  value  distributions  for  the  five 
variables  studied  as  well  as  the  corresponding  three  parameter  curve  fits. 
Tables  of  the  curve  fit  parameters  are  shown  in  Appendix  0.  Correlation 
coefficients  for  the  curve  fits  were  typically  0.98  or  higher  indicating  an 
excellent  fit  to  the  data. 

Figure  13  shows  the  single  sub-panel  pressure  data  from  the  North  Chukchi 
Winter  83  data  set  as  an  example.  Also  plotted  is  the  Gumbel  distribution 
computed  graphically  for  386  of  these  events  in  the  1983  report  describing 
this  data  collection  [1].  As  one  can  see,  the  three  parameter  curve  fit 
indicates  a  Frechet  (Type  2)  extreme  value  distribution,  due  to  the  upward 
curvature.  The  curve  fit  of  reference  [1]  is  sligntly  below  the  full  data  set 
because  only  a  portion  of  the  data  was  used  in  the  initial  analysis. 
Extrapolation  to  longer  return  periods  using  the  reference  [1]  curve  fit  would 
under-predict  the  extreme  pressure  substantially  compared  with  the  three 
parameter  curve  fit. 

Figures  14  and  15  show  a  comparison  of  the  highest  average  pressure  over 
one  sub-panel  and  the  highest  force  on  the  entire  panel  extreme  value 
distributions,  respectively,  for  different  geographical  areas  in  the  Alaskan 
Arctic.  Figures  16  and  17  show  corresponding  distributions  for  different  ice 
conditions.  The  relative  increase  in  the  severity  of  the  ice  loads  with 
increasing  latitude  and  increasing  severity  of  ice  conditions  is  apparent  in 
the  figures. 
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HIGHEST  AVERAGE  SINGLE 
SUB-PANEL  PRESSURE 


Figure  13 

EXAMPLE  OF  THE  3  PARAMETER  CURVE  FIT  OF  THE 
NORTH  CHUKCHI  WINTER  83  DATA  SET 


?? 


PROBABILITY  OF  NON-EXCEEDANCE 


Figure  14 

COMPARISON  OF  SINGLE  SUB-PANEL  PRESSURE  FOR 
DIFFERENT  GEOGRAPHICAL  AREAS 
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Figure  15 

COMPARISON  OF  TOTAL  PANEL  FORCE  FOR 
DIFFERENT  GEOGRAPHICAL  AREAS 
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Figure  16 

COMPARISON  OF  SINGLE  SUB-PANEL  PRESSURE  FOR 
DIFFERENT  ICE  CONDITIONS 


Figure  17 

COMPARISON  OF  TOTAL  PANEL  FORCE  FOR 
DIFFERENT  ICE  CONDITIONS 


4.3  Application  to  a  Design.  Procedure 


Section  4.1  described  the  asymptotic  nature  of  the  pressure  versus  area 
curve,  for  a  single  instant  of  time,  observed  from  the  measured  data.  That 
section  also  showed  that  the  slope  of  the  pressure  asymptote  appears  to  fall 
within  a  narrow  range  of  values  related  to  the  crushing  strength  versus  strain 
rate  behavior  of  the  impacted  ice.  If  the  slope  for  the  pressure  asymptote  is 
defined  by  the  type  of  ice,  then  the  extreme  value  distribution  for  single 
sub-panel  pressure  can  be  used  to  estimate  the  highest  expected  pressure  in  a 
given  return  period  for  one  contact  area  (1.63  ft  2  or  .15  m 2)  and  thus 
define  the  location  of  that  asymptote  on  the  pressure-area  curve.  The  force 
asymptote  can  be  similarly  determined  from  the  extreme  value  distribution  of 
force  using  the  same  return  period.  Since  average  impact  pressure  equals  the 
impact  force  divided  by  the  contact  area,  the  force  asymptote  plots  as  a  45 
degree  line  on  a  log-log  pressure-area  curve  (see  Figure  13). 

But  which  force  distribution  should  be  used?  The  extreme  value  distribu¬ 
tions  described  in  Section  4.2  and  given  in  Appendix  C  are  a  characterization 
of  the  forces  over  very  specific  areas  ana,  additionally,  very  specific  shaped 
areas.  The  total  panel  force  is  measured  over  the  entire  instrumented  panel; 
an  area  of  98  ft  2  or  9.1  m2.  The  panel  dimensions  are  7.3  feet  high  by  13.3 
feet  long  (2.24  by  4.07  m).  The  force  on  a  frame  data  used  in  the  distribu¬ 
tion  are  the  highest  force  for  each  event  computed  from  the  pressure  versus 
length  along  a  frame  described  in  the  previous  section.  This  force  is  con¬ 
sidered  to  act  over  a  vertical  strip  of  hull  plating  16  inches  (400  mm)  wide 
and  up  to  the  height  of  the  panel  high.  Similarly,  the  force  on  a  stringer  is 
assumed  to  act  over  a  horizontal  strip  along  the  hull  that  is  14.7  inches 
(375  mm)  wide  and  up  to  the  length  of  the  panel  in  length.  For  local  ice 
impact  load  development,  l.e.  the  determination  of  loads  for  plating,  frames 
and  stringers,  the  force  on  a  frame  for  transversely  framed  ships  and  the 
force  on  a  stringer  for  longitudinally  framed  ships  should  be  used.  The  force 
distribution  must  be  consistent  with  the  loaded  area  of  the  scantlings  for 
which  it  will  be  used. 

The  distributions  of  total  force  on  the  panel  and  highest  average 
pressure  over  four  sub-panels  also  given  in  Appendix  C  are  included  for 
completeness.  The  total  force  distribution  indicates  the  magnitude  that  ice 
forces  can  reach  for  bow  contact  areas  up  to  98  ft  2  (9.1  m2).  The  panel  did 
not  measure  total  bow  force,  however,  since  some  of  the  shell  plating  was 
obviously  loaded  outside  the  instrumented  panel.  Large  area  loads  such  as 
those  that  might  be  useful  for  design  of  girders,  decks  and  bulkheads  are 
therefore  best  determined  either  from  global  load  measurements  or  analytical 
models  that  estimate  total  bow  force.  The  measured  data  indicate  that  a 
contact  area  corresponding  to  four  sub-panels  is  about  where  the  intersection 
of  the  force  and  pressure  asymptotes  occur  (see  Figure  13).  This  is  the  part 
of  the  pressure-area  curve  where  the  actual  pressure  deviates  most  from  the 
asymptotes.  The  four  sub-panel  distributions  give  an  indication  of  the 
magnitude  of  this  deviation,  therefore. 
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To  summarize,  the  extreme  value  distributions  of  the  measured  data  can  be 
used  to  develop  a  pressure  versus  area  curve  that  describes  the  highest 
expected  ice  impact  pressure  in  a  given  return  period  for  the  range  of  impact 
areas  associated  with  local  hull  scantlings.  The  data  presented  in  Section  4.2 
is  only  appropriate  to  hull  forms  of  similar  size  and  shape  to  that  of  the 
POLAR  Class  and  framing  systems  of  similar  spacing  to  that  described  here. 
Application  of  this  approach  to  other  ship  designs  will  be  discussed  in  the 
following  section. 


AREA  (Square  Feet) 

I _ i  i  j  i  i _ i _ i_x_i _ i _ I  i  i  i  i  i  1  ) 

0.1  05  1  5  10 

AREA  (Square  Meters) 


Figure  18 

EXAMPLE  OF  HOW  A  PRESSURE-AREA  CURVE 
CAN  BE  GENERATED  FROM  THE 
STATISTICAL  DISTRIBUTIONS 
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5.  RECOMMENDATIONS  FOR  LOCAL  ICE  LOAD  DESIGN  CRITERIA 


5.1  Bow  Structure.  Load.  Criteria 


It  is  the  opinion  of  the  authors  that  two  conditions  or  return  periods 
should  be  considered  when  developing  the  ice  loads  for  design;  the  normal 
operating  condition  (loads  in  the  range  of  one  to  three  year  return  period) 
and  the  survival  condition  (lifetime  loads).  Normal  operating  loads  should 
cause  no  deterioration  in  the  ship's  operating  performance  while  survival 
loads  may  cause  some  loss  of  performance  but  not  catastrophic  failure. 

To  develop  the  loads  for  each  loading  condition,  the  expected  number  of 
impacts  must  be  estimated  for  the  time  period  by  conducting  an  operational 
assessment  of  the  ship  in  the  ice  conditions  in  which  it  will  operate.  Table  3 
gives  a  summary  of  the  frequency  of  impacts  for  the  collected  data  to  assist 
in  this  estimate.  The  reciprocal  of  the  number  of  impacts  expected  is  the 
probability  of  non -exceedance  used  to  enter  the  figures  in  Appendix  C  for  the 
loads.  The  single  sub-panel  pressure  establishes  one  point  on  the  pressure 
asymptote  (see  Figure  18).  The  maximum  force  on  a  frame  or  stringer  distribu¬ 
tions  establish  the  force  asymptote,  depending  on  whether  the  ship  is  trans¬ 
versely  or  longitudinally  framed,  respectively.  It  should  be  noted  that  the 
normal  operating  condition  will  be  associated  with  a  relatively  small  number 
of  impacts  which  fall  within  or  just  beyond  (slight  extrapolation  from)  the 
measured  data.  This  is  not  usually  the  case  with  survival  loads;  the  number 
of  impacts  may  require  a  large  extrapolation  from  the  existing  data  base  of 
measured  loads.  There  is  typically  a  much  higher  confidence  in  the  normal 
operating  loads,  therefore. 

An  example  will  be  presented  to  help  illustrate  the  proposed  procedure. 
Assume  that  an  icebreaker  is  being  designed  for  operation  in  the  Northern 
Bering  Sea  and  is  expected  to  be  underway  there  for  two  winter  months  out  of 
the  year.  It  is  expected  to  operate  12  hours  per  day  during  this  time.  The 
annual  number  of  impacts  is  estimated  to  be: 

8.2  impacts/hr  X  12  hr/day  X  60  days/yr  =  5904  impacts/yr 

See  Table  3  for  measured  impacts  per  hour  in  different  operating  areas. 

If  the  normal  operating  loads  are  taken  as  those  expected  annually,  the 
probability  to  enter  the  distributions  for  N  Bering  Winter  83  (Appendix  C)  is 
one  minus  the  reciprocal  of  5904  or  .99983.  The  graphs  of  the  distributions 
can  be  used  directly  or,  more  accurately,  the  equation  for  the  three  parameter 
curve  fit  can  be  used: 


Result  =  [l-(-ln(Probabil  ity))*-]  (A2/C)  +  A1 


The  coefficients 

from  Tables  D1  (single 

sub-panel 

pressure)  and 

1  D4  (force  on  a 

frame)  are: 

C 

A1 

A 2 

RESULT 

Pressure 

.026 

289  (1.99) 

84  (0.58) 

942  (6.49) 

psi  (MPa) 

Force  on  a  Frame 

-.239 

36  (0.36) 

11  (0.11) 

356  (3.55) 

LT  (MN) 
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These  results  are  shown  graphically  in  Figure  19.  The  authors  recommend  using 
the  average  pressure  over  an  area  equal  to  the  frame  spacing  squared  as  the 
design  pressure  for  plating.  For  a  16  inch  (400  mm)  frame  spacing,  the 
pressure  is  slightly  less  than  the  result  shown  above  since  the  sub-panel 
measurement  area  was  slightly  smaller  than  the  frame  spacing  squared.  The 
results  can  be  scaled  by  the  factor  [(16  x  16) /( 144  x  1.63)]-0.2  =  .983  to 
obtain  the  plating  design  over  the  frame  spacing  squared  (926  psi  or  6.38  MPa 
for  this  example).  Design  pressures  for  frame  design  can  be  taken  from  Figure 
19  for  this  example  using  an  area  of  the  frame  spacing  times  the  loaded  length 
along  the  frame.  A  loaded  length  equal  to  the  frame  span  will  result  in  a  low 
uniform  pressure  over  the  entire  span  while  choosing  a  shorter  loaded  length 
will  result  in  a  higher  uniform  pressure  over  a  shorter  length,  presumably  the 
limiting  design  condition. 

The  measurement  panel  that  recorded  the  data  presented  in  the  previous 
chapters  was  located  in  the  bow  of  the  POLAR  SEA.  Load  criteria  based  on  the 
measured  data  can  only  be  developed  for  the  bow  therefore.  Extension  of  these 
loads  to  other  areas  of  the  ship  will  be  discussed  in  the  next  section. 


TABLE  3 


FREQUENCY  OF 

OCCURRENCE 

OF  IMPACT 

EVENTS  FOR  THE 

MEASURED 

DATA 

DATA  SET 

TYPICAL 
THRESHOLD 
(  uc) 

AVERAGE 

IMPACT 

FREQUENCY 

(events/hr) 

RECORDING 

TIME 

(hrs) 

TOTAL 

ELAPSED 

TIME 

(hrs) 

NUMBER 

OF 

IMPACTS 

Beaufort  Summer  32 

250 

3.2 

52.2 

314.4 

167 

S  8ering  Winter  83 

75 

10.5 

16.5 

29.0 

173 

N  Bering  Winter  33 

120 

8.2 

29.5 

48.5 

241 

S  Chukchi  Winter  83 

120 

4.4 

68.0 

206.5 

299 

N  Chukchi  Winter  83 

150 

3.6 

143.0 

617.0 

513 

Antarctic  Summer  84 

100 

21.0 

15.0 

15.0 

309 

Beaufort  Summer  84 

150 

4.9 

68.8 

325.5 

337 
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AREA  (Square  Meters) 


Figure  19 

PRESSURE-AREA  CURVE  FOR  THE  EXAMPLE 
(Northern  Bering  Sea  Icebreaker) 


5.2  Extension  of  the  Criteria,  to.  Other  Areas,  of.  the.  Ship 


Measurement  of  hull -ice  impact  loads  has  concentrated  on  the  bow  of 
icebreaking  ships  since  this  is  where  the  highest  local  loads  occur.  The  bow 
area  is  normally  considered  to  extend  to  the  point  of  maximum  beam  though  some 
reduction  in  impact  pressures  is  expected  near  the  shoulders  due  to  the 
relatively  low  angles  of  incidence  with  the  ice.  Section  5.3  will  discuss  the 
effect  of  hull  shape  on  the  local  ice  loads.  This  section  will  deal  with  the 
areas  aft  of  the  forwardmost  point  of  maximum  beam  which  are  generalized  as 
the  amidship  area  and  the  stern  area.  Specific  ice  impact  loads  have  not  as 
yet  been  measured  in  these  areas  of  an  icebreaking  hull  form.  Guidelines  that 
exist  are  based  on  theory  or  experience  or  both. 

The  stern  area  will  be  addressed  first  since  it  can  more  directly  be 
related  to  bow  loads.  Previous  sections  stated  that  the  pressure  asymptote  of 
the  pressure-area  curve  appears  to  be  largely  independent  or  very  weakly  a 
function  of  impact  speed,  however  the  force  asymptote  is  linearly  related  to 
impact  speed.  One  would  therefore  expect  that  if  the  maximum  astern  speed  was 
some  percentage  of  the  maximum  forward  speed  in  heavy  ice  conditions,  say  30 
to  50  percent,  then  tne  force  on  a  frame  or  stringer  values  would  be  reduced 
by  a  similar  factor  relative  to  the  bow  forces.  This  results  in  a  shift  of 
the  force  asymptote  to  the  left  on  the  pressure-area  curve  as  shown  in  Figure 
19.  The  pressure  asymptote  remains  the  same  indicating  that  average  impact 
pressur_s  over  very  small  areas  are  expected  to  be  equivalent  to  pressures 
over  the  same  area  at  the  bow.  This  probably  means  that  frame  design  loads 
will  be  reduced  more  than  plating  design  loads  relative  to  the  bow  (presuming 
that  the  bow  plating  is  determined  by  the  pressure  asymptote).  An  astern 
shape  similar  to  the  bow  shape  has  been  assumed.  While  this  may  be  a  good 
assumption  for  a  conventional  icebreaker  like  the  POLAR  Class,  transom-sterned 
ships  such  as  many  of  the  icebreaking  supply  boats  that  have  been  built 
recently  will  probably  require  additional  corrections  for  local  hull  angles. 

For  the  amidship  area,  the  problem  is  more  complex.  The  ice  impact  speed 
normal  to  the  hull  in  the  midship  area  is  small  for  straight  ahead 
icebreaking.  While  the  normal  velocity  can  increase  in  turns  or  when  the  ship 
is  maneuvering,  the  limiting  load  is  almost  certainly  due  to  a  pressured  ice 
condition.  Pressured  ice  conditions  can  occur  when  the  ship  is  moving, however, 
the  greatest  large  area  average  pressures  are  most  likely  when  the  ship  is 
beset.  Several  theoretical  solutions  have  been  advanced  for  the  load  per  unit 
length  that  can  be  developed  in  an  ice  sheet  under  pressure.  These  are 
grouped  in  two  categories  (summarized  in  References  16,  17  )  known  as  ridge 

building  forces,  the  load  that  causes  failure  in  the  ice  sheet,  and  limiting 

driving  force,  the  load  that  can  be  developed  due  to  wind  and  current  loads. 
The  maximum  force  per  unit  length  that  exists  in  an  ice  sheet  is  the  minimum 
of  the  two  solutions;  that  is,  if  the  driving  force  causing  the  pressure 
exceeds  the  load  carrying  capacity  of  the  ice  sheet,  then  the  ice  sheet  will 

fail  limiting  the  load  to  the  failure  value.  Several  of  the  solutions  are 

shown  in  Figure  20  for  the  load  on  a  multi-year  floe  in  a  three  foot  thick 
pressured  ice  sheet.  It  should  be  noted  that  driving  force  is  limited  by  the 
first-year  ice  cover,  the  weakest  part  of  an  ice  field  that  may  contain 
multi-year  floes  of  much  greater  thickness.  The  same  loads  could  be  seen  by  a 
ship  beset  in  these  conditions. 
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Figure  20  shows  that  loads  range  from  20  LT/ft  (0.66  MN/m)  to  over  100  LT/ft 
(3.3  MN/m)  depending  on  the  loaded  length  and  the  thickness  of  the  ice  cover. 
But  these  loads  are  global  in  nature,  those  that  are  expected  over  the  length 
of  the  waterline.  Local  impact  pressures  recommended  by  various  classification 
societies  for  the  midship  area  are  approximately  sixty  percent  of  bow  design 
pressures.  Canadian  Arctic  Shipping  Pollution  Prevention  Rules  (CASPPR)  [18] 
recommend  63  percent,  the  USSR  Rules  [19]  recommend  60  percent  and  the  new 
American  Bureau  of  Shipping  (ABS)  Rules  [20]  recommend  50  to  60  percent  for 
icebreakers  depending  on  ice  class. 
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DRIVING  FORCE 


MN/m  LT/ft 


MULTIYEAR  FLOE  DIAMETER 


Figure  20 

ESTIMATES  OF  LIMITED  DRIVING  FORCE  ON  A  MULTIYEAR  FLOE 
SURROUNDED  BY  3  FOOT  THICK 
FIRST  YEAR  ICE 
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5.3  Areas,  for  Improving  the.  Load.  Criteria 


Certainly  one  area  for  improving  ice  load  criteria  is  measurement  of 
local  ice  loads  in  locations  other  than  at  the  bow.  Section  5.2  discussed  ways 
to  apply  the  data  base  of  measured  bow  loads  to  other  areas  of  the  hull  which 
are  based  on  theoretical  development  and  experience.  Measured  loads, 
especially  in  the  midship  area,  would  not  only  verify  these  methods  but  also 
provide  a  better  understanding  of  the  ice-structure  interaction  as  well.  The 
slow  speed  impact  loads  that  might  be  experienced  by  the  midbody  of  a  ship  in 
the  beset  condition  also  have  application  to  offshore  structures  and, 
conversely,  measured  data  on  offshore  structures  might  be  useful  in  studying 
the  ship  problem. 

Two  other  areas  of  improvement  must  be  considered  to  properly  address  the 
full  range  of  icebreaking  ship  sizes  and  hull  impact  locations.  The  first  is 
the  aspect  of  the  effect  of  local  hull  angles  on  the  resulting  loads.  Both 
the  A3S  and  Russian  Rules  incorporate  a  theoretical  solution  of  the  effect  of 
local  hull  angles.  Figure  9.2  of  the  new  ABS  Rules  [20]  is  reproduced  in 
Figure  21.  The  figure  shows  the  variation  in  non-dimensional  force  with 
changing  waterline  half-angle  (a)  and  local  section  angle  (  g)  measured  from 
the  vertical.  The  values  of  these  variables  at  the  location  of  the 
measurement  panel  on  the  bow  of  POLAR  SEA  were  a  =  30  degrees  and  g  =  54 
degrees.  The  resulting  non-dimensional  load  factor  is  0.9.  One  can  obtain  an 
estimate  of  the  expected  loads  at  a  location  with  different  hull  angles  by 
reading  the  load  factor  for  the  location  from  Figure  21,  dividing  the  factor 
by  .9,  and  then  multiplting  the  result  by  the  statistically  derived  pressure 
of  Section  4  appropriate  to  the  ice  conditions  being  considered.  Experimental 
verification  of  the  underlying  theory  would  increase  one's  confidence  in  the 
validity  of  the  predicted  results,  however. 

The  second  area  is  the  effect  of  the  ship's  displacement  and  power  on  the 
resulting  loads.  Johansson  [21]  proposed  a  relationship  for  ice  impact 
pressure  (P)  as  a  function  of  the  ship's  power  (N)  and  displacement  (a)  as: 

P  =  P0  +  c(N  A)b 

where  P0  and  c  are  constants  and  b  is  0.5  in  his  work  developing  the  Finnish- 
Swedish  Rules  [22]  and  later  modified  the  exponent  b  to  0.33  for  Arctic  LNG 
tankers  [23],  In  his  discussion  of  the  latter,  Tunik  [24]  recommended  the 
expression : 

P  =  P0  +  c  NO. 18  a0.05 

based  on  the  work  of  Kurdyumov  [25].  More  recently,  in  his  work  developing  the 
new  ABS  Rules,  Tunik  has  adopted  an  expression  of  the  form: 

P  =  c  N0-2  a0-15 
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As  more  data  is  collected  from  ships  of  different  sizes  and  with 
different  available  power,  it  will  be  important  to  try  to  validate  these 
functional  relationships.  The  task  is  very  difficult  however.  Variations  in 
ice  conditions  between  different  sets  of  measured  loads  can  easily  mask  the 
effects  of  power  and  displacement.  It  is  important  to  note  that  care  must  be 
used  in  scaling  measured  ice  loads  by  these  relationships  as  well,  since  the 
total  available  power  is  not  always  employed  when  the  data  are  being 
collected.  This  is  particularly  true  in  the  case  of  the  POLAR  SEA  due  to  the 
flexibility  of  her  propulsion  plant.  Each  of  her  three  shafts  can  be  driven 
in  diesel -electric  mode  up  to  6000  HP  or  gas  turbine  mode  up  to  20,000+  HP. 
She  is  often  operated  in  some  combined  configuration  which  could  lead  to 
erroneous  results  if  scaling  is  based  on  total  installed  horsepower  of  60,000 
HP. 


Figure  21 

COEFFICIENTS  FOR  NON-DIMENSIONAL  PRESSURE  VARIATION 

WITH  LOCAL  HULL  ANGLES 
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6.  RECOMMENDATIONS  FOR  RESPONSE  CRITERIA 


6.1  Plating  Response 


A  statistical  description  of  the  ice  pressures  and  forces  has  been 
adopted.  The  previous  section  presented  extreme  value  distributions  of  the 
loads,  extreme  pressures  or  forces  versus  the  probability  of  non-exceedance. 
For  a  new  design,  the  number  of  impacts  expected  in  a  given  return  period  is 
estimated  from  an  operational  assessment  of  the  ship's  intended  mission. 

Annual  or  several  year  loads  might  be  considered  operational;  that  is,  those 
that  are  resisted  with  no  deterioration  in  the  ship's  performance.  A  second 
return  period,  the  desired  lifetime  of  the  ship,  can  also  be  considered.  Loads 
for  this  time  period  represent  survival  loads;  the  ship  should  resist  these 
with  some  probable  reduction  in  performance  but  no  catastrophic  failure. 
Obviously,  response  criteria  for  the  two  loading  conditions  will  vary  with  the 
ship  and  its  intended  use.  A  typical  approach,  however,  might  be  to  allow  a 
small  amount  of  permanent  set  for  the  plating  and  only  elastic  response  of  the 
framing  under  normal  operating  loads.  Much  more  permanent  set  and  even  some 
plasticity  in  the  framing  could  be  considered  under  survival  conditions. 

A  variety  of  theories  and  response  criteria  are  available  in  the 
literature  that  could  be  considered  for  plastic  response  of  plates.  Figure  22 
shows  a  comparison  of  plate  response  equations  for  the  range  from  purely 
elastic  response  to  rupture.  The  curves  in  the  figure  are  computed  for  high 
tensile  steel  assuming  a  yield  strength  of  50,000  psi  (345  MPa),  an  ultimate 
strength  of  85,000  psi  (586  MPa)  and  a  16  in  (400  mm)  frame  spacing.  The 
extreme  left  side  of  the  figure  (curve  14)  shows  an  estimate  of  the  load  to 
cause  rupture  assuming  a  strain  of  14  percent  [30],  The  right  side  of  the 
figure  (curve  1)  shows  the  load  to  cause  the  yield  stress  assuming  purely 
elastic  response.  Between  are  plotted  a  number  of  equations  for  plastic 
response  to  various  levels  of  permanent  set  and,  to  the  right  side  of  the 
figure,  the  current  response  criteria  of  ABS  (curve  2),  the  USSR  (curve  3), 
and  Canada  (curve  4).  It  should  be  noted  that  all  three  of  these  response 
criteria  use  an  equation  of  the  form: 

t  =  C  s  /P/ 

where  is  the  yield  stress,  P  is  the  uniform  pressure,  s  is  the  frame 
spacing  or  plate  span,  and  C  is  a  coefficient  that  falls  between  the  values  of 
.707  for  an  elastic  response  (curve  1)  and  .5  for  formation  of  three  plastic 
hinges  in  a  rigid-perfectly  plastic  model  of  an  infinitely  long  plate  with 
fixed  end  conditions  (curve  7).  The  ABS  and  USSR  rules  add  an  additional  0.236 
in  (6  mm)  to  the  resulting  thickness  (The  value  is  less  for  other  than  the 
highest  ice  class)  which  accounts  for  their  shift  toward  the  elastic  response 
curve. 
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Figure  22 

COMPARISON  OF  EQUATIONS  DESCRIBING  PLATE  RESPONSE 

UP  TO  RUPTURE 


(1)  Elastic  response,  fixed  end  conditions,  infinitely  long  plate  [26] . 

(2)  ABS  Rules,  0.236  in  (6  mm)  corrosion  allowance  included  [20], 

(3)  USSR  Rules,  0.236  in  (6  mm)  corrosion  allowance  included  [19] . 

(4)  CASPP  Rules,  no  corrosion  allowance  required  [18] . 

(5)  Clarkson,  elastic-plastic,  fixed  end  conditions,  infinitely  long 
plate,  permanent  set  of  the  same  order  as  welding  deflections  [27] . 

(6)  Hughes,  elastic-plastic,  clamped,  ends  free  to  pull  in,  aspect  ratio  of 
0.133  [28]. 

(7)  Johansson  &  Jones,  perfectly  plastic,  fixed  ends,  infinitely  long  plate, 
formation  of  three  plastic  hinges  [21,  29] . 

(8)  Chiu,  Haciski  and  Hirsimaki,  STAGS  finite  element  model,  fixed  end 
conditions,  infinitely  long  plate,  permanent  set  is  0.3  percent  of 
span  [30] . 

(9)  Chiu,  Haciski  and  Hirsimaki,  STAGS  finite  element  model,  fixed  end 
conditions,  infinitely  long  plate,  permanent  set  is  0.5  percent  of 
span  [30] . 

(10)  Jones,  perfectly  plastic,  fixed  ends,  infinitely  long  plate, 
permanent  set  equal  to  plate  thickness  [29] . 

(11)  Jones,  perfectly  plastic  using  the  ultimate  stress,  fixed  ends, 
infinitely  long  plate,  permanent  set  equal  to  10  percent  of  the 
thickness  [29] . 

(12)  Membrane  response,  edge  springs,  infinitely  long  plate,  permanent  set  of 
10  percent  of  plate  span  (Appendix  E) . 

(13)  Chiu,  Haciski  and  Hirsimaki,  membrane,  pinned  end  conditions,  infinitely 
long  plate,  permanent  set  is  10  percent  of  span  [30] . 

(14)  Chiu,  Haciski  and  Hirsimaki,  membrane,  pinned  end  conditions,  infinitely 
long  plate,  permanent  set  is  23  percent  of  span,  estimate  of  rupture 
using  a  strain  of  14  percent  [30] . 
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Also  shown  ara  Clarkson's  (curve  5)  and  Hughes'  (curve  6)  elastic-plastic 
solutions.  Curve  5  assumes  an  infinitely  long  plate  with  permanent  set  similar 
to  that  experienced  during  fabrication  process.  Hughes'  equation  allows  for 
edge  pull-in,  which  he  and  Ratlaff  and  Kennedy  [31,  32]  state  is  important  for 
stout  (small  span  to  thickness  ratio)  plates  typical  of  icebreakers.  Hughes' 
equation  includes  the  effect  of  plate  aspect  ratio  which  can  also  be 
significant  as  plate  length  becomes  smaller.  Curve  6  uses  an  aspect  ratio  of 
.133  or  a  plate  length  of  10  ft  (3  m)  for  this  example,  since  permanent  set  is 
included  in  the  equation  only  if  the  aspect  ratio  is  non-zero.  A  permanent  set 
equal  to  0.5  percent  of  the  span  has  been  used  for  curve  6.  Chiu  et  al . 
recommend  a  maximum  acceptable  permanent  set  of  1  percent  of  the  span  [33]. 
Further,  they  state  that  a  reasonable  criterion  for  permanent  set  is  the 
fabrication  fairing  criteria  which  is  in  the  range  of  0.3  to  0.5  percent  of 
the  plate  span  for  icebreakers.  The  same  paper  also  computes  the  load  versus 
permanent  set  for  a  variety  of  plate  thicknesses  that  might  be  used  in 
icebreakers  using  a  finite  element  model  with  non-linear  material  properties. 

A  two-dimensional  model,  where  the  plate  has  pinned  restraints  at  the  frame 
locations,  is  used  and  the  loading  conditions  include  uniform  loads  over  all 
and  alternating  frame  bays.  For  the  case  where  all  frame  bays  are  loaded 
(fixed  end  conditions),  results  are  presented  here  in  curves  8  and  9  for 
permanent  set  equal  to  0.3  and  0.5  percent  of  the  plate  span,  respectively. 

Coburn,  et  al .  [34 ]}  suggest  a  response  criteria  of  permanent  set  equal  to 
the  pate  thickness  for  extreme  loads  on  icebreakers  using  Jones'  equation  [29] 
developed  from  rigid-perfectly  plastic  analysis  of  a  plate  with  fixea  end 
conditions.  For  an  infinitely  long  plate,  the  equation  takes  the  form  of  that 
presented  above  in  this  section  with  C  equal  to  .354  (curve  10  in  Figure  22). 
Though  curve  10  compares  well  with  the  finite  element  results  of  curves  3  and 
9,  the  amounts  of  assumed  permanent  set  are  quite  different.  Curves  8  and  9 
use  0.3  and  0.5  percent  of  the  plate  span  which  is  approximately  3  to  16 
percent  of  the  plate  thickness  for  the  range  of  thicknesses  considered. 
D'Olivera  [35]  states  that  Jones  equation  is  conservative  for  stout  plates  and 
suggests  using  the  ultimate  stress  or  the  average  between  the  ultimate  stress 
and  the  yield  stress  to  give  a  more  reasonable  result.  If  Jones'  equation  is 
used  with  the  ultimate  stress  (85,000  psi  or  586  MPa)  and  permanent  set  equal 
to  ten  percent  of  the  thickness,  the  result  (the  dashed  line,  curve  11,  of 
Figure  22)  falls  within  tne  finite  element  results  of  Chiu  et  al .  Jones' 
equation  has  the  advantage  of  being  simpler  to  use  in  design  studies  since  it 
is  a  straight-forward,  explicit  equation: 


P/ct  =  4  (t/s)2  (1  +  w/t)  2  w  <  t 


where  w  is  the  amount  of  permanent  set. 
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Appendix  E  presents  a  development  of  a  plate  response  equation  assuming 
pure  membrane  response  with  edge  pull-in.  The  edge  pull-in  is  resisted  by 
springs  equivalent  in  stiffness  to  the  surrounding  plating  and  permanent  set 
of  ten  percent  of  the  plate  span  is  assumed  (curve  12  of  Figure  22). 

Appendix  t  snows  tnat  stout  piui.es  uevelu^  edge  hinges  at  small  amounts  of 
permanent  set  using  the  work  of  Hughes.  For  the  relatively  large  amount  of 
permanent  set  assumed  in  curve  12,  the  plate  will  exhibit  membrane  behavior  and 
edge  pull-in  will  be  resisted  by  the  surrounding  plate.  This  equation  is 
compared  in  the  figure  to  the  case  of  pure  membrane  response  with  fixed  end 
conditions  and  the  same  permanent  set  (curve  13)  [30]. 

Figure  22  illustrates  the  wide  range  of  reponse  equations  and  criteria 
that  could  be  applied  to  the  icebreaker  design  problem*  It  dees  seem  that 
with  a  better  understanding  of  ice  impact  loaas,  a  somewhat  less  conservative 
response  criteria  could  be  used  resulting  in  substantial  weight  savings  in  the 
hull  structure.  There  is  a  lack  of  test  data  however  on  stout  plates  which 
could  be  used  to  validate  the  various  theories.  Some  lata  is  presented  oy 
Chiu,  et  al.  [30.  33].  The  authors  feel  that  aaditiona!  test  data  on  stout 
plates  is  requi  red,  however,  to  fully  validate  these  met.nods. 

Additionally,  certain  equations  presented  here  are  appropriate  for  small 
amounts  of  permanent  set  such  as  Hughes  work  and  the  finite  element  work  of 
Chiu  et  al .  Small  permanent  set  (0.5  percent  or  less  of  the  plate  span)  is 
the  most  appropriate  response  criteria  for  normal  operating  loads  where 
permanent  deflections  are  held  to  normal  building  fairness  tolerances.  Extreme 
or  lifetime  loads  might  be  used  with  a  much  less  conservative  response  criteria 
such  as  10  percent  of  the  plate  span  since  there  is  still  a  great  amount  of 
load  capacity  in  the  plate  before  rupture.  Jones'  equation  or  tnese  of 
Appendix  E  are  intended  to  be  used  with  larger  amounts  of  permanent  set. 
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6.2  Frame  Response  Criteria 


While  various  amounts  of  permanent  set  have  been  considered  for  plating 
(lesser  amounts  for  normal  operating  loads  and  greater  amounts  for  extreme 
loads),  framing  should  generally  be  designed  elastically.  When  one  considers 
plastic  design  of  the  plating,  the  assumption  is  that  the  edges  of  the  plate 
do  not  deflect.  Standard  practice  is  to  design  framing  to  a  maximum  stress 
both  in  shear  and  bending  which  is  determined  by  applying  a  factor  of  safety 
to  both  the  yield  stress  in  bending  and  the  maximum  shear  stress  at  yield 
[36].  Some  additional  specific  considerations  are  necessary  for  icebreakers. 
The  heavily  loaded  hull  plating  and  framing  in  icebreakers  requires  specific 
attention  to  a  framing  system  that  is  forgiving  (not  prone  to  buckling  such  as 
a  trussed  support  system)  and  resistant  to  local  crippling  failure  of  the 
frame  webs  and  tripping  of  the  frames.  The  latter  requires  the  use  of 
tripping  brackets  and  aligning  the  frames  perpendicular  to  the  hull  plating  as 
much  as  the  structural  arrangement  permits. 

Framing  tends  to  be  governed  by  the  shear  loads,  at  least  in  high  Arctic 
icebreakers.  The  frame  span  can  then  be  determined  by  the  bending  moment  after 
the  frame  depth  as  been  determined  by  shear  force  requirements.  As  shown  in 
Section  4,  the  highest  average  pressures  on  a  frame  quickly  reach  the  force 
asymptote  as  the  span  of  the  frame  is  increased,  indicating  that  a  short  frame 
carries  essentially  the  same  load  as  a  longer  one.  Analytical  methods  such  as 
a  finite  element  analysis  of  a  typical  frame  can  refine  the  response 
predictions  and  therefore  produce  a  more  refined  structure  than  is  possible 
with  simple  beam  theory. 

It  seems  prudent  to  use  the  existing  practice  of  elastic  fr.  design 
with  a  factor  of  safety  for  normal  operating  loads.  Extreme  loads  could  be 
considered  with  smaller  or  no  factor  of  safety  or  even  small  amounts  of 
permanent  set  as  long  as  the  framing  system  can  resist  these  loads  in  shear 
and  bending  without  catastrophic  collapse  in  some  other  mode  of  failure.  It 
should  be  noted  that  the  thick  plating  typically  employed  in  icebreakers 
provides  substantial  increases  in  the  plastic  section  modulus  of  frames  over 
the  elastic  one  unlike  conventional  ship  structure.  An  excellent  paper  by 
Varsta  [37]  treats  plastic  response  of  an  icebreaker  frame  and  provides 
insight  into  design  for  other  modes  of  failure  besides  shear  and  bending. 
Plastic  design  of  framing  can  only  be  used  with  confidence  after  a  better 
understanding  of  the  strength  of  the  combined  system  of  framing  and  hull 
plating  is  developed.  Finite  element  models  of  grillage  systems  for 
icebreakers  may  help  in  providing  a  better  understanding  of  this  problem. 
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7.  SUMMARY  OF  CONCLUSIONS  AND  RECOMMENDATIONS 


The  data  collection  and  analysis  efforts  of  this  study  have,  over  the 
past  five  years,  produced  a  solid  data  base  of  measured  loads  and  procedure  to 
use  the  data  for  design.  Measured  loads  have  been  gathered  in  all  operating 
areas  of  the  Polar  Class  including  thick  level  ice  in  the  Antarctic  and  both 
first-year  and  multi-year  ice  in  the  Arctic  (2039  recorded  ice  impact 
time-histories).  This  study  was  the  first  to  produce  detailed  results  of  both 
the  spatial  and  temporal  variation  of  local  ice  impact  pressures,  and  has 
V  greatly  improved  understanding  of  the  ship/ice  interaction  process.  Specific 
"conclusions  from, the  study  are  as  follows: 

'  >'  ■  '  f  y  /<■*:•••  V  ''I  , 

•  Both  peak  force  and  peak  pressures  during  an  impact  increase  with  ice 
severity  (ice  thickness  and  ice  strength).  :  - 

•  In  the  Arctic,  operation  at  higher  latitudes  increases  ice  severity  and 
therefore  ice  loads.  •  ') 

•  Peak  pressure  during  an  impact  appears  to  be  only  weakly  dependent  on 
impact  speed,  ano  no  dependancy  was  discernable  from  the  measured  data. 

•  PeaK  force  during  an  impact  does  increase  with  impact  speed  and  a  linear 

relationship  between  them  appears  reasonable.  '  ) 

•  Total  force  increases  with  increasing  contact  area  but  average  pressure 
decreases  with  increasing  contact  area.  (  i?) 

-s  /  -  :a> 

•  The  average  pressure  distribution  on  the/hull  within  the  impact  zone  at 
an  instant  in  time  is  asymptotic  to  a  line  of  constant  force  at  large 
areas  and  to  a  line  proportional  to  a  negative  fractional  power  of  area 
(in  the  range  of  -0.2  to  -0.3)  at  small  areas.  The  latter  asymptote 
appears  related  to  the  triaxial  crushing  strength  versus  strain  rate 
dependancy  of  the  ice.  ' 

•  Extreme  value  distributions  of  ice  force  and  pressure  for  the  most  severe 
ice  conditions  show  Frechet  type  or  upward  curving  distributions. 
Intermediate  ice  conditions  follow  a  Gumbel  or  linear  type  distribution. 
Only  the  distributions  for  data  recorded  in  light  first-year  ice 
conditions  appear  to  be  bounded  or  Weibul  type  distributions. 

•  It  is  possible  to  develop  a  design  pressure  versus  area  curve  for  local 
loads  based  on  the  extreme  value  distributions  of  the  measured  loads  and 
an  operational  assessment  of  the  ship's  location.  The  design  curve  is 
best  suited  to  ships  similar  in  size  and  shape  to  the  Polar  Class, 
however,  extension  to  ships  of  other  sizes  and  shapes  is  possible  using 
theoretical  and  empirical  methods  for  scaling  the  loads. 

•  The  statistical  design  approach  is  a  rational  one  given  the  randomness  of 
ice  impact  geometries  and  ice  properties. 

•  Two  return  period  or  load  cases  should  be  considered  with  the  statistical 
design  approach;  a  one  to  three  year  normal  operating  load  and  a  lifetime 
survival  load.  The  ship  should  resist  the  normal  operating  load  with  no 
deterioration  in  performance  and  resist  survival  loads  with  no 
catastrophic  failure. 
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•  Plastic  plating  response  criteria  should  be  considered  to  reduce  hull 
weight  given  the  high  local  loads.  Small  amounts  of  permanent  set 
similar  to  the  fabrication  fairing  criteria  can  be  considered  for  normal 
operating  loads  while  considerable  amounts  of  permanent  deformation  can 
be  considered  for  survival  loads. 

•  Plastic  design  of  frames  in  shear  or  bending  should  be  considered  very 
cautiously  and  only  if  all  other  failure  mechanisms  (tripping,  web 
crippling,  etc.)  are  carefully  considered.  Permanent  set  in  frames  is 
probably  only  acceptable  for  survival  loads,  and  elastic  design  for 
normal  operating  loads  is  more  prudent. 


The  authors  feel  that  tnis  study  provides  a  good  basis  by  documenting  and 
providing  understanding  of  local  ice  impact  loads.  It  is  recommended  that,  if 
additional  research  is  done  in  this  area,  the  focus  of  chat  work  be  directed 
toward  gathering  full-scale  data  on  ships  of  other  sizes  and  shapes  ana  also 
collecting  data  at  different  areas  other  than  the  bow.  A  systematic  test  of 
load  versus  permanent  set  for  small  span  to  thickness  plates  would  be  useful 
in  validating  the  various  response  therories  appropriate  to  icebreakers. 
Finally,  testing  or  analysis  of  thick  plate  grillage  systems  would  provide 
insight  to  the  amount  of  plastic  deformation  that  can  be  considered  in  frame 
response. 
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APPENDIX  A 

SUMMARY  OF  MEASURED  DATA  RANKED  BY 
SINGLE  SUB-PANEL  PRESSURE 


KEY: 

PM1  -  Maximum  single  sub-panel  pressure  (psi ) 
PAI  -  Average  pressure  over  the  contact  area  at 
the  time  of  peak  pressure  (psi) 

A1  -  Contact  area  at  the  time  of  peak  pressure 
(sub-panels) 

FI  -  Total  panel  force  at  the  time  of  peak 
pressure  (LT) 

PM2  -  Maximum  single  sub-panel  force  (psi) 

PA2  -  Average  pressure  over  the  contact  area  at 
the  time  of  peak  force  (psi) 

A2  -  Contact  area  at  the  time  of  peak  force 
(sub-panels) 

F2  -  Peak  total  panel  force  (LT) 

VEL  -  Ship  velocity  at  impact  (kts) 


A-l 


TABLE  1 


SUMMARY  OF  DEPLOYMENTS,  DATA  SETS  AND  ICE  CONDITIONS 


TITLE  &  DATE 

LOCATION 

ICE 

TYPE 

NO  OF 
EVENTS 

Beaufort  Summer  82 

Sep  28  -  Oct  16 

100-130  nm  north  of  Prudhoe  Bay 
in  the  Alaskan  Beaufort  Sea 

MY 

167 

S  Bering  Winter  83 

Mar  24  -  Mar  26 

Transit  from  St. Paul  Is.  to  the 
west  end  of  St. Lawrence  Is.  in 
the  Bering  Sea 

FY 

173 

N  Bering  Winter  83 

Mar  27  -  Mar  28 

Transit  from  St.  Lawrence  Is.  to 
the  Bering  Strait  in  the  Bering  Sea 

FY 

241 

S  Chukchi  Winter  83 

Mar  29  -Apr  2 

Apr  28  -  May  2 

Transit  from  the  Bering  Strait  to 
Point  Hope  in  the  Chukchi  Sea  and 
return 

FY ,  MY 

299 

N  Chukchi  Winter  83 

Apr  3  -  Apr  27 

Round  trip  transit  Point  Hope  to 
Wainwright  in  the  Chukchi  Sea, 
operation  off  Wainwright 

FY ,  MY 

513 

Antarctic  Summer  84 

Jan  9  -  Jan  13 

McMurdo  Sound,  break-in  to  McMurdo 
Base 

FY 

309 

Beaufort  Summer  84 

Nov  1 8  -  Dec  1 

Operation  between  Barter  Is.  and 
Barrow  in  the  Beaufort  Sea,  transit 
through  the  Chukchi  Sea  to  the 

Bering  Strait 

FY ,  MY 

337 

SUBSETS  OF  KNOWN  MULTIYEAR  EVENTS 

N  Chukchi  Winter  83  MY 
Apr  3  -  Apr  27 

North  Chukchi  Sea  off  Wainwright 

MY 

67 

Beaufort  Summer  84  MY 
Nov  1 2  -  Doc  1 

Beaut  or t  and  Chukchi  Seas 

MY 

32 
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BEAUFORT  SUMMER  1982 


DATE 

time 

PM1 

PA1 

A1 

FI 

PM2 

PA2 

a: 

F2 

UEL 

03 

OCT 

1982 

17 

39 

45 

193 

39 

8 

33 

1 83 

33 

e 

33 

0.0 

03 

OCT 

1982 

20 

44 

43 

195 

72 

6 

4= 

153 

81 

5 

51 

0.0 

03 

OCT 

1982 

17 

31 

54 

1 95 

132 

3 

42 

195 

133 

3 

42 

0.0 

03 

OCT 

1982 

17 

32 

24 

195 

132 

3 

42 

195 

132 

3 

42 

0.0 

03 

OCT 

1982 

21 

15 

44 

203 

29 

33 

too 

203 

29 

33 

100 

0.0 

02 

OCT 

1982 

22 

c: 

35 

208 

28 

45 

132 

174 

33 

45 

159 

0.0 

03 

OCT 

1982 

17 

30 

54 

213 

44 

10 

45 

21 1 

44 

10 

45 

0.0 

02 

OCT 

1982 

21 

37 

51 

225 

55 

12 

02 

1 55 

28 

28 

82 

0.0 

02 

OCT 

19e2 

21 

52 

03 

235 

40 

1 1 

45 

150 

38 

17 

63 

0.0 

03 

OCT 

1982 

20 

47 

29 

242 

242 

1 

25 

242 

55 

'5 

•?  s 

0.0 

02 

OCT 

1982 

19 

30 

45 

243 

61 

6 

38 

238 

53 

7 

39 

0.0 

03 

OCT 

1982 

20 

20 

15 

24S 

31 

1  1 

35 

237 

116 

4 

48 

0.0 

02 

OCT 

1982 

19 

21 

04 

250 

127 

4 

53 

153 

50 

18 

94 

0.0 

03 

OCT 

1982 

22 

24 

33 

250 

51 

13 

70 

250 

51 

13 

70 

0.0 

03 

OCT 

1932 

21 

IS 

50 

253 

85 

S 

43 

253 

85 

5 

45 

0.0 

01 

OCT 

1982 

1 1 

43 

07 

250 

119 

4 

50 

250 

119 

4 

50 

0.0 

03 

OCT 

1982 

20 

29 

50 

250 

162 

2 

38 

245 

205 

2 

43 

0.0 

03 

OCT 

1982 

20 

45 

42 

255 

83 

1 1 

95 

255 

53 

11 

95 

0.0 

03 

OCT 

1982 

17 

43 

53 

275 

100 

5 

52 

252 

35 

15 

57 

0.0 

03 

OCT 

1982 

22 

27 

24 

292 

292 

1 

31 

292 

292 

1 

31 

0.0 

01 

OCT 

1982 

12 

40 

45 

297 

63 

5 

33 

278 

157 

S 

82 

0.0 

07 

OCT 

1982 

17 

52 

51 

298 

74 

1 1 

85 

269 

72 

13 

98 

0.0 

03 

OCT 

1982 

20 

43 

23 

311 

78 

5 

41 

1 55 

93 

6 

62 

0.0 

07 

OCT 

1982 

17 

53 

57 

316 

95 

10 

101 

291 

93 

1 1 

113 

0.0 

04 

OCT 

1932 

00 

05 

53 

324 

103 

S 

54 

291 

85 

6 

54 

0.0 

01 

OCT 

1982 

12 

37 

52 

325 

139 

4 

53 

325 

139 

4 

53 

0.0 

02 

OCT 

1982 

19 

35 

34 

327 

129 

4 

54 

327 

129 

4 

54 

0.0 

02 

OCT 

1932 

19 

50 

48 

327 

245 

4 

103 

325 

45 

31 

150 

0.0 

01 

OCT 

1982 

12 

39 

38 

332 

198 

2 

42 

137 

50 

8 

50 

0.0 

07 

OCT 

1982 

17 

55 

27 

333 

104 

4 

44 

237 

65 

8 

55 

0.0 

03 

OCT 

1982 

17 

30 

24 

334 

82 

8 

63 

334 

82 

8 

69 

0.0 

03 

OCT 

1932 

20 

45 

12 

350 

ua 

4 

52 

550 

U8 

4 

62 

0.0 

03 

OCT 

1982 

18 

03 

34 

355 

90 

6 

57 

348 

81 

8 

53 

0.0 

07 

OCT 

1932 

18 

07 

S3 

353 

323 

2 

63 

353 

323 

2 

63 

0.0 

02. 

OCT 

1982 

20 

05 

45 

352 

193 

3 

61 

255 

25 

31 

81 

0.0 

03 

OCT 

1982 

18 

03 

04 

364 

55 

10 

301 

49 

21 

108 

0.0 

03 

OCT 

1962 

21 

1 1 

50 

354 

354 

1 

38 

279 

104 

6 

65 

0.0 

01 

OCT 

1382 

21 

17 

S3 

355 

22 

54 

125 

312 

23 

S3 

129 

0.0 

03 

OCT 

1982 

28 

47 

370 

113 

12 

1  -46 

370 

1  IS 

12 

US 

0.0 

02 

OCT 

1932 

19 

44 

01 

371 

43 

20 

101 

320 

55 

19 

112 

0.0 

02 

OCT 

1982 

20 

12 

32 

374 

299 

2 

63 

291 

•>  -> 

*- 

35 

81 

0.0 

01 

OCT 

1  ao? 

*  WW4. 

1  1 

SO 

05 

375 

278 

2 

53 

355 

293 

2 

61 

0.0 

01 

OCT 

1932 

12 

22 

12 

378 

76 

9 

72 

204 

50 

19 

100 

0.0 

07 

OCT 

1982 

17 

52 

21 

373 

63 

10 

55 

293 

81 

10 

85 

0.0 

01 

OCT 

1  S32 

12 

25 

34 

390 

34 

14 

50 

390 

34 

14 

50 

0.0 

01 

CCT 

1982 

21 

23 

29 

3=2 

50 

43 

226 

377 

S3 

42 

234 

0.0 

02 

OCT 

1  982 

19 

53 

03 

395 

395 

1 

41 

315 

13 

30 

57 

0.0 

01 

OCT 

1932 

12 

20 

25 

404 

72 

7 

S3 

274 

97 

10 

102 

0.0 

07 

OCT 

1932 

18 

55 

17 

405 

71 

1  I 

82 

405 

71 

1 1 

32 

0.0 

07 

OCT 

1982 

13 

55 

419 

102 

7 

75 

371 

I  OS 

7 

78 

0.0 
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BEAUFORT  SUMMER  1982 


OATS 

TIME 

PM1 

PA1 

A! 

FI 

PM2 

PA2 

A2 

F  2 

OEL 

01 

OCT 

1932 

12 

39 

02 

420 

125 

6 

79 

420 

125 

6 

79 

0.0 

02 

OCT 

1952 

18 

33 

34 

421 

64 

10 

67 

408 

80 

8 

67 

0.0 

05 

OCT 

1952 

19 

45 

31 

425 

68 

1 3 

93 

326 

59 

15 

S3 

0.0 

03 

OCT 

1902 

22 

25 

18 

425 

135 

10 

142 

421 

1 36 

10 

143 

0.0 

03 

OCT 

1932 

21 

13 

S3 

433 

103 

5 

57 

198 

52 

13 

71 

0.0 

U 

OCT 

1982 

1 1 

59 

00 

434 

118 

9 

1 1 1 

419 

102 

1  1 

118 

0.0 

07 

OCT 

1982 

18 

35 

32 

436 

50 

17 

89 

389 

Et 

15 

SS 

0.0 

07 

OCT 

1932 

18 

29 

22 

437 

77 

1 1 

89 

437 

77 

1  1 

es 

0.0 

02 

OCT 

1  982 

20 

30 

09 

439 

41 

43 

185 

271 

42 

45 

203 

0.0 

02 

OCT 

1982 

20 

25 

04 

440 

76 

1 1 

83 

339 

45 

22 

104 

0.0 

02 

OCT 

1982 

20 

35 

1  1 

445 

63 

25 

165 

445 

63 

25 

165 

0.0 

07 

OCT 

1932 

18 

29 

SI 

448 

226 

4 

95 

443 

226 

4 

95 

0.0 

13 

OCT 

1932 

19 

18 

Ci. 

452 

142 

7 

104 

389 

144 

9 

1 35 

0.0 

07 

OCT 

1982 

18 

07 

24 

452 

120 

14 

176 

452 

120 

14 

175 

0.0 

02 

OCT 

1932 

18 

05 

20 

453 

27 

31 

88 

335 

115 

1 1 

1 33 

0.0 

07 

OCT 

1982 

18 

19 

14 

455 

48 

16 

81 

280 

71 

15 

112 

0.0 

01 

OCT 

1982 

1 1 

39 

47 

465 

194 

4 

81 

397 

121 

7 

89 

0.0 

01 

OCT 

1982 

12 

22 

39 

472 

105 

20 

222 

407 

102 

21 

225 

0.0 

03 

OCT 

1982 

18 

02 

35 

475 

91 

9 

86 

431 

67 

19 

134 

0.0 

07 

OCT 

1982 

17 

55 

57 

475 

63 

9 

59 

288 

1 55 

4 

65 

0.0 

14 

OCT 

1932 

08 

17 

38 

481 

43 

48 

217 

481 

43 

43 

217 

0.0 

10 

OCT 

1982- 

23 

14 

33 

482 

100 

16 

163 

437 

102 

16 

171 

0.0 

12 

OCT 

1982 

19 

S3 

00 

483 

92 

9 

87 

374 

1  1 1 

13 

151 

0.0 

03 

OCT 

1982 

19 

42 

32 

494 

69 

14 

101 

494 

69 

U 

101 

0.0 

13 

OCT 

1982 

20 

15 

1 1 

505 

368 

3 

116 

505 

368 

3 

,116 

0.0 

12 

OCT 

1982 

19 

09 

51 

510 

100 

7 

73 

510 

100 

7 

73 

0.0 

07 

OCT 

1982 

18 

13 

34 

51 1 

41 

22 

95 

469 

38 

24 

96 

0.0 

10 

OCT 

1  Qon 

1  JU4. 

18 

24 

24 

514 

190 

5 

100 

542 

90 

12 

113 

0.0 

14 

OCT 

1982 

1  I 

33 

42 

514 

86 

9 

81 

514 

86 

9 

81 

0.0 

12 

OCT 

1982 

18 

41 

49 

516 

1  14 

10 

120 

516 

114 

10 

129 

0.0 

14 

OCT 

1982 

1  1 

46 

05 

524 

46 

19 

92 

522 

46 

19 

92 

0.0 

03 

OCT 

1982 

20 

31 

59 

525 

93 

1  1 

107 

436 

140 

1 1 

162 

0.9 

12 

OCT 

1982 

16 

46 

28 

525 

74 

20 

155 

248 

101 

15 

159 

0.0 

10 

OCT 

1982 

15 

S3 

18 

535 

91 

8 

76 

360 

63 

12 

79 

0.0 

.01 

OCT 

1982 

12 

51 

37 

542 

542 

1 

57 

542 

542 

1 

57 

0.0 

1 1 

OCT 

1982 

00 

19 

37 

544 

1  94 

9 

183 

475 

242 

8 

203 

0.0 

03 

OCT 

1982 

20 

29 

06 

546 

171 

5 

90 

402 

87 

10 

SI- 

0.0 

07 

OCT 

1932 

18 

23 

52 

546 

124 

8 

104 

239 

75 

17 

134 

0.0 

03 

OCT 

1932 

21 

12 

59 

550 

121 

5 

63 

550 

121 

5 

53 

0.0 

03 

OCT 

1932 

20 

47 

00 

551 

77 

34 

275 

550 

83 

33 

237 

0.0 

02 

OCT 

1982 

19 

00 

25 

553 

139 

8 

117 

484 

127 

10 

1 33 

0.0 

03 

OCT 

1982 

21 

10 

32 

557 

93 

13 

127 

397 

33 

14 

131 

0.0 

12 

OCT 

1932 

IS 

53 

L.L. 

550 

124 

10 

130 

3 1 6 

1  IS 

1  1 

1 33 

0.0 

12 

OCT 

1332 

16 

41 

20 

553 

94 

14 

138 

444 

122 

12 

154 

0.0 

10 

OCT 

1982 

13 

43 

51 

559 

182 

8 

153 

299 

125 

16 

210 

0.0 

12 

OCT 

1992 

17 

16 

31 

573 

96 

17 

171 

500 

94 

21 

207 

0.0 

14 

OCT 

1982 

1  1 

26 

4S 

531 

257 

3 

81 

514 

59 

U 

87 

0.0 

14 

OCT 

1  982 

13 

30 

39 

532 

134 

7 

135 

505 

213 

7 

1 60 

0.0 

12 

OCT 

1982 

IS 

57 

53 

533 

95 

9 

90 

329 

61 

23 

179 

0.0 

U 

OCT 

1932 

1  4 

01 

53 

533 

94 

13 

123 

373 

1  1  9 

15 

137 

0.0 

A-4 


BEAUFORT  SUMMER  1982 


DATE 

TIME 

PM1 

PA1 

A1 

FI 

PM2 

Pn2 

A2 

F2 

VEL 

12 

OCT 

1882 

18 

10 

54 

534 

180 

6 

101 

504 

161 

7 

lie 

0.0 

12 

OCT 

1  8S2 

20 

05 

43 

584 

58 

23 

140 

552 

S3 

38 

21  1 

0.0 

13 

OCT 

I  882 

IS 

17 

49 

557 

1 55 

27 

385 

514 

154 

24 

388 

0.0 

10 

OCT 

1832 

1  5 

37 

38 

600 

55 

14 

97 

432 

80 

15 

126 

0.0 

10 

OCT 

1882 

15 

55 

31 

607 

75 

22 

173 

607 

75 

22 

173 

0.0 

11 

OCT 

1882 

00 

34 

610 

1  Eo 

7 

115 

574 

165 

8 

155 

0.0 

13 

OCT 

1882 

19 

53 

57 

618 

72 

1 1 

S3 

610 

72 

1  1 

83 

0.0 

10 

OCT 

1  S32 

16 

48 

14 

630 

130 

13 

177 

630 

130 

13 

1 77 

0.0 

13 

OCT 

1882 

17 

20 

07 

631 

1  18 

1  1 

134 

480 

94 

14 

138 

0.0 

07 

OCT 

1932 

19 

39 

55 

633 

181 

6 

114 

653 

181 

6 

1  14 

0.0 

14 

OCT 

1  982 

13 

31 

09 

654 

C“? 

W  A* 

14 

1 55 

612 

65 

15 

135 

0.0 

10 

OCT 

1832 

15 

SS 

01 

535 

129 

7 

95 

475 

142 

7 

104 

0.0 

0t 

OCT 

1932 

12 

37 

25 

638 

1  10 

7 

81 

612 

134 

6 

84 

0.0 

14 

OCT 

1  832 

08 

31 

00 

633 

27 

50 

142 

541 

63 

48 

317 

0.0 

13 

OCT 

1  882 

20 

10 

27 

642 

104 

18 

196 

342 

114 

17 

203 

0.0 

12 

OCT 

1882 

14 

51 

47 

645 

288 

3 

91 

503 

67 

21 

U3 

0.0 

08 

OCT 

1982 

10 

27 

50 

632 

£5 

21 

167 

632 

es 

21 

187 

0.0 

10 

OCT 

isez 

17 

05 

13 

655 

79 

19 

157 

GEE 

79 

19 

1 57 

0.0 

02 

OCT 

1  9S2 

20 

43 

21 

£55 

83 

35 

505 

451 

87 

35 

319 

0.0 

10 

OCT 

1982 

16 

38 

45 

653 

183 

7 

145 

653 

198 

7 

145 

0.0 

01 

OCT 

1982 

12 

23 

07 

653 

153 

12 

193 

650 

138 

16 

232 

0.0 

12 

OCT 

1982 

14 

39 

52 

678 

120 

12 

151 

655 

116 

13 

153 

0.0 

14 

OCT 

1982 

13 

45 

17 

634 

53 

16 

89 

684 

53 

15 

89 

0.0 

10 

OCT 

1982 

17 

05 

42 

637 

153 

6 

99 

637 

153 

6 

99 

0.0 

1 1 

OCT 

1982 

00 

12 

22 

'  687 

24S 

6 

154 

682 

218 

7 

150 

0.0 

03 

OCT 

1962 

09 

00 

44 

685 

472 

"> 

QQ 

w 

333 

49 

24 

123 

0.0 

03 

OCT 

1982 

20 

48 

16 

702 

138 

16 

232 

627 

112 

20 

235 

0.0 

02 

OCT 

1982 

17 

51 

42 

714 

115 

8 

97 

714 

115 

8 

97 

0.0 

07 

OCT 

1982 

20 

43 

46 

718 

155 

8 

130 

718 

155 

8 

130 

0.0 

10 

OCT 

1982 

23 

20 

04 

725 

124 

12 

1  Eo 

699 

1  14 

14 

157 

0.0 

14 

OCT 

1932 

13 

48 

05 

727 

1  1  4 

9 

108 

595 

U2 

9 

134 

0.0 

01 

OCT 

1982 

12 

0S 

27 

733 

94 

14 

138 

689 

123 

1 1 

143 

0.0 

10 

OCT 

1982 

23 

19 

30 

759 

438 

c 

92 

335 

122 

10 

128 

0.0 

03 

OCT 

1982 

07 

43 

40 

763 

172 

1 1 

193 

537 

79 

35 

290 

0.0 

12 

OCT 

1582 

14 

37 

49 

767 

72 

21 

153 

680 

72 

22 

166 

0.0 

14 

OCT 

1982 

12 

48 

54 

782 

95 

28 

279 

712 

101 

27 

266 

0.0 

07 

OCT 

1  932 

1  9 

4? 

23 

794 

1  96 

12 

247 

749 

1 03 

2S 

270 

0.0 

07 

OCT 

1  832 

23 

02 

29 

803 

125 

1 1 

144 

451 

1 55 

\  3 

21 1 

0.0 

03 

OCT 

1392 

07 

31 

42 

817 

353 

4 

150 

622 

163 

14 

247 

0.0 

10 

OCT 

1582 

16 

54 

25 

818 

243 

5 

130 

670 

335 

4 

162 

0.0 

H 

OCT 

1  232 

00 

21 

52 

819 

213 

5 

112 

792 

221 

5 

1  IS 

0.0 

03 

OCT 

1882 

00 

35 

10 

873 

175 

8 

147 

757 

188 

8 

153 

0.0 

08 

OCT 

1962 

09 

54 

07 

877 

145 

19 

289 

781 

81 

36 

30S 

0.0 

10 

OCT 

1592 

17 

58 

27 

878 

124 

10 

130 

878 

124 

10 

130 

0.0 

12 

OCT 

1532 

16 

2 4 

33 

889 

ISO 

16 

252 

889 

150 

l 5 

252 

0.0 

10 

OCT 

1  222 

15 

23 

02 

SOI 

45 

32 

151 

295 

57 

4  *"» 

**  fm 

251 

0.0 

07 

OCT 

1  ap? 

17 

c  a 

31 

920 

223 

8 

187 

539 

100 

19 

139 

0.0 

07 

OCT 

1  2S2 

18 

2S 

20 

523 

67 

29 

204 

523 

67 

29 

204 

0.0 

01 

OCT 

1  582 

12 

20 

53 

551 

230 

7 

335 

S3 

19 

195 

0.0 

13 

OCT 

I  rfCi 

1  3 

57 

54 

5=3 

1  53 

9 

154 

947 

133 

I  1 

153 

0.0 

A-5 


BEAUFORT  SUMMER  1982 


OATS 

TIME 

pm 

PA  1 

A1 

FI 

PM2 

PA2 

A2 

F2 

WEL 

10 

OCT 

1382 

15:44 

41 

380 

330 

8 

277 

608 

183 

20 

384 

0.0 

10 

OCT 

1382 

16:35 

14 

1010 

172 

13 

235 

381 

173 

14 

254 

0.0 

02 

OCT 

1532 

20:10 

03 

1013 

157 

12 

138 

306 

146 

14 

214 

0.0 

12 

OCT 

1382 

18:58 

17 

1015 

254 

11 

233 

862 

200 

17 

357 

0.0 

07 

OCT 

1382 

(8:00 

00 

1023 

2?? 

7 

153 

878 

210 

8 

176 

0.0 

07 

OCT 

1382 

18:48 

1 1 

1030 

174 

8 

145 

7SS 

138 

3 

187 

0.0 

12 

OCT 

1382 

17:07 

44 

1053 

518 

3 

483 

1053 

518 

3 

483 

0.0 

10 

OCT 

1382 

13:20 

07 

1033 

308 

7 

225 

1033 

308 

7 

225 

0.0 

10 

OCT 

I3S2 

18:41 

15 

1103 

115 

16 

133 

1103 

115 

16 

133 

0.0 

13 

OCT 

1382 

20:14 

41 

1115 

135 

10 

206 

775 

201 

11 

232 

0.0 

88 

NO  V 

1382 

88:88 

88 

1140 

306 

5 

161 

317 

337 

4 

167 

0.0 

14 

OCT 

(382 

11:48 

28 

1155 

212 

11 

245 

752 

207 

20 

434 

0.0 

14 

OCT 

1382 

11:30 

12 

1206 

136 

15 

214 

1205 

135 

15 

214 

0.0 

01 

OCT 

1382 

12:21 

26 

1453 

183 

14 

263 

1 365 

175 

15 

275 

0.0 

10 

OCT 

1582 

18:38 

IS 

1464 

335 

12 

422 

1365 

344 

12 

433 

0.0 

07 

OCT 

1382 

23:30 

23 

1433 

334 

10 

413 

1433 

334 

10 

413 

0.0 

14 

OCT 

1332 

1 1 :37 

33 

1617 

235 

IS 

435 

1617 

235 

16 

435 

0.0 

A-6 


SOUTH  BERING  WINTER  1983 


OATE 

TIME 

PM1  1 

PA1 

26 

MAR 

1283 

11:57:53 

26 

10 

26 

MAR 

1283 

12:2:17 

36 

36 

26 

MAR 

1  2S3 

1  1  : 55: 4-7 

37 

16 

26 

MAR 

1283 

1 1:57:46 

39 

6 

26 

MAR 

1283 

12:4:28 

40 

4 

26 

MAR 

1233 

12:7:42 

42 

16 

26 

MAR 

1283 

1 1 :53:35 

43 

30 

26 

MAR 

1283 

1 1 :S4: 1 

47 

6 

26 

MAR 

1283 

12:2:25 

43 

10 

26 

MAR 

1983 

1 1 :S4:20 

50 

12 

26 

MAR 

1983 

11:54:12 

51 

27 

26 

MAR 

1983 

12:14:32 

52 

8 

26 

MAR 

1983 

1 1 :55:1 

54 

8 

26 

MAR 

1983 

1 1 :54:30 

62 

7 

26 

MAR 

1983 

1 1 :53:5o 

68 

13 

26 

MAR 

1983 

19:31:44 

113 

21 

26 

MAR 

1983 

19:51 :23 

117 

18 

26 

MAR 

1983 

12:51 :7 

121 

21 

26 

MAR 

1983 

19:46:13 

123 

25 

26 

MAR 

1983 

19:0:10 

135 

26 

MAR 

1283 

20:32:26 

135 

25 

26 

MAR 

1983 

18:59:43 

137 

49 

26 

MAR 

198.5 

: S : 4S: 1 4 

137 

17 

26 

MAR 

1983 

19:45:52 

139 

19 

25 

MAR 

1283 

19:9:27 

140 

IS 

26 

MAR 

1283 

18:55:39 

141 

74 

26 

MAR 

1983 

19:51:17 

141 

18 

26 

MAR 

ises 

20:32:43 

143 

19 

26 

MAR 

1983 

18:46:6 

144 

17 

25 

MAR 

1983 

19:45:5 

145 

17 

26 

MAR 

1283 

19: 10: 12 

148 

27 

26 

MAR 

1233 

19:40:1 

149 

18 

26 

MAR 

1983 

19:51:1 

149 

16 

26 

MAR 

1983 

10:11:17 

152 

20 

25 

MAR 

1933 

19:50:24 

152 

22 

26 

MAR 

1933 

19:45:27 

153 

17 

26 

MAR 

1983 

18:39:23 

154 

44 

26 

MAR 

1983 

18:47:53 

154 

30 

26 

MAR 

1983 

18:46:39 

154 

20 

26 

MAR 

1  933 

19:50:55 

154 

UU 

26 

MAR 

1283 

18:51 :29 

154 

53 

26 

MAR 

1983 

1 :42:29 

157 

18 

25 

MAR 

1283 

18:59:57 

157 

19 

26 

MAR 

1983 

19:10:1 

157 

27 

26 

MAR 

1283 

o :  j  *  *:4 

160 

17 

26 

MAR 

1983 

2:55:43 

161 

20 

26 

MAR 

1283 

18:45:51 

161 

IS 

26 

MAR 

1283 

19:46:7 

152 

15 

26 

MAR 

1283 

19:51 :42 

153 

20 

26 

MAR 

1983 

18:SS:2S 

154 

47 

Al 

FI 

PM2  1 

=A2 

A2 

F2 

UEL 

IS 

15 

5 

45 

22 

0.0 

t 

14 

20 

4 

40 

18 

0.0 

3 

IS 

21 

S 

43 

21 

0.0 

42 

29 

30 

7 

41 

31 

0.0 

35 

14 

26 

9 

15 

25 

0.0 

7 

25 

37 

9 

15 

31 

C.0 

3 

26 

45 

9 

38 

36 

0.0 

38 

27 

44 

8 

33 

29 

0.0 

8 

22 

49 

10 

8 

22 

0.0 

9 

17 

35 

7 

38 

29 

0.0 

2 

IS 

28 

5 

48 

23 

0.0 

18 

21 

44 

7 

33 

24 

0.0 

47 

38 

42 

9 

46 

44 

0.0 

41 

32 

43 

8 

49 

43 

0.0 

8 

17 

47 

1 1 

9 

19 

0.0 

21 

59 

90 

17 

33 

69 

0.0 

48 

98 

92 

20 

48 

101 

0.0 

45 

101 

105 

24 

47 

118 

0.0 

45 

118 

123 

25 

45 

118 

0.0 

48 

1  I  1 

32 

23 

49 

1  17 

0.0 

47 

125 

89 

26 

50 

158 

0.0 

4 

35 

125 

59 

3 

43 

0.0, 

S3 

93 

92 

18 

50 

94 

0.0 

46 

90 

86 

21 

43 

96 

0.0 

22 

42 

89 

18 

10 

48 

0.0 

5 

57 

141 

60 

10 

81 

0.0 

45 

87 

89 

31 

12 

93 

0.0 

52 

104 

143 

19 

52 

104 

0.0 

12 

29 

82 

14 

21 

36 

0.0 

44 

79 

92 

19 

48 

95 

0.0 

10 

43 

77 

19 

27 

64 

0.0 

41 

78 

102 

18 

46 

86 

0.0 

45 

74 

100 

22 

46 

106 

0.0 

20 

49 

99 

21 

20 

52 

0.0 

45 

105 

92 

24 

45 

114 

0.0 

5! 

90 

92 

17 

51 

90 

0.0 

8 

59 

154 

44 

8 

59 

0.0 

43 

139 

94 

35 

44 

166 

0.0 

39 

93 

154 

20 

39 

93 

0.0 

44 

104 

154 

22 

44 

104 

0.0 

7 

127 

154 

53 

7 

127 

0.0 

59 

114 

109 

22 

55 

129 

0.0 

27 

53 

1  15 

18 

30 

61 

0.0 

18 

59 

1  13 

23 

29 

91 

0.0 

56 

98 

65 

IS 

59 

1  10 

0.0 

53 

120 

S2 

24 

59 

145 

0.0 

12 

50 

161 

18 

12 

50 

0.0 

48 

77 

92 

18 

47 

92 

0.0 

49 

(01 

163 

20 

49 

101 

0.0 

4 

37 

124 

30 

10 

43 

0.0 

A-  7 


SOUTH  BERING  WINTER  1983 


DATE 

TIME 

PM1 

PA1 

A 1 

FI 

PM2  1 

PA2 

A2 

F2 

UEL 

26 

MAR 

1983 

19:15:25 

164 

164 

l 

38 

88 

19 

10 

42 

0.0 

26 

MAR 

1983 

19:52:13 

164 

18 

46 

87 

162 

35 

13 

88 

0.0 

25 

MAR 

1983 

10:17:33 

165 

33 

6 

35 

152 

99 

L 

36 

0.0 

25 

MAR 

1983 

2j  :  6  •  o 

1 65 

IS 

51 

82 

123 

•n 

cu 

50 

1  17 

0.0 

26 

MAR 

1983 

10:18:5 

166 

14 

42 

60 

146 

23 

30 

75 

0.0 

25 

MAR 

1983 

17:12:0 

167 

35 

12 

52 

162 

39 

1 1 

55 

0.0 

25 

MAR 

1983 

23:29:14 

167 

IS 

5! 

78 

167 

15 

51 

78 

0.0 

26 

MAR 

1983 

19:39:47 

169 

15 

46 

72 

103 

21 

45 

100 

0.0 

25 

MAR 

1983 

18:35:28 

170 

23 

16 

44 

90 

39 

8 

50 

0.0 

26 

MAR 

1983 

10:18:17 

172 

24 

39 

98 

172 

24 

39 

98 

0.0 

26 

MAR 

1983 

1 9:45:45 

172 

18 

48 

92 

1  12 

22 

49 

115 

0.0 

26 

MAR 

1983 

3:20:9 

174 

19 

59 

120 

174 

19 

59 

120 

0.0 

26 

MAR 

1983 

19:51 :36 

174 

18 

45 

87 

93 

26 

43 

117 

0.0 

26 

MAR 

1983 

18:55:45 

176 

60 

3 

42 

76 

22 

14 

50 

0.0 

26 

MAR 

1983 

19:49:33 

176 

25 

46 

121 

168 

26 

45 

122 

0.0 

26 

MAR 

1983 

0:44:7 

177 

16 

56 

93 

138 

16 

55 

96 

0.0 

26 

MAR 

1983 

0:48:47 

177 

17 

56 

102 

130 

24 

53 

134 

0.0 

25 

MAR 

1983 

23:4:59 

179 

18 

53 

102 

148 

27 

49 

138 

0.0 

26 

MAR 

1983 

10:18:25 

179 

19 

38 

81 

1 55 

23 

37 

89 

0.0 

26 

MAR 

1983 

2:53:3 

180 

32 

54 

184 

175 

32 

55 

1 85 

0.0 

25 

MAR 

1983 

22:S7:37 

181 

22 

52 

119 

155 

23 

52 

127 

0.0 

25 

MAR 

1983 

2:58:38 

182 

19 

53 

115 

182 

19 

58 

IIS 

0.0 

25 

MAR 

1983 

22:17:58 

183 

19 

52 

101 

146 

32 

46 

155 

0.0 

ZS 

MAR 

1983 

18:55:29 

183 

S3 

7 

48 

147 

82 

5 

55 

0.0 

26 

MAR 

1983 

0:54:53 

184 

33 

50 

172 

184 

33 

50 

172 

0.0 

26 

MAR 

1983 

10: 1 : 12 

184 

39 

7 

45 

121 

15 

13 

47 

0.0 

26 

MAR 

1983 

10:4:48 

184 

32 

20 

77 

142 

24 

32 

86 

0.0 

25 

MAR 

1983 

21:40:48 

185 

13 

54 

73 

161 

lbl 

1 

130 

0.0 

25 

MAR 

1983 

22:57:18 

187 

21 

52 

1  17 

187 

21 

52 

117 

0.0 

26 

MAR 

1983 

2:44:25 

187 

27 

53 

165 

175 

30 

57 

177 

0.0 

25 

MAR 

1983 

10:11:37 

187 

29 

1  1 

48 

84 

15 

36 

62 

0.0 

26 

MAR 

1983 

19:49:39 

188 

24 

46 

118 

172 

25 

49 

127 

0.0 

25 

MAR 

1993 

23:5:48 

191 

19 

55 

107 

1 65 

32 

49 

162 

0.0 

25 

MAR 

1983 

23:3: 27 

193 

20 

55 

115 

143 

29 

48 

146 

0.0 

26 

MAR 

1983 

19:49:27 

193 

1  9 

31 

89 

155 

31 

26 

1 1  1 

0.0 

26 

MAR 

1983 

2:54:5 

194 

17 

57 

100 

1  36 

19 

55 

1  10 

0.0 

ZS 

MAR 

1983 

22:29:22 

197 

12 

57 

63 

197 

12 

57 

69 

0.0 

26 

MAR 

1983 

0:54:54 

197 

18 

59 

112 

ISO 

22 

55 

124 

0.0 

25 

MAR 

1983 

22:11:18 

198 

21 

49 

1  10 

135 

28 

40 

128 

0.0 

26 

MAR 

1983 

0:57:44 

198 

19 

55 

1  10 

182 

21 

54 

1  17 

0.0 

26 

MAR 

1983 

1  9:45:41 

198 

21 

41 

108 

117 

22 

5 1 

121 

0.0 

25 

MAR 

1983 

16:54:6 

199 

40 

12 

67 

128 

26 

37 

106 

0.0 

26 

MAR 

1983 

0:54:25 

200 

33 

51 

177 

193 

34 

50 

179 

0.0 

25 

MAR 

1983 

23:11:57 

201 

25 

55 

144 

150 

29 

50 

154 

0.0 

25 

MAR 

1983 

22:29:7 

202 

21 

S3 

IIS 

162 

21 

54 

1  16 

0.0 

26 

MAR 

1983 

1 :8 : 43 

202 

24 

53 

146 

189 

29 

55 

165 

0.0 

25 

MAR 

1983 

2:49:27 

203 

28 

53 

172 

177 

32 

57 

194 

0.0 

26 

MAR 

1983 

2:54:49 

203 

18 

57 

1  10 

202 

19 

57 

112 

0.0 

26 

MAR 

1983 

2:43:55 

204 

17 

57 

103 

95 

19 

55 

1  1  3 

0.0 

25 

MAR 

1983 

23:3:51 

205 

19 

52 

106 

197 

25 

52 

136 

0.0 

A- 3 


SOUTH  BERING  WINTER  1983 


DATE 

TIME 

PM1 

PA1 

A1 

F! 

PMZ 

PA2 

A2 

F2 

VEL 

25 

MAR 

1583 

22:2:42 

205 

18 

51 

95 

100 

18 

54 

103 

0.0 

25 

MAR 

1983 

21:49:12 

207 

50 

7 

55 

123 

30 

27 

97 

0.0 

25 

MAR 

1983 

2:50:59 

207 

22 

53 

1 36 

191 

25 

55 

145 

0.0 

25 

MAR 

1983 

19:49:12 

207 

17 

46 

83 

143 

22 

46 

107 

0.0 

25 

MAR 

1983 

17: 1 1 :S5 

209 

17 

35 

66 

115 

18 

36 

69 

0.0 

25 

MAR 

1983 

22:43:20 

210 

30 

45 

141 

210 

30 

45 

141 

0.0 

25 

MAR 

1983 

22:55:57 

21  1 

12 

55 

69 

1 35 

15 

52 

82 

0.0 

25 

MAR 

1983 

n.ct  «  •?*? 

4.  -  -  1  •  4-- 

21 1 

19 

53 

116 

21  1 

19 

53 

1  16 

0.0 

25 

MAR 

1933 

10:14:53 

21  1 

16 

33 

65 

157 

21 

29 

70 

0.0 

25 

MAR 

1983 

22:48:31 

214 

17 

49 

89 

99 

22 

45 

106 

0.0 

25 

MAR 

1983 

10:12:43 

214 

14 

44 

65 

201 

15 

43 

63 

0.0 

25 

MAR 

1983 

22:18:26 

216 

19 

53 

112 

154 

20 

52 

1  14 

0.0 

26 

MAR 

1983 

1 :5:42 

216 

25 

55 

146 

175 

25 

56 

149 

0.0 

26 

MAR 

1983 

10:8:45 

216 

IIS 

2 

33 

99 

52 

8 

53 

0.0 

25 

MAR 

1983 

17:17:33 

220 

55 

5 

47 

220 

56 

5 

47 

0.0 

26 

MAR 

1983 

1:16:2 

220 

18 

57 

1  10 

14S 

20 

60 

128 

0.0 

25 

MAR 

1983 

18:38:53 

220 

28 

42 

125 

186 

40 

31 

141 

0.0 

26 

MAR 

1  983 

13:55:24 

222 

34 

9 

45 

126 

30 

10 

45 

0.0 

25 

MAR 

1983  ' 

10:4:27 

223 

61 

4 

46 

106 

27 

17 

60 

0.0 

25 

MAR 

1983 

22:13:17 

226 

23 

49 

119 

216 

26 

48 

128 

0.0 

25 

MAR 

1983 

21 :S0: 1 

228 

27 

SI 

147 

228 

29 

50 

153 

0.0 

26 

MAR 

1983 

10:8:57 

228 

IS 

26 

45 

146 

87 

4 

68 

0.0 

25 

MAR 

1983 

2:52:25 

229 

18 

57 

108 

222 

>8 

55 

108 

0.0 

25 

MAR 

1983 

2:53:38 

231 

19 

57 

113 

178 

20 

55 

115 

0.0 

26 

MAR 

1933 

10:13:55 

251 

37 

21 

91 

231 

37 

21 

91 

0.0 

25 

MAR 

1933 

22:57:24 

234 

IS 

53 

85 

234 

15 

53 

85 

0.0 

26 

MAR 

1983 

10:12:17 

234 

34 

40 

145 

198 

37 

37 

14B 

0.0 

25 

MAR 

1983 

22:11 :30 

238 

26 

S0 

137 

1 55 

29 

51 

1 55 

0.0 

26 

MAR 

1983 

6:10:4 

241 

76 

5 

63 

17! 

25 

21 

70 

0.0 

25 

MAR 

1983 

23:11 :S1 

249 

38 

53 

211 

242 

39 

52 

212 

0.0 

25 

MAR 

1933 

17:40:59 

250 

52 

8 

69 

91 

15 

53 

81 

0.0 

26 

MAR 

1933 

2:43:31 

250 

49 

48 

247 

250 

49 

48 

247 

0.0 

26 

MAR 

1983 

5:31 : 41 

250 

44 

13 

65 

167 

40 

17 

76 

0.0 

25 

MAR 

1983 

23:5:25 

2S2 

24 

51 

130 

191 

24 

51 

1 30 

0.0 

26 

MAR 

1983 

5:5:24 

252 

33 

28 

107 

163 

37 

27 

116 

0.0 

25 

MAR 

1983 

22:29:51 

261 

32 

52 

173 

255 

31 

S3 

174 

0.0 

26 

MAR 

1933 

3:11:30 

261 

28 

53 

173 

261 

28 

53 

173 

0.0 

25 

MAR 

1983 

23:4:54 

262 

21 

54 

1  17 

174 

35 

43 

180 

0.0 

26 

MAR 

1983 

2:51  : 52 

263 

19 

57 

1  I  1 

253 

19 

57 

1  1  1 

0.0 

26 

MAR 

1983 

2:51 : 17 

268 

22 

57 

134 

1  33 

23 

S3 

141 

0.0 

26 

MAR 

1983 

2:49:0 

270 

20 

53 

121 

184 

22 

55 

126 

0.0 

26 

MAR 

1983 

19:59:22 

272 

23 

53 

126 

272 

23 

53 

125 

0.0 

26 

MAR 

1983 

19:11:29 

277 

24 

35 

89 

138 

24 

37 

96 

0.0 

25 

MAR 

1983 

22:10:51 

279 

17 

53 

95 

267 

16 

55 

97 

0.0 

25 

MAR 

1983 

19:45:21 

280 

20 

49 

103 

96 

19 

53 

105 

0.0 

25 

MAR 

1983 

20:42:30 

288 

22 

44 

103 

282 

31 

40 

130 

0.0 

25 

MAR 

1583 

22:28:53 

299 

16 

55 

90 

186 

21 

55 

121 

0.0 

25 

MAR 

1983 

3:7:50 

289 

24 

55 

1  w 

n  -»  « 

7  « 
w  i 

53 

t  og 

0 . 0 

25 

MAR 

1983 

17:41:32 

291 

27 

45 

129 

249 

32 

36 

131 

0.0 

25 

MAR 

1983 

2:46:24 

292 

20 

53 

123 

292 

20 

S3 

1  23 

0.0 

A-9 


SOUTH  BERING  WINTER  1983 


DATE 

TIME 

PM  1 

PA  1 

A1 

F! 

PM2 

PA2 

A2 

F2 

UEL 

2S 

MAR 

1983 

2:4i:8 

294 

22 

53 

1 35 

253 

24 

59 

147 

0.0 

25 

MAR 

1983 

2:45:54 

295 

24 

59 

146 

250 

29 

57 

171 

0.0 

26 

MAR 

1983 

2:43:36 

296 

43 

48 

217 

208 

39 

57 

253 

0.0 

26 

MAR 

1933 

3:8:9 

295 

26 

54 

146 

295 

25 

54 

146 

0.0 

26 

MAR 

1985 

19:44:39 

303 

20 

45 

97 

239 

28 

45 

133 

0.0 

26 

MAR 

1983 

10:17:25 

304 

27 

26 

82 

217 

28 

44 

133 

0.0 

26 

MAR 

1983 

0:43:39 

312 

24 

S3 

132 

125 

24 

55 

142 

0.0 

26 

MAR 

1983 

19:39:42 

327 

20 

41 

92 

319 

44 

22 

117 

0.0 

26 

MAR 

1983 

10:4:32 

335 

63 

20 

141 

160 

42 

29 

143 

0.0 

26 

MAR 

1983 

1 :55 : 8 

335 

31 

53 

188 

536 

31 

58 

138 

0.0 

26 

MAR 

1983 

19:10:40 

537 

37 

17 

94 

337 

37 

17 

94 

0.0 

26 

MAR 

1983 

6:18:39 

345 

35 

31 

121 

325 

33 

35 

141 

0.0 

26 

MAR 

1983 

2:55:53 

550 

25 

58 

153 

350 

25 

53 

153 

0.0 

26 

MAR 

1983 

19:51 : 12 

355 

32 

29 

113 

198 

25 

42 

1  15 

0.0 

25 

MAR 

1983 

5:7:53 

363 

21 

22 

50 

353 

21 

22 

50 

0.0 

26 

MAR 

1983 

.10:1 :38 

367 

87 

5 

53 

189 

120 

4 

59 

0.0 

25 

MAR 

1983 

22:2:30 

369 

19 

48 

95 

182 

18 

55 

103 

0.0 

26 

MAR 

1983 

6:9:51 

374 

37 

22 

93 

374 

37 

22 

93 

0.0 

26 

MAR 

1983 

3:19:43 

385 

24 

53 

144 

385 

24 

53 

144 

0.0 

25 

MAR 

1983 

23:29:5 

409 

28 

50 

148 

327 

29 

c-? 

w«» 

161 

0.0 

26 

MAR 

1983 

6:26:9 

431 

41 

23 

145 

431 

41 

23 

145 

0.0 

25 

MAR 

1983 

21  :49:3 

483 

42 

49 

216 

281 

45 

50 

236 

0.0 

25 

MAR 

1983 

2:52:14 

1137 

47 

51 

253 

717 

49 

54 

280 

0.0 

A- 10 


NORTH  BERING  WINTER  1983 


DATE  TIME  PM1  PA1  A!  FI  PM2  PA2  A2  F2  VEL 

27  MAR  1283  7:32:39  1U  19  S9  118  IH  19  59  118  0.0 

27  MAR  1983  7:32:22  124  29  7  G3  32  15  41  E9  0.0 

27  MAR  1983  7:34:44  152  20  53  119  83  20  58  121  0.0 

27  MAR  1963  7:34:S3  153  19  57  111  25  19  57  113  0.0 

27  MAR  1983  7:19:13  160  20  55  115  113  22  52  121  0.0 

27  MAR  1983  7:3S:S4  JS2  21  48  103  114  20  54  112  0.0 

27  MAR  1283  7:39:34  165  21  59  133  125  23  57  136  0.0 

27  MAR  1 98"  15:15:45  165  85  8  79  164  83  8  82  0.0 

27  MAR  19S3  7:23:34  163  20  58  119  132  21  55  121  0.0 

27  MAR  1983  6:59=50  170  19  35  81  137  114  3  93  0.0 

27  MAR  1583  6:S9*.36  175  16  54  88  173  16  S3  89  0.0 

27  MAR  1983  15:6:50  179  111  2  27  169  125  2  30  0.0 

27  MAR  1283  14:49:33  190  100  3  41  162  102  3  44  0.0 

27  MAR  1983  14:54:7  191  24  11  39  108  26  1 1  40  0.0 

27  MAR  1283  7:28:23  192  31  55  177  183  59  49  204  0.0 

27  MAR  1983  14:54:0  195  47  6  45  147  39  8  47  0.0 

26  MAR  1933  23:22:7  193  18  39  94  122  3!  45  155  0.0 

27  MAR  1983  1:3:33  200  20  57  120  193  39  46  197  0.0 

27  MAR  1283  15:13:54  200  81  7  63  200  81  7  63  0.0 

27  MAR  1283  7:14:58  203  21  57  128  203  21  57  128  0.0 

27  MAR  1983  7:3:7  204  16  50  86  102  17  49  25  0.0 

27  MAR  1283  6:24:41  207  20  59  126  153  25  55  149  0.0 

27  MAR  1983  14:35:46  209  28  41  121  209  28  41  121  0.0 

27  MAR  1 S83  1:24:37  211  23  55  131  191  23  53  133  0.0 

27  MAR  1933  7:3:15  216  24  55  136  135  26  52  143  0.0 

2 (  MAR  1983  0:19:38  218  20  48  100  218  20  48  100  0.0 

27  MAR  1983  0:46:28  220  21  54  119  186  25  51  133  0.0 

27  MAR  1983  0:1:25  222  23  49  123  U9  31  42  U3  0.0 

27  MAR  1983  22=51:34  229  78  7  69  229  78  7  69  0.0 

27  MAR  1983  0:19:32  231  41  49  209  231  41  49  209  0.0 

27  MAR  1983  7:24:5  232  27  56  157  232  27  56  157  0.0 

27  MAR  1983  14:47:52  236  43  9  63  136  86  5  71  0.0 

27  MAR  1933  21:9:48  236  54  5  54  164  41  10  70  0.0 

-27  MAR  1983  21 : S3 : 5 1  236  25  23  63  123  22  26  64  0.0 

27  MAR  1983  15:15:33  238  46  16  81  157  67  18  131  0.0 

27  MAR  1983  15=58:4  239  82  8  77  222  9S  8  90  0.0 

27  MAR  1983  15:16=50  245  101  10  113  196  105  10  117  0.0 

27  MAR  1983  21:10:9  245  46  16  83  237  50  19  102  0.0 

27  MAR  1983  0:55=38  246  21  54  121  184  49  46  238  0.0 

27  MAR  1985  1S:53=36  247  87  10  93  247  87  10  93  0.0 

27  MAR  1983  0:8:0  249  20  49  102  249  20  49  102  0.0 

27  MAR  1983  15=6:0  252  49  7  43  189  107  4  60  0.0 

27  MAR  1983  13  =  38:22  253  50  6  45  194  38  8  47  0.0 

28  MAR  1983  14:4:27  253  49  14  80  253  49  14  80  0.0 

27  MAR  1983  15:16  =  13  254  55  21  U0  254  55  21  MO  0.0 

27  MAR  1983  0:48:29  257  26  52  U1  204  42  49  214  0.0 

27  MAR  1983  15:14:40  259  65  1  I  63  259  66  1 1  83  0.0 

27  MAR  1983  21:45=29  262  114  4  91  201  90  9  127  0.0 

27  MAR  1983  0:48:4  255  25  52  134  265  25  52  134  3.0 

27  MAR  1983  13:52=38  265  65  8  53  151  41  16  74  0.0 

A- 11 


NORTH  BERING  WINTER  1983 


JATE 

TIME 

PM1 

PA! 

A1 

FI 

PM2 

PA2 

A2 

F2 

VEL 

28 

MAR 

1983 

0:13:41 

265 

24 

50 

130 

15= 

27 

46 

1 36 

0.0 

27 

MAR 

1983 

1S:7: IS 

266 

83 

6 

60 

265 

83 

6 

60 

0.0 

28 

MAR 

1953 

0:7:28 

256 

75 

8 

100 

264 

53 

1 3 

108 

0.0 

28 

MAR 

1983 

14:3:33 

267 

63 

17 

129 

267 

68 

17 

129 

0.0 

27 

MAR 

1983 

15:29:0 

263 

81 

18 

158 

22! 

70 

n 

165 

0.0 

27 

MAR 

1983 

7:23:29 

259 

42 

51 

222 

236 

42 

51 

227 

0.0 

27 

MAR 

1983 

22:27:42 

269 

31 

32 

107 

178 

36 

33 

125 

C.0 

28 

MAR 

1983 

0:32:35 

269 

25 

56 

U8 

134 

33 

54 

223 

0.0 

27 

MAR 

1983 

15:13:40 

271 

92 

5 

50 

271 

92 

5 

50 

0.0 

28 

MAR 

1  983 

21 : 1 : 43 

271 

67 

1 1 

80 

271 

67 

1 1 

80 

0.0 

27 

MAR 

1983 

0:14:55 

272 

15 

40 

63 

203 

27 

42 

120 

0.0 

27 

MAR 

1983 

0:48:12 

272 

25 

53 

143 

205 

31 

54 

175 

0.0 

27 

MAR 

1983 

18:22:56 

273 

44 

8 

43 

251 

62 

17 

123 

0.0 

28 

MAR 

1983 

22:17:0 

274 

157 

4 

61 

274 

157 

4 

81 

0.^ 

27 

MAR 

1983 

17:20:S4 

275 

65 

5 

42 

ie7 

46 

9 

43 

0.0 

27 

MAR 

1983 

15:29:15 

276 

212 

2 

49 

188 

127 

4 

62 

0.0 

27 

MAR 

1983 

22:10:8 

276 

55 

22 

157 

255 

69 

16 

170 

0.0 

28 

MAR 

1983 

0:6:39 

276 

125 

3 

76 

234 

82 

5 

79 

0.0 

28 

MAR 

1983 

21 :36:S6 

277 

88 

4 

42 

72 

20 

21 

55 

0.0 

28 

MAR 

1983 

0:20:27 

278 

29 

48 

153 

278 

29 

43 

153 

0.0 

27 

MAR 

1983 

15:12:26 

279 

106 

1  1 

137 

279 

106 

1 1 

137 

0.0 

27 

MAR 

1983 

17:30:21 

279 

108 

13 

152 

279 

108 

13 

152 

0.0 

27 

MAR 

1983 

23:22:33 

279 

18 

37 

84 

137 

22 

45 

109 

0.0 

28 

MAR 

1983 

13:50:7 

279 

30 

33 

117 

279 

30 

33 

117 

0.0 

28 

MAR 

1983 

0:12:26 

280 

24 

52 

131 

280 

24 

52 

131 

0.0 

27 

MAR 

1983 

0:55:22 

281 

25 

57 

147 

•VMr 

45 

51 

242 

0.0 

28 

MAR 

1983 

0:25:59 

281 

22 

S3 

125 

119 

25 

55 

145 

0.0 

27 

MAR 

1983 

IS: 13:23 

282 

57 

7 

44 

282 

57 

7 

44 

0.0 

27 

MAR 

1933 

21:59:30 

282 

77 

4 

50 

147 

26 

10 

57 

0.0 

27 

MAR 

1983 

17:43:16 

283 

151 

4 

71 

233 

151 

4 

71 

0.0 

28 

MAR 

1993 

0:34:8 

2e3 

26 

54 

145 

126 

28 

56 

162 

0.0 

27 

MAR 

1933 

1 :24:46 

284 

27 

SI 

142 

243 

30 

43 

153 

0.0 

27 

MAR 

1983 

16:55:49 

284 

40 

14 

71 

113 

36 

21 

84 

0.0 

28 

MAR 

1933 

1 1 :S0:S4 

284 

31 

25 

86 

284 

31 

25 

85 

0.0 

27 

MAR 

1993 

21 :51 :7 

285 

37 

16 

67 

281 

46 

17 

90 

0.0 

29 

MAR 

1933 

13:26:0 

285 

31 

41 

131 

235 

31 

4! 

131 

0.0 

27 

MAR 

1983 

15:12:50 

236 

106 

5 

64 

200 

51 

13 

78 

0.0 

28 

MAR 

1983 

7:45:55 

236 

21 

50 

108 

195 

23 

49 

121 

0.0 

28 

MAR 

1983 

7:55: 46 

287 

38 

52 

208 

287 

41 

49 

210 

0.0 

27 

MAR 

1983 

21:49:35 

283 

36 

24 

97 

288 

35 

24 

97 

0.0 

28 

MAR 

1983 

0:13:30 

289 

24 

44 

I  18 

239 

24 

44 

118 

0.0 

23 

MAR 

1983 

0:14:45 

291 

30 

43 

153 

269 

30 

49 

159 

0.0 

29 

MAR 

1933 

0:43:11 

291 

Jb 

55 

187 

208 

34 

54 

193 

0.0 

29 

MAR 

1983 

14:2:54 

291 

40 

1  3 

57 

126 

64 

10 

75 

0.0 

23 

MAR 

1  933 

0:27:17 

25 

54 

US 

208 

25 

53 

1 45 

0.0 

27 

MAR 

1983 

17:3:40 

293 

59 

9 

61 

1  86 

64 

6 

55 

0.0 

27 

MAR 

1933 

18:26:2 

293 

82 

4 

4J. 

145 

31 

18 

6! 

0.0 

27 

MAR 

1983 

13:37:29 

233 

71 

10 

73 

154 

75 

10 

85 

0.0 

29 

MAR 

i  583 

11:46:9 

294 

35 

20 

79 

132 

35 

25 

102 

0.0 

27 

M  .R 

1  383 

0:1:6 

298 

29 

51 

1  52 

205 

45 

157 

0.0 

A-12 


NORTH  BERING  WINTER  1983 


OATE 

TIME 

PM  1 

PA1 

At 

F! 

PM2 

PA2 

A2 

F2 

OEL 

27 

MAR 

1S63 

15:11 :58 

298 

51 

15 

84 

285 

55 

19 

1  1G 

0.0 

28 

MAR 

1  SS3 

0:39:30 

296 

25 

54 

143 

296 

25 

54 

143 

0.0 

27 

MAR 

1983 

1:3:43 

297 

25 

57 

149 

297 

25 

57 

1 49 

0.0 

27 

MAR 

1  383 

0:45=31 

298 

35 

47 

172 

230 

33 

49 

173 

0.0 

28 

MAR 

1983 

22:14:9 

300 

10S 

4 

82 

115 

27 

33 

92 

0.0 

27 

MAR 

1983 

15:12:8 

301 

77 

5 

49 

209 

SI 

8 

90 

0.0 

27 

MAR 

1983 

21  : 29 : 35 

301 

132 

3 

47 

250 

199 

2 

48 

0.0 

28 

MAR 

19e3 

21 :36:50 

30! 

128 

4 

61 

91 

26 

32 

91 

0.0 

27 

MAR 

1983 

7:28:34 

306 

32 

57 

191 

237 

44 

51 

241 

0.0 

27 

MAR 

1  S83 

15:8:5 

306 

119 

5 

40 

185 

97 

5 

53 

0.0 

28 

MAR 

1983 

7:48:28 

306 

34 

35 

123 

232 

35 

38 

145 

0.0 

28 

MAR 

1983 

13:13:9 

310 

39 

50 

206 

253 

39 

5! 

208 

0.0 

27 

MAR 

1983 

15:12:14 

3!  I 

220 

L. 

65 

1  32 

67 

8 

74 

0.0 

28 

MAR 

1983 

17:22:27 

312 

51 

8 

63 

312 

51 

8 

63 

0.0 

28 

MAR 

1983 

17:27:52 

313 

55 

12 

80 

233 

21 

35 

82 

0.0 

27 

MAR 

1983 

18:51 :39 

314 

41 

16 

72 

3 1  4 

41 

16 

72 

0.0 

28 

MAR 

1933 

8:13:24 

314 

109 

4 

51 

219 

196 

3 

70 

0.0 

27 

MAR 

1933 

14:47:45 

315 

19 

44 

85 

3 1  5 

19 

44 

85 

0.0 

28 

MAR 

1983 

17:44:47 

315 

59 

12 

108 

219 

40 

26 

113 

0.0 

28 

MAR 

1983 

22:13:3 

318 

20 

39 

80 

124 

40 

19 

87 

0.0 

27 

MAR 

1983 

23:9:24 

519 

37 

50 

195 

319 

37 

50 

195 

0.0 

28 

MAR 

1983 

22:13:54 

320 

23 

30 

75 

257 

47 

12 

77 

0.0 

27 

MAR 

1983 

21:30:19 

321 

46 

27 

1 55 

321 

46 

27 

156 

0.0 

27 

MAR 

1983 

22:29:27 

7^*1 

Ji.- 

101 

4 

85 

322 

101 

4 

85 

0.0 

27 

MAR 

1983 

21:44:15 

323 

57 

49 

191 

323 

37 

49 

191 

0.0 

28 

MAR 

1983 

8:57:41 

323 

33 

27 

97 

323 

33 

27 

97 

0.0 

27 

MAR 

1983 

7:2:57 

324 

21 

55 

124 

324 

21 

56 

124 

0.0 

27 

MAR 

1983 

22:17:11 

326 

36 

17 

72 

180 

30 

21 

74 

0.0 

28 

MAR 

1983 

0:14:15 

326 

35 

48 

182 

227 

39 

45 

192 

0.0 

23 

MAR 

1  983 

13: 13:0 

326 

42 

45 

197 

248 

39 

53 

205 

0.0 

28 

MAR 

1983 

13:25:3 

329 

29 

35 

108 

329 

29 

35 

108 

0.0 

28 

MAR 

1983 

13:55:0 

331 

122 

3 

45 

223 

48 

10 

83 

0.0 

28 

MAR 

1983 

17:44:52 

333 

55 

19 

113 

315 

68 

16 

129 

0.0 

28 

MAR 

1983 

8:53:  15 

J  j4 

52 

1  1 

6! 

124 

24 

26 

86 

0.0 

•27 

MAR 

1933 

17:44:44 

355 

27 

30 

90 

135 

25 

57 

109 

0.0 

23 

MAR 

1983 

0:13:52 

338 

43 

43 

198 

333 

43 

43 

1  98 

0.0 

28 

MAR 

1983 

17:57:9 

339 

61 

8 

53 

219 

38 

23 

96 

0.0 

27 

MAR 

1983 

18:21:41 

3-0 

45 

1  4 

75 

172 

147 

3 

89 

0.0 

27 

MAR 

1  983 

21:6:51 

340 

175 

4. 

49 

327 

327 

1 

53 

0.0 

27 

MAR 

1983 

22:19:23 

310 

53 

8 

91 

338 

73 

1 3 

113 

0.0 

28 

MAR 

1  933 

1 3 : 34: 30 

3-10 

27 

40 

1  14 

277 

30 

39 

120 

0.0 

27 

MAR 

1383 

22:18:28 

3-4 

43 

10 

54 

209 

55 

9 

60 

0.0 

27 

MAR 

1983 

20:33:41 

3*15 

44 

19 

90 

155 

40 

** 

94 

0.0 

28 

MAR 

1983 

19:7:5 

347 

24 

33 

89 

347 

24 

33 

89 

0.0 

23 

MAR 

1983 

17:39:13 

343 

48 

15 

95 

321 

93 

8 

100 

0.0 

27 

MAR 

1  583 

20:18:20 

552 

n 

Ci- 

29 

54 

165 

S3 

9 

78 

0.0 

28 

MAR 

1983 

8:57:51 

353 

54 

9 

55 

353 

54 

9 

65 

0.0 

27 

MAR 

1  983 

17:7:26 

35A 

57 

2i 

1 13 

251 

74 

23 

184 

0.0 

23 

MAR 

1933 

0:39:4 

355 

I  *1 

3 

1 13 

151 

25 

53 

159 

0.0 

23 

MAR 

1  333 

13:12:9 

395 

38 

53 

1  95 

291 

42 

13 

210 

0.0 

A- 1 3 


NORTH  BERING  WINTER  1933 


OATE 

TIME 

PM  I 

PA  1 

A 1 

FI 

PM2 

PA2 

A2 

F2 

UEL 

28 

MAR 

1983 

0:14:50 

355 

S3 

12 

128 

338 

27 

45 

128 

0.0 

27 

MAR 

1983 

0:43:21 

353 

34 

47 

175 

358 

34 

47 

176 

0.0 

27 

MAR 

1983 

7:3:51 

350 

37 

55 

212 

324 

37 

55 

214 

0.0 

28 

MAR 

1983 

22:8:55 

363 

32 

25 

95 

360 

32 

25 

56 

0.0 

27 

MAR 

1  983 

21  :4S:S 

3c  1 

s: 

1  1 

69 

229 

44 

21 

108 

0.0 

28 

MAR 

1  983 

1  ~  •  C-?  •  1C 

361 

E0 

9 

78 

35 1 

60 

9 

73 

0.0 

27 

MAR 

1  933 

22:54:23 

352 

98 

9 

113 

239 

94 

10 

119 

0.0 

27 

MAR 

1983 

1:3:12 

0  0  w 

37 

54 

21  1 

226 

45 

55 

257 

0.0 

27 

MAR 

1983 

21 : 10:3 

3S3 

202 

c 

54 

125 

45 

12 

61 

0.0 

28 

MAR 

1983 

17:25:39 

4C  4 

78 

8 

115 

353 

78 

8 

1  16 

0.0 

27 

MAR 

1983 

21  :54: 1 1 

357 

51 

20 

177 

348 

S3 

30 

207 

0.0 

27 

MAR 

1983 

22:25:0 

357 

39 

17 

104 

357 

39 

17 

104 

0.0 

28 

MAR 

1933 

0:32:23 

353 

38 

55 

220 

247 

42 

53 

231 

0.0 

27 

MAR 

1983 

18:29:27 

359 

120 

7 

115 

369 

120 

7 

116 

0.0 

27 

MAR 

1 9S3 

20:28:42 

369 

33 

40 

138 

359 

33 

40 

138 

0.0 

27 

MAR 

1933 

23:21 : 35 

359 

27 

31 

100 

369 

27 

31 

100 

0.0 

27 

MAR 

1983 

1  :4: 1 

370 

24 

54 

136 

370 

24 

54 

1 35 

0.0 

28 

MAR 

1983 

14:3:5 

370 

104 

8 

97 

370 

104 

8 

97 

0.0 

2S 

MAR 

1983 

23:22:13 

372 

26 

34 

115 

159 

48 

21 

121 

0.0 

28 

MAR 

1983 

13:55:33 

373 

71 

12 

93 

327 

57 

12 

100 

0.0 

27 

MAR 

1983 

21:13:34 

374 

S3 

15 

92 

374 

S3 

15 

92 

0.0 

27 

MAR 

1  983 

22:28:15 

374 

95 

5 

75 

372 

129 

5 

97 

0.0 

27 

MAR 

1933 

20:31 :43 

375 

82 

1 3 

121 

293 

77 

14 

122 

0.0 

28 

MAR 

1983 

22:13:15 

375 

23 

36 

94 

375 

23 

36 

94 

0.0 

27 

MAR 

1983 

15:53:0 

578 

40 

27 

115 

357 

55 

20 

121 

0..0 

27 

MAR 

1993 

0:48: 18 

379 

23 

55 

130 

182 

25 

55 

142 

0.0 

27 

MAR 

1983 

7:44:9 

379 

25 

58 

155 

379 

25 

58 

155 

0.0 

28 

MAR 

1933 

0:23:14 

330 

22 

51 

120 

1  13 

27 

59 

163 

0.0 

28 

MAR 

1  983 

13:54:6 

380 

178 

5 

59 

580 

178 

3 

59 

0.0 

27 

MAR 

1  983 

15:13:31 

381 

80 

12 

106 

275 

77 

15 

108 

0.0 

28 

MAR 

1983 

14:5:11 

334 

54 

19 

117 

334 

54 

19 

117 

0.0 

28 

MAR 

1983 

13:15:53 

335 

*1 

*.«- 

40 

95 

298 

24 

40 

100 

0.0 

27 

MAR 

1983 

1  : 3 : 5 1 

387 

36 

50 

186 

312 

40 

47 

199 

0.0 

23 

MAR 

1983 

22:8:41 

387 

25 

32 

86 

307 

30 

33 

107 

0.0 

27 

MAR 

1983 

21:17:49 

394 

50 

25 

148 

191 

53 

25 

1 55 

0.0 

23 

MAR 

1983 

22:13:20 

396 

25 

41 

109 

257 

31 

40 

130 

0.0 

28 

MAR 

1933 

8:42:42 

393 

45 

13 

67 

367 

99 

5 

63 

0.0 

23 

MAR 

1983 

0:22:2 

400 

51 

43 

262 

233 

55 

47 

292 

0.0 

28 

MAR 

1993 

17:42:19 

400 

41 

19 

98 

363 

37 

21 

99 

0.0 

28 

MAR 

1983 

U:4:4 

403 

80 

13 

1  14 

403 

80 

13 

1  U 

0.0 

27 

MAR 

1983 

18: 18:3 

404 

71 

13 

1  46 

252 

79 

17 

150 

0.0 

27 

MAR 

1983 

17:21:1 

407 

95 

5 

64 

21  1 

74 

8 

75 

0.0 

28 

MAR 

1  983 

17:38:41 

408 

31 

36 

125 

396 

40 

33 

U8 

0.0 

28 

MAR 

1933 

22:13:33 

409 

87 

10 

152 

287 

113 

8 

1  So 

0.0 

27 

MAR 

1933 

21  : 40 : 45 

410 

91 

10 

107 

410 

91 

10 

107 

0.0 

27 

MAR 

1  933 

17:28:34 

413 

413 

1 

59 

232 

125 

10 

147 

0.0 

27 

MAR 

1  983 

21:48:46 

415 

21 

35 

81 

185 

30 

36 

1  13 

0.0 

23 

MAR 

1  993 

11:45:39 

415 

35 

29 

120 

4  1  6 

35 

29 

120 

0.0 

23 

MAR 

1  933 

17:52:57 

415 

47 

17 

99 

253 

29 

35 

106 

0.0 

23 

MAR 

1  293 

13:15:34 

4  ^ 

**  fc.  i- 

2a 

51 

1  23 

1  87 

36 

47 

130 

0.0 

A- 14 


NORTH  BERING  WINTER  1983 


DATE 

TIME 

PM1 

PA1 

A1 

FI 

PM2 

PA2 

A2 

F2 

UEL 

28 

MAR 

1983 

0:38:58 

423 

30 

54 

173 

425 

30 

54 

173 

0.0 

27 

MAR 

1983 

23:8:25 

424 

41 

48 

206 

424 

41 

48 

206 

0.0 

27 

MAR 

1983 

7:28:41 

425 

39 

50 

207 

408 

46 

51 

245 

0.0 

28 

MAR 

1983 

13:23:7 

425 

42 

46 

204 

35 1 

43 

50 

227 

0.0 

28 

MAR 

1983 

22:17:28 

426 

54 

15 

97 

426 

54 

IS 

97 

0.0 

27 

MAR 

1983 

17:14:20 

427 

70 

7 

63 

240 

75 

7 

65 

0.0 

28 

MAR 

1983 

0:31:11 

428 

34 

53 

192 

428 

34 

53 

■192 

0.0 

28 

MAR 

1983 

13:18:47 

429 

21 

45 

101 

429 

21 

45 

101 

0.0 

27 

MAR 

1983 

15:16:19 

430 

58 

19 

125 

395 

62 

24 

165 

0.0 

28 

MAR 

1983 

19:47:25 

430 

52 

47 

255 

430 

52 

47 

256 

0.0 

28 

MAR 

1983 

7:45:47 

441 

36 

44 

168 

322 

33 

S3 

186 

0.0 

27 

MAR 

1S83 

21:40:47 

443 

56 

23 

191 

443 

55 

23 

191 

0.0 

27 

MAR 

1983 

17:45:33 

446 

33 

52 

18! 

446 

33 

52 

181 

0.0 

28 

MAR 

1983 

13:27:3 

452 

30 

46 

144 

452 

30 

45 

144 

0.0 

27 

MAR 

1983 

14:59:37 

453 

45 

34 

162 

4  S3 

45 

34 

162 

0.0 

28 

MAR 

1983 

18:5:38 

459 

38 

34 

140 

387 

57 

22 

148 

0.0 

27  MAR 

19e3 

16:0:0 

462 

>46 

4 

80 

413 

306 

c 

87 

0.0 

27 

MAR 

1983 

17:7:4 

472 

62 

13 

96 

154 

74 

1 S 

130 

0.0 

28 

MAR 

1983 

0:26:21 

473 

37 

5! 

193 

473 

37 

51 

198 

0.0 

28 

MAR 

1985 

19:9:5 

475 

43 

26 

127 

21  1 

60 

23 

162 

0.0 

28 

MAR 

1983 

0:9:27 

478 

29 

55 

163 

478 

29 

55 

163 

0.0 

27 

MAR 

1983 

17:4:19 

484 

80 

12 

166 

484 

80 

12 

1 65 

0.0 

27 

MAR 

1983 

15:12:32 

490 

96 

13 

147 

421 

87 

26 

245 

0.0 

27 

MAR 

1983 

0:48:42 

491 

31 

52 

163 

558 

36 

49 

198 

0.0 

28 

MAR 

1983 

13:55:50 

506 

77 

12 

101 

481 

115 

8 

105 

0.0 

27 

MAR 

1983 

7:28:7 

507 

32 

56 

189 

474 

32 

58 

196 

0.0 

28 

MAR 

1983 

17:26:23 

507 

26 

45 

121 

262 

32 

44 

146 

0.0 

27 

MAR 

1983 

20:25:39 

508 

144 

9 

145 

508 

1  4-4 

9 

U5 

0.0 

27 

MAR 

1993 

17:26:51 

516 

123 

5 

78 

196 

180 

3 

80 

0.0 

27 

MAR 

1983 

17:26:29 

535 

109 

31 

359 

535 

109 

31 

359 

0.0 

27 

MAR 

1983 

20:53:3 

557 

64 

10 

75 

123 

56 

19 

79 

0.0 

28 

MAR 

1983 

11:53:33 

550 

38 

40 

162 

461 

84 

17 

167 

0.0 

28 

MAR 

1933 

0:26:30 

557 

54 

53 

303 

414 

68 

50 

557 

0.0 

28 

MAR 

1983 

14:4:56 

563 

114 

7 

89 

553 

114 

7 

89 

0.0 

28 

MAR 

1983 

14:32:51 

578 

83 

1  1 

107 

555 

83 

12 

1  10 

0.0 

27 

MAR 

1983 

1  :3:2S 

534 

50 

52 

277 

261 

S3 

54 

301 

0.0 

28 

MAR 

1983 

14:6:29 

539 

57 

39 

240 

539 

57 

39 

240 

0.0 

27 

MAR 

19e3 

23:43:28 

598 

55 

4) 

242 

598 

55 

41 

242 

0.0 

28 

MAR 

1983 

19:41 : 14 

599 

106 

9 

123 

599 

106 

9 

123 

0.0 

28 

MAR 

1983 

0:27:39 

723 

56 

55 

327 

5)7 

57 

55 

330 

0.0 

28 

MAR 

1983 

7:54:47 

745 

61 

43 

288 

650 

82 

40 

359 

0.0 

SOUTH  CHUKCHI  WINTER  1983 


DATE 

TIME 

PM1 

PA  I 

A1 

FI 

PM2 

PA2 

A2 

F2 

UEL 

50 

MAR 

I9S3 

1  1  : 46:36 

198 

16 

45 

77 

80 

21 

37 

86 

0.0 

31 

MAR 

1963 

22:46:15 

205 

47 

7 

37 

102 

102 

1 

43 

0.0 

1 

APR 

1  9S3 

8:3:27 

205 

43 

15 

70 

114 

45 

21 

103 

0.0 

1 

APR 

1963 

8:16:22 

210 

46 

7 

42 

143 

77 

5 

53 

0.0 

1 

APR 

1963 

8:14:28 

213 

82 

4 

42 

1  10 

45 

6 

44 

0.0 

31 

MAR 

1953 

22:47:21 

216 

63 

5 

42 

165 

34 

15 

55 

0.0 

31 

MAR 

1583 

22:52:1 

217 

72 

4 

34 

195 

195 

1 

46 

0.0 

31 

MAR 

1983 

22:52:8 

217 

95 

3 

34 

190 

190 

1 

85 

0.0 

31 

MAR 

1983 

23:28:3 

220 

65 

8 

63 

117 

50 

120 

0.0 

1 

APR 

1963 

1  :45:37 

220 

32 

38 

133 

210 

96 

10 

134 

0.0 

1 

APR 

15e3 

8:19:53 

220 

147 

2 

53 

193 

70 

12 

93 

0.0 

31 

MAR 

1963 

22:52:16 

222 

222 

1 

31 

222 

2*7 

1 

31 

0.0 

31 

MAR 

1963 

15:54:46 

223 

25 

21 

53 

223 

25 

21 

S3 

0.0 

1 

APR 

1983 

8:21 :3S 

225 

94 

12 

128 

225 

94 

12 

123 

0.0 

30 

MAR 

1983 

1 1 : 43: 29 

226 

67 

4 

92 

226 

67 

4 

0"< 

0.0 

31 

MAR 

1983 

22:59:44 

226 

226 

1 

30 

215 

215 

1 

32 

0.0 

31 

MAR 

1983 

22:59:49 

226 

226 

1 

31 

218 

218 

1 

33 

0.0 

1 

APR 

1983 

8:19:10 

226 

50 

13 

70 

143 

53 

14 

80 

0.0 

31 

MAR 

1983 

22:59:25 

227 

227 

1 

29 

225 

225 

1 

30 

0.0 

1 

APR 

1983 

0:40:13 

230 

35 

30 

115 

209 

31 

34 

118 

0.0 

1 

APR 

1933 

8:14:20 

230 

56 

9 

59 

165 

59 

13 

85 

0.0 

29 

MAR 

1983 

16:3:20 

232 

66 

12 

93 

204 

59 

20 

135 

0.0 

1 

APR 

1983 

0:42:47 

232 

99 

3 

45 

232 

99 

3 

45 

0.0 

30 

MAR 

1983 

10:30:38 

233 

68 

S 

44 

194 

41 

21 

97 

0.0 

31 

MAR 

1983 

22:43:49 

233 

218 

3 

75 

233 

218 

3 

75 

.  0.0 

31 

MAR 

1983 

22:55:3 

233 

233 

1 

30 

213 

213 

1 

35 

0.0 

30, 

MAR 

1983 

8:31:14 

234 

44 

7 

59 

133 

51 

7 

66 

0.0 

1 

APR 

1983 

8:1  l  :7 

235 

50 

9 

55 

190 

43 

10 

56 

0.0 

29 

MAR 

1983 

17:15:34 

237 

24 

43 

109 

231 

24 

42 

112 

0.0 

31 

MAR 

1983 

22:45:34 

237 

49 

16 

84 

190 

53 

14 

86 

0.0 

31 

MAR 

1933 

22:48:34 

237 

69 

7 

53 

237 

69 

7 

53 

0.0 

31 

MAR 

1983 

22:52:22 

237 

237 

1 

39 

191 

75 

3 

46 

0.0 

1 

APR 

1983 

0:42:16 

243 

87 

7 

74 

U9 

66 

10 

82 

0.0 

30 

APR 

1983 

4:25:1 

243 

52 

10 

65 

124 

28 

23 

74 

5.3 

30 

APR 

1963 

5:15:15 

243 

97 

9 

152 

243 

97 

9 

152 

3.7 

1 

APR 

1983 

3:6:34 

245 

34 

49 

175 

245 

34 

49 

175 

0.0 

30 

MAR 

1983 

i 1 :55:S0 

245 

25 

40 

104 

230 

29 

34 

107 

0.0 

31 

MAR 

1983 

21 : 1 :36 

252 

31 

54 

177 

182 

35 

59 

218 

0.0 

31 

MAR 

1983 

22:44:32 

253 

53 

7 

53 

23! 

71 

5 

56 

0.0 

29 

MAR 

1983 

15:23:43 

254 

61 

5 

39 

249 

37 

10 

42 

0.0 

30 

MAR 

1983 

19:12:36 

254 

71 

27 

201 

254- 

71 

27 

201 

0.0 

31 

MAR 

1983 

13:22:43 

254 

72 

5 

43 

252 

72 

5 

43 

0.0 

30 

APR 

1983 

3:20: 16 

254 

36 

20 

83 

176 

34 

21 

85 

2.5 

30 

MAR 

1983 

7:21:13 

255 

23 

39 

94 

209 

25 

37 

99 

0.0 

30 

MAR 

1983 

10:41:14 

255 

86 

12 

118 

255 

86 

12 

118 

0.0 

31 

MAR 

1983 

22:45:47 

255 

100 

4 

63 

155 

67 

14 

105 

0.0 

1 

APR 

1993 

3:45:27 

255 

28 

55 

167 

205 

30 

54 

183 

0.0 

1 

APR 

1983 

4:9:24 

255 

27 

54 

1 55 

255 

27 

54 

155 

0.0 

1 

APR 

1983 

8  : 22 : -2 

295 

73 

7 

103 

255 

73 

7 

108 

0.0 

23 

MAR 

1983 

1  6  :  -IS  :  S3 

257 

40 

9 

67 

108 

15 

44 

70 

0.0 

SOUTH  CHUKCHI  WINTER  1983 


OATE 

TIME 

PM1  PA  1 

31 

MAR 

1993 

10:48:2 

258 

4S 

1 

APR 

1  SS3 

8:31  :7 

253 

55 

31 

MAR 

1993 

22:48:7 

260  i 

103 

1 

APR 

1983 

J  :47:S1 

260 

23 

1 

APR 

iges 

8:23:18 

260 

101 

30 

APR 

1933 

11:54:25 

260 

88 

30 

MAR 

1933 

9:32:0 

261 

74 

31 

MAR 

1983 

22:39:18 

261 

82 

30 

APR 

1983 

8:40:23 

261 

50 

30 

MAR 

1983 

12:26:5 

263 

49 

30 

MAR 

1963 

19:6:58 

265 

66 

31 

MAR 

1983 

11:19:10 

265 

82 

30 

APR 

1983 

7:21:30 

267 

55 

2 

APR 

1983 

13:23:39 

263 

62 

30 

MAR 

1983 

8:44:9 

269 

37 

1 

APR 

1983 

0:36:30 

269 

57 

30 

APR 

1983 

8:14:18 

270 

42. 

30 

APR 

1983 

11:12:16 

270 

67 

30 

APR 

1933 

4:59:47 

271 

56 

23 

MAR 

1983 

17:18:1 

272 

19 

1 

AFR 

1983 

8:24:32 

272 

57 

30 

APR 

1993 

5:24:51 

272 

27 

50 

AFR 

1983 

7:25:15 

273 

33 

30 

AFR 

1985 

8:59:52 

273 

72 

1 

APR 

1983 

8:21 : 42 

274 

104 

31 

MAR 

1983 

9:52:38 

275 

83 

51 

MAR 

1983 

22:47:10 

275 

121 

1 

APR 

1983 

3:44:30 

275 

26 

30 

APR 

1983 

1:22:55 

275 

39 

30 

APR 

1983 

6:13:28 

275 

55 

50 

AFR 

1983 

3:2:4 

276 

55 

30 

MAR 

1983 

12:26:51 

277 

70 

31 

MAR 

1983 

9:56:32 

278 

92 

29 

MAR 

1983 

13:17:49 

279 

1 03 

29 

MAR 

1983 

15:59:54 

279 

93 

29 

MAR 

1933 

16:6:23 

280 

26 

31 

MAR 

1983 

10:29:31 

280 

54 

50 

APR 

1  983 

3:13:3 

280 

280 

30 

MAR 

1933 

10:40:0 

282 

79 

1 

APR 

1983 

2:9:13 

283 

27 

30 

MAR 

1983 

16:57:39 

284 

73 

1 

APR 

1983 

4:53:23 

284 

53 

20 

APR 

1983 

12:18:52 

234 

1  1  9 

31 

MAR 

1933 

13:10:36 

295 

39 

29 

MAR 

1983 

15:22:55 

286 

65 

23 

!  MAR 

1983 

16:6:53 

235 

43 

30 

I  MAR 

1983 

8:44:17 

286 

1  1  9 

31 

MAR 

1933 

10:25:49 

296 

50 

31 

MAR 

1933 

10:50:35 

237 

241 

2C 

)  MAR 

1963 

17:2:  10 

238 

93 

A1 

FI 

PM2  PA2 

A2 

F2 

UEL 

16 

79 

231 

120 

5 

94 

0.0 

8 

81 

216 

55 

IS 

92 

0.0 

3 

44 

197 

33 

18 

67 

0.0 

47 

113 

128 

26 

50 

133 

0.0 

5 

53 

190 

65 

U 

95 

0.0 

6 

78 

1 88 

77 

8 

55 

8.4 

<6 

133 

231 

106 

1  1 

1 37 

0.0 

5 

50 

34 

41 

12 

58 

0.0 

16 

100 

261 

50 

16 

100 

1.4 

8 

64 

263 

49 

8 

64 

0.0 

5 

43 

266 

65 

5 

43 

0.0 

5 

48 

253 

47 

13 

74 

0.0 

9 

59 

2-2 

87 

10 

100 

2.9 

20 

145 

253 

71 

20 

1 92 

0.0 

13 

62 

202 

23 

39 

95 

0.0 

14 

97 

263 

57 

14 

97 

0.0 

18 

97 

252 

1 33 

4 

99 

2.9 

8 

82 

“215 

67 

20 

143 

10.3 

1  1 

75 

271 

55 

1  1 

75 

4.2 

45 

93 

272 

19 

45 

95 

0.0 

14 

87 

272 

57 

14 

87 

0.0 

29 

87 

cue 

23 

39 

121 

4.4 

15 

55 

198 

CC 

25 

59 

7.  1 

15 

105 

243 

71 

12 

117 

7.7 

6 

90 

265 

95 

8 

97 

0.0 

4 

37 

275 

63 

4 

37 

0.0 

5 

81 

275 

121 

6 

81 

0.0 

53 

1  44 

275 

26 

53 

144 

0.0 

10 

48 

133 

30 

23 

76 

6.7 

8 

Go 

275 

55 

8 

65 

3.8 

10 

72 

186 

20 

50 

106 

9.8 

10 

82 

277 

70 

10 

82 

0.0 

6 

66 

209 

55 

12 

91 

0.0 

5 

75 

225 

112 

S 

82 

0.0 

13 

U9 

279 

98 

13 

149 

0.0 

43 

120 

260 

26 

43 

120 

0.0 

20 

1  13 

280 

54 

20 

1  13 

0.0 

1 

50 

230 

n-> 

c 

23 

71 

5.2 

8 

73 

232 

79 

8 

73 

0.0 

53 

153 

234 

32 

55 

185 

0.0 

9 

70 

125 

60 

13 

82 

0.0 

8 

64 

240 

27 

33 

107 

0.0 

3 

41 

1  U 

114 

1 

44 

6.5 

5 

43 

285 

89 

5 

43 

0.0 

9 

75 

213 

50 

17 

91 

0.0 

14 

72 

257 

«•  c 

35 

85 

0.0 

3 

51 

286 

1  19 

3 

61 

0.0 

1  4 

74 

236 

53 

14 

74 

0.0 

3 

87 

287 

241 

3 

87 

0.0 

7 

65 

232 

63 

1  1 

79 

0.0 

A- 1 7 


SOUTH  CHUKCHI  WINTER  1983 


DATE 

TIME 

PM1 

PA1 

A1 

FI 

PM2 

PA2 

A2 

F2 

UEL 

30 

MAR 

1933 

12:28:40 

290 

104 

4 

69 

202 

67 

8 

78 

0.0 

31 

MAR 

1  SS3 

10:49:3 

290 

95 

7 

79 

155 

50 

1 1 

107 

0.0 

30 

APR 

1933 

10:53:10 

290 

48 

8 

SI 

255 

51 

10 

6! 

4.8 

30 

APR 

1  933 

14:43:53 

290 

87 

IS 

141 

269 

77 

13 

148 

.5 

30 

MAR 

1983 

18:5:40 

291 

1 55 

3 

6! 

214 

63 

12 

79 

0.0 

31 

MAR 

1963 

13:23:5 

291 

109 

4 

43 

284 

109 

4 

51 

0.0 

1 

APR 

1983 

8:22:57 

293 

69 

1  1 

86 

249 

175 

4 

90 

0.0 

30 

APR 

1983 

1:14:53 

253 

54 

9 

55 

284 

27 

27 

80 

4.1 

30 

APR 

1983 

7:52:14 

293 

46 

20 

102 

2S3 

46 

20 

1 02 

2.2 

30 

APR 

1933 

7:54:53 

294 

45 

15 

107 

294 

45 

15 

107 

3. 1 

30 

APR 

1983 

12:17:53 

294 

65 

9 

66 

261 

87 

10 

95 

8.0 

23 

MAR 

1 SS3 

13:17:25 

295 

16 

46 

73 

213 

25 

35 

97 

0.0 

31 

MAR 

1933 

13:22:38 

255 

74 

5 

42 

293 

75 

5 

45 

0.0 

31 

MAR 

1983 

13:19:22 

296 

38 

19 

78 

295 

38 

19 

78 

0.0 

31 

MAR 

1983 

22:48:22 

296 

48 

1  1 

57 

253 

52 

17 

97 

0.0 

2 

APR 

1983 

22:45:31 

296 

84 

8 

76 

251 

83 

8 

73 

0.0 

30 

APR 

1983 

6:31:33 

256 

85 

5 

70 

254 

74 

5 

81 

3 . 6 

31 

MAR 

1983 

13:22:34 

297 

63 

5 

4 1 

295 

63 

5 

42 

0.0 

1 

APR 

1983 

9:53:2 

297 

117 

9 

120 

297 

1  17 

9 

120 

0.0 

29 

MAR 

1983 

17:22:49 

298 

29 

42 

135 

292 

31 

42 

148 

0.0 

30 

MAR 

1983 

6:24:39 

298 

55 

IS 

124 

161 

62 

19 

151 

0.0 

30 

MAR 

1983 

10:40:12 

299 

96 

10 

107 

259 

96 

10 

107 

0.0 

29 

MAR 

1983 

.7:53:3 

300 

55 

8 

63 

500 

55 

8 

63 

0.0 

2 

APR 

1983 

21 : 34:23 

300 

30 

39 

124 

229 

30 

43 

137 

0.0 

29 

MAR 

1983 

17:19:36 

301 

25 

46 

126 

234 

27 

46 

157 

0.0 

2 

APR 

1983 

21:33:46 

301 

50 

36 

1  17 

296 

59 

31 

150 

0.0 

30 

APR 

1983 

6:11:2 

301 

43 

24 

125 

301 

43 

24 

125 

8.2 

31 

MAR 

1933 

22:47:27 

302 

44 

16 

76 

302 

44 

16 

76 

0.0 

1 

APR 

1983 

8:21:53 

302 

68 

25 

185 

254 

92 

24 

237 

0.0 

30 

MAR 

1  S83 

8:54:1 1 

503 

21 

48 

107 

303 

21 

48 

107 

0.0 

31 

MAR 

1933 

13:22:43 

303 

80 

5 

47 

303 

80 

5 

47 

0.0 

31 

MAR 

1983 

21:9:32 

304 

47 

57 

282 

243 

47 

57 

283 

0.0 

31 

MAR 

1933 

22:55:8 

304 

60 

12 

102 

280 

71 

13 

121 

0.0 

2 

APR 

1933 

22:32:17 

304 

48 

21 

108 

222 

57 

19 

117 

0.0 

31 

MAR 

1983 

9:31 :55 

505 

47 

29 

U3 

305 

47 

29 

143 

0.0 

30 

APR 

1983 

1  :S2: 16 

305 

41 

1  1 

49 

222 

50 

9 

50 

8.5 

31 

MAR 

1983 

12:52:32 

306 

98 

6 

67 

153 

89 

10 

99 

0.0 

2 

APR 

1983 

15:35:43 

306 

7*1 

44 

U7 

306 

32 

44 

U7 

0.0 

il 

APR 

1933 

21:35:11 

306 

45 

40 

191 

231 

46 

41 

198 

0.0 

30 

APR 

1933 

3:3:55 

306 

84 

5 

64 

207 

108 

4 

70 

4.0 

30 

APR 

1983 

3:27:3 

306 

45 

30 

14  4 

306 

45 

30 

1  44 

7.3 

29 

MAR 

1993 

15:53:55 

308 

25 

30 

105 

224 

49 

22 

1  19 

0.0 

30 

MAR 

1983 

17:0:5 

309 

174 

C 

43 

309 

174 

(■ 

43 

0.0 

29 

MAR 

1  283 

16:0:4 

310 

44 

12 

70 

310 

44 

12 

70 

0.0 

30 

APR 

1933 

7:11:14 

310 

40 

16 

78 

310 

40 

16 

73 

S.3 

31 

MAR 

1993 

13:22:29 

31  1 

1  13 

3 

40 

309 

112 

3 

40 

0.0 

1 

APR 

1983 

1:31:14 

31  1 

177 

119 

251 

174 

n 

C 

127 

0.0 

1 

APR 

1383 

7:1:19 

31  1 

31 

77 

31  1 

31 

77 

0.0 

30 

APR 

1993 

14:11:43 

313 

54 

13 

80 

3 1  3 

54 

13 

80 

.4 

30 

APR 

1  933 

1 : 57 : 4R 

316 

50 

20 
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APPENDIX  B 


USE  OF  EXTREME  VALUE  STATISTICS  IN  ENGINEERING  ANALYSIS 

Random  phenomenae,  such  as  ice  loads,  wind  speed  and 
wave  heights  can  be  quantified  using  statistical  techniques. 
Probability  density  functions  describe  the  chances  of  any  value 
occurino  in  a  given  sample  rather  than  predicting  a  certain  value. 
Most  statistical  methods  deal  with  the  analysis  of  the  behaviour 
of  the  random  variable  around  the  mean  (i.e.  Traffic  Volumes, 
Mortality  Rates,  I.Q.).  However,  certain  classes  of  problems 
demand  that  the  rare  events  be  investigated.  This  is  particularly 
true  in  engineering  where  the  rare  event  is  normally  the  design 
event  (i.e.  fastest  wind,  highest  wave,  largest  load).  The 
statistical  theory  of  extreme  values  has  been  developed  to 
describe  the  probabilities  of  these  rare  events. 


Gumbel  [1]  stated:  "The 


statis' 


theory  o 


VJ  UlllU  ^  .A.  w  w  V.  -  —  —  —  —  — - -  -  > 

extreme  values  is  to  explain  observed  extremes...  and  to  fcrecas 
(future)  extremes".  Extreme  value  functions  naturally  depend  on 
the  initial  or  underlying  distributions.  Figure  1  illustrates 
this  with  an  initial  distribution  Fx(x)  the  random  variable  X,  and 
the  distribution  of  the  extremes  Gy  (yft)  where; 

N 


FX(X) 


»  probability  that  any  X  will  be  less  than  x 


Gy  (yn)  =  probability  that  any  Y^  will  be  less 


than  y 


where  X  «■  random  variable, 

x  =•  a  specific  value  which  X  might  take 

Y  -  Max  of  N  samples  of  X  (also  a  random 

N  variable) 

v  -  a  specific  value  which  Y  might  take 


por  y  to  be  less  than  y  ,  all  N  X  must  be  less  than  yn< 

follows  that; 


GY(yn)  -  ifx(*>1  -  yn 
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rr  t~t, 


Thus,  if  the  initial  distribution  of  a  Function  is  known  exactly, 
then  the  extreme  value  function  can  be  derived  exactly.  in 
engineering  practice  this  is  not  usually  possible.  Small  errors 
ip^  estimates  of  the  initial  distribution  are  compounded  (to  the 
N”  power)  in  the  extremes.  It  is  precisely  due  to  this  problem 
that  statistical  theory  of  extreme  values  was  developed.  The 
early  works  on  extreme  value  statistics  are  due  to  Vcn  Kises 
(1922),  Frechet  (1927)  and  Fischer  and  Tippett  (1928).  Gumbels' 
landmark  paper  and  book  [1,2]  give  both  the  mathematical  deriva¬ 
tion  and  engineering  uses  of  extreme  distributions .  There  are 
three  main  types  of  extreme  value  functions,  (see  Table  1)  origin¬ 
ally  classified  by  Fischer  and  Tippett,  and  which  have  also  come 
to  be  known  by  the  mathematicians  who  studied  them;  Type  I 
(Gumbel),  Type  II  (Frechet)  and  Type  III  (Weibull).  All  three 
distributions  have  a  largest  value  and  smallest  value  form.  Only 
the  largest  value  forms  will  be  discussed  here. 

Asymptotic  Distributions 

Th*  reasons  that  three  specific  extreme  value 
distributions  have  been  defined  is  that  they  are  all  asymptotic 
distributions.  For  most  situations  of  practical  importance,  for  a 
wide  range  of  initial  distributions,  the  extremes  as  defined  in 
Equation  (1)  will  asymptotically  approach  one  of  the  three 
distributions,  as  the  value  of  N  grows  large. 

Type  I  (Gumbel) 

The  Type  I  distribution  arises  when  dealing  with 
unlimited  random  variables  whose  initial  distributions  tail  off 
"in  an  exponential  manner";  as  for  example, 

F(x)  -  1  -  e'Xx  (2) 


which  is  the  negative  exponential  distribution.  The  normal  and 
gamma  distribution  are  cf  this  general  type.  The  Type  I,  defining 
the  largest  of  many  independent  random  variables  with  a  common 
exponential  type  of  upper  tail,  is  written  as; 


—  oc(  y-u ) 

ysr> 

—  a 

=*  e  —  °°  <  y  <  °° 

(3) 

where 

a 

=  measure  of  dispersion 

u 

=  mode  of  y 
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Type  II  (Frechet) 

The  Type  II  also  represents  the  largest  of  many  random 
variables,  but  which  are  limited  on  the  left  to  zero  (i.e. 
positive  only)  and  have  a  tail  which  falls  off  as: 

Fx)  -1-6  (I)"  x  >  0  ,4) 


The  Cauchy  ana  log  normal  follow  this  form.  The  Type  I  and  Type 
II  have  the  same  relationship  as  the  normal  and  log  normal  distri¬ 
butions.  If  a  variable  is  log  normally  distributed  then  its 
logarithm  is  normally  distributed,  and  if  a  variable  is  Type  II 
then  its  logarithm  is  Type  I.  This  can  be  seen  in  the  probability 
paper  used  to  plot  Type  I  and  Type  II.  Both  Type  I  and  II  use  the 
same  scale  for  probability,  but  while  Type  I  uses  a  linear  verti¬ 
cal  scale.  Type  II  uses  a  logarithmic  scale  (see  Figures  2  and  3). 
Note  also  that  if  x  in  Equation  4  is  replaced  with  lnx,  then  equa¬ 
tion  4  can  be  transformed  into  equation  2.  All  this  leads  to  the 
Type  II  distribution  with  the  form; 


v  k 

F(y)  -  e-(-)  y  >  0  (5) 

y 

where  \>  -  a  parameter  between  the  mode  and  median 

k  -  a  measure  of  dispersion. 


Type  III  (Weibull) 

The  Type  III  distribution  describes  the  largest  of  many 
random  variables  which  are  limited  to  some  maximum  value  w,  and 
fall  off  in  a  manner  such  that  near  w; 

F ( x )  -  1  -  c ( w  -  x)k  x  <  w  k  >  0  (6) 


The  distribution  of  the  largest  values  is: 
k 

-fw-y) 

F(y)  -  e  lw-uj  y  <  w  (7) 
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where  w  -  maximum  possible  value  of  y 
u  -  mode  of  y 
k  -  measure  of  dispersion 


Weibull  [1934]  actually  used  the  distribution  to  study  smalles 
rather  than  largest  values.  He  studied  material  strength  i 
tension  and  fatigue.  While  his  justification  for  the  distribution 
was  purely  empirical,  Freudenthal  [1951,  1956]  has  given  physical 
reasons  for  its  practicability. 
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Examde 


The  maximum  annual  mean  hourly  wind  speed  for 


London  airport  for 

the 

years 

1939-1961  inclusive  are  as 

follows ; 

36,  37,  45, 

50, 

39,  33, 

37,  35,  41,  52, 

41,  58, 

39,  46,  39, 

42, 

45,  36, 

55,  32,  43,  34, 

39 

To  plot  the  extreme  values  the  data  must  be  ranked 
as  follows: 


RANK 

m 

VELOCITY 

(MPH) 

P  ( v )  =  1  -  m 

n  +  1 

Y  =  - 1  n  (  - 1  n  (  P  )  ) 

1 

58 

.958 

3.15 

2 

55 

.917 

2.44 

3 

52 

.875 

2.01 

4 

50 

.833 

1.70 

5 

46 

.792 

1.46 

6 

45 

.750 

1.25 

7 

45 

.708 

1.06 

8 

43 

.667 

.90 

9 

42 

.625 

.75 

10 

41 

.583 

.62 

11 

41 

.542 

.49 

12 

39 

.500 

.37 

13 

39 

.458 

.25 

14 

39 

.417 

.13 

15 

37 

.375 

.02 

16 

37 

.333 

-.09 

17 

37 

.292 

-.21 

18 

36 

.250 

-.33 

19 

36 

.208 

-.45 

20 

35 

.167 

-.58 

21 

34 

.125 

-.73 

22 

33 

.083 

-.91 

23 

32 

.042 

-1.16 

v  =  41.39  a  -  6.87 
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The  data  can  be  plotted  against  any  of  the  three  distributions. 
This  is  best  done  by  first  linearizing  the  equations  as  follows: 


Tvoe  I 


-a(  x-u ) 

F ( x )  -  e~e 


(3) 


can  be  converted  to: 


-ln(-ln( F(x) ) )  ■  ox  -  au 


which  is  linear  in  -ln(-ln( F(x) ) )  and  x.  The  plotting  paper  used 
in  Figure  2  is  specially  constructed  so  that  data  fitting  Equation 
3  will  be  a  straight  line. 


Type  I I 


(5) 


can  be  converted  to: 


-ln(-ln( F(x ) ) )  -  klnx  -  klnv 


which  is  linear  in  -ln(-ln(Fx) ) )  and  Inx.  The  same  kind  of 
plotting  paper  can  be  used  to  plot  Type  II  if  log  x  is  used. 
Figure  3  uses  a  logarithmic  scale  instead. 

Type  III 


F(  x ) 


k 


(7) 


can  be  converted  to: 


-ln( -ln( F( x ) ) )  -  kln(w-u)  -  kln(w-x) 
which  is  again  linear.  Figure  4  plots  w-x  on  a  logarithmic  scale. 


B-8 


RETURN  PERIOD 


e 

E 


cn 

ci 


cn 


in 


(3 


>~ 

|_ 


CQ 

cr 

DQ 

o 

C V 
CL 


in 

(V 


EtnSinsinsin 

ld  ld  ui  ui  v  n~  cn 


C HdW )  Q33dS  QNIM 


s  in 

cn  cm 


B-9 


FIGURE  EV.2  TYPE  I  EXTREME  VALUE  PLOT  OF  WIND  SPEEDS 


RETURN  PERIOD 


[I] 


FIGURE  EV.3  TYPE  II  EXTREME  VALUE  PLOT  OF  WIND  SPEEDS 


RETURN  PERIOD 


(  KdW ) a 3 jcS  aNIM-SZ 
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FIGURE  EV.4  TYPE  III  EXTREME  VALUE  PLOT  OF  EXTREME  WINDS 


COMPARISON  OF  DIFFERENT  DISTRIBUTIONS 


This  example  illustrates  a  very  important  point 
concerning  extreme  value  statistics.  All  three  distributions  fit 
the  data  reasonably  well  over  the  range  of  the  data.  If  one  were 
interested  in  predictions  within  the  range  of  probabilities  given, 
the  choice  would  be  inconsequential.  The  problem  arises  when 
extrapolating  to  longer  time  spans  than  the  available  data.  The 
100  year  and  1000  year  expected  extreme  winds  are: 


100  yr.  1000  yr . 

Type  I 
Type  II 
Type  III 


63 

75 

62 


75 

105 

73.5 


When  extrapolating  2,  3  or  even  4  orders  of  magnitude  (i.e.  from  2 
ship  weeks  of  ice  loads  data  to  50,000  ship-weeks  of  operations) 
the  problems  of  divergence  of  the  distributions  is  even  greater. 

The  choice  of  the  distribution  must  rely  on  the 
applicability  of  assumptions  governing  the  distributions  to  the 
particular  process  governing  the  data.  The  key  distinguishing 
aspects  of  the  distributions  are: 

Type  I  -  unlimited  additive  ( averaging ) process 

Type  II  -  unlimited  multiplicative  (sequential) 
process 

Type  III  -  limited  variable 

The  extreme  value  method  is  a  valuable  predictive  tool,  but  care 
must  be  taken  to  get  sufficient  data  and  to  use  the  method  wisely. 
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THE  THREE  DISTRIBUTIONS  COMBINED 

It  has  been  shown  (Gnedenko  1943)  that  the  three  Fisher 
Tippett  distributions  are  the  only  stable  asymptotic  forms 
for  extreme  values  (except  for  improper  distribution). 
Further,  Jenkinson  (1955)  shewed  that  ail  three  forms  could 
be  expressed  as  special  cases  of  a  more  general  distribution. 
For  the  purposes  of  clarity  this  most  general  distribution 
will  be  called  the  Jenkinson  distribution.  It  is  written  as 
follows : 


11/C 

-  1  -  C  (x  -  Al)  (8) 

F ( x )  =  e  L  A2 

This  reduces  to  a  type  I,  II  or  III  depending  on  the  value  of 
C,  as  follows; 


Type  I ;  C  «•  0 


-1/A2  (x  -  Al) 
-e 

F  ( x )  -  e 


where  1/A2  -  a  ,  Al  -  U 


Type  II;  C  <  0 


-1/C 

-f  Al/C _ 1 

F  (  x  )  =  e  Lx-  ( Al  -  A2/C)  J 

where  X  .  =•  Al  -  A2/C  (lower  limit,  usually  0) 

mi  n 

-1/C  -  k 


A2/C 


V 


Type  III: 


C  >  0 


F(x)  -  e 


C/A2  ( (A2/C  +  Al)  -  x] 


1/C 


where  X  =  W  =  Al  +  A2/C  (uDDer  limit) 

max  / 


1/C 


Al  =*  U 


These  interrelationships  allow  fitting  data  to  equation  8  and 
then  determining  both  the  type  and  parameters  of  the  extreme 
value  equation. 
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APPENDIX  C 


TABLE  I 


SUMMARY  OF  DEPLOYMENTS,  DATA  SETS  AND  ICE  CONDITIONS 


TITLE  &  DATE 

LOCATION 

ICE 

TYPE 

NO  OF 
EVENTS 

Beau  1  or t  Summer  82 

Sep  28  -  Oct  16 

100-150  nm  north  of  Prudhoe  Bay 
in  the  Alaskan  Beaufort  Sea 

MY 

167 

S  Bering  Winter  83 

Mar  29  -  Mar  26 

Transit  from  St. Paul  Is.  to  the 
west  end  of  St. Lawrence  Is.  in 
the  Bering  Sea 

FY 

1  73 

N  Bering  Winter  83 

Mar  27  -  Mar  28 

Transit  from  St.  Lawrence  Is.  to 
the  Bering  Strait  in  the  Bering  Sea 

FY 

291 

S  Chukchi  Winter  83 

Mar  29  -Apr  2 

Apr  28  -  May  2 

Transit  from  the  Bering  Strait  to 
Point  Hope  in  the  Chukchi  Sea  and 
return 

FY  ,MY 

299 

N  Chukchi  Winter  83 

Apr  3  -  Apr  27 

Round  trip  transit  Point  Hope  to 
W’ainwrighc  in  the  Chukchi  Sea, 
operation  off  Wainwright 

FY ,  MY 

513 

Antarctic  Summer  89 

Jan  9  -  Jan  13 

McMurdo  Sound,  break-in  to  McMurdo 

Base 

FY 

309 

Beau i  or  l  Summe r  89 

Nov  18  -  Dec  1 

Operation  between  Barter  Is.  and 
Barrow  in  the  Beaufort  Sea,  transit 
through  the  Chukchi  Sea  to  the 

Bering  Strait 

FY ,  MY 

337 

SUBSETS  OF  KNOWN  MULTIYEAR  EVENTS 

N  Cliukch  i  Winter  8  5  MY 
Apr  3  -  A p  r  27 

North  Chukchi  Sea  of!  Wainwright 

MY 

6  7 

Beau)  ■  r t  Summe r  89  MY 

Beaufort  and  Chukchi  Seas 

MY 

32 

Nov  I  2  -  Dec  I 
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TABLE  2 


SUMMARY  OF  COMBINED  DATA  SETS 


TITLE 

COMBINED  FROM 

ICE 

TYPE 

NO  OF 
EVENTS 

Known  Multiyear 

Beaufort  Summer  82 

N  Chukchi  Winter  83  MY 

Beaufort  Summer  84  MY 

MY 

266 

Heavy  Mixed  FY  f.  MY 

Beaufort  Summer  82 

S  Chukchi  Sea  83  N  Chukchi  Sea  83 
Beaufort  Summer  84 

FY ,  MY 

1017 

Known  First  Year 

S  Bering  Winter  83 

N  Bering  Winter  83 

Antarctic  Summer  84 

FY 

723 

Summer  B<-nufort  Sea 

Beaufort  Summer  82 

Beaufort  Summer  84 

mostly  MY 

504 

W  inter  Chukchi  Sea 

S  Chukchi  Winter  83 

N  Chukchi  Winter  83 

FY ,  MY 

398 

W  i  n  1 1  r  Be  ring  Sea 

S  Bering  Winter  83 

N  Bering  Winter  83 

FY 

812 
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APPENDIX  D 


SUMMARY  OF  TABLES  OF  THE  COEFFICIENTS  OF  THE  3  PARAMETER 
CURVE  FIT  TO  THE  EXTREME  VALUE  DISTRIBUTIONS 


TABLE  D-I 

REGRESSION  COEFFICIENTS  FOR  THE  3  PARAMETER 
EXTREME  VALUE  DISTRIBUTION  FOR 
HIGHEST  AVERAGE  PRESSURE  ON  A  SINGLE  SUB-PANEL 


REGRESSION  COEFFICIENTS 
NUMBER  A1  A2 


TITLE 

OF  EVENTS 

C 

MPa 

psi 

MPa 

psi 

Beaufort  Summer  82 

167 

-.063 

3.10 

450 

1.43 

207 

S  Bering  Winter  83 

173 

-.071 

1.27 

184 

.35 

51 

N  8ering  Winter  83 

241 

.026 

1.99 

289 

.58 

84 

S  Chukchi  Winter  83 

299 

-.218 

2.01 

291 

.43 

62 

N  Chukchi  Winter  83 

513 

-.198 

2.50 

363 

.71 

103 

Antarctic  Summer  84 

310 

.042 

1.90 

276 

.52 

75 

Beaufort  Summer  84 

337 

.000 

1.97 

286 

.83 

121 

N  Chukchi  Winter  83  MY 

67 

-.236 

2.95 

428 

.95 

138 

Beaufort  Summer  84  MY 

32 

-.247 

1.96 

284 

.92 

134 

Known  Multi-Year 

266 

.000 

2.89 

419 

1.52 

221 

Heavy  Mixed  FY  and  MY 

1017 

-.128 

2.37 

343 

.92 

133 

Known  First-Year 

723 

.042 

1.71 

248 

.63 

91 

Summer  Beaufort  Sea 

504 

-.116 

2.22 

322 

1.06 

154 

Winter  Chukchi  Sea 

812 

-.193 

2.28 

330 

.65 

94 

Winter  Bering  Sea 

398 

.000 

1.56 

226 

.66 

96 

TABLE  D-2 


REGRESSION  COEFFICIENTS  FOR  THE  3  PARAMETER 
EXTREME  VALUE  DISTRIBUTION  FOR 
HIGHEST  AVERAGE  PRESSURE  ON  FOUR  SUB-PANELS 

REGRESSION  COEFFICIENTS 


NUMBER 

A1 

A2 

TITLE 

OF  EVENTS 

C 

MPa 

psi 

MPa 

psi 

Beaufort  Summer  82 

148 

-.204 

1.31 

190 

.60 

87 

S  Bering  Winter  83 

154 

.020 

.56 

81 

.21 

31 

N  Bering  Winter  83 

228 

-.092 

.86 

125 

.29 

42 

S  Chukchi  Winter  83 

271 

-.104 

.88 

128 

.28 

40 

N  Chukchi  Winter  83 

460 

-.200 

1.14 

165 

.33 

48 

Antarctic  Summer  84 

283 

.101 

.88 

127 

.23 

34 

Beaufort  Summer  84 

289 

.055 

.86 

124 

.43 

62 

N, Chukchi  Winter  83  MY 

63 

-.322 

1.31 

190 

.46 

67 

Beailfort  Summer  84  MY 

24 

.071 

.88 

128 

.52 

75 

Known  Multi-Year 

266 

.001 

1.21 

175 

.81 

118 

Heavy  Mixed  FY  and  MY 

897 

-.168 

1.05 

152 

.41 

60 

Known  First-Year 

665 

.005 

.78 

113 

.28 

40 

Summer  Beaufort  Sea 

-.136 

.89 

129 

.53 

77 

Winter  Chukchi  Sea 

-.149 

1.05 

153 

.35 

51 

Winter  Bering  Sea 

-.043 

.72 

104 

.33 

48 

0-2 


TABLE  D-3 


REGRESSION  COEFFICIENTS  FOR  THE  3  PARAMETER 
EXTREME  VALUE  DISTRIBUTION  FOR 
HIGHEST  TOTAL  FORCE  ON  THE  PANEL 


REGRESSION  COEFFICIENTS 


NUMBER 

A1 

A2 

TITLE 

OF  EVENTS 

C 

MN 

LT 

MN 

LT 

Beaufort  Summer  82 

167 

-.133 

1.12 

112 

.62 

62 

S  Bering  Winter  83 

173 

.130 

.96 

96 

.38 

38 

N  Bering  Winter  83 

241 

-.048 

1.03 

103 

.45 

45 

S  Chukchi  Winter  83 

299 

-.052 

.87 

87 

.37 

37 

N  Chukchi  Winter  83 

513 

-.016 

1.23 

123 

.63 

63 

Antarctic  Summer  84 

310 

-.005 

.93 

93 

.30 

30 

Beaufort  Summer  84 

337 

.000 

.77 

77 

.46 

46 

N  Chukchi  Winter  83  MY 

67 

-.163 

1.40 

140 

.60 

60 

Beaufort  Summer  84  MY 

32 

-.028 

.73 

73 

.52 

52 

Known  Multi-Year 

266 

.001 

1.08 

108 

.75 

75 

Heavy  Mixed  FY  and  MY 

1017 

-.043 

1.04 

104 

.60 

60 

Known  First-Year 

723 

-.048 

.96 

96 

.36 

36 

Summer  Beaufort  Sea 

504 

-.114 

.83 

83 

.50 

50 

Winter  Chukchi  Sea 

812 

.001 

1.06 

106 

.63 

63 

Winter  Bering  Sea 

398 

.001 

.95 

95 

.45 

45 

D-3 


TABLE  0-4 


REGRESSION  COEFFICIENTS  FOR  THE  3  PARAMETER 
EXTREME  VALUE  DISTRIBUTION  FOR 
HIGHEST  TOTAL  FORCE  ON  A  FRAME 


REGRESSION  COEFFICIENTS 


TITLE 

NUMBER 

A1 

A2 

OF  EVENTS 

C 

MN 

LT 

MN 

LT 

Beaufort  Summer  82 

167 

.001 

.56 

56 

.37 

37 

S  Bering  Winter  83 

173 

-.197 

.22 

22 

.07 

7 

N  Bering  Winter  83 

241 

-.239 

.36 

36 

.11 

11 

S  Chukchi  Winter  83 

299 

-.084 

.39 

39 

.10 

10 

N  Chukchi  Winter  83 

513 

-.127 

.46 

46 

.16 

16 

Antarctic  Summer  84 

310 

.028 

.33 

33 

.10 

10 

Beaufort  Summer  84 

337 

.005 

.30 

30 

.13 

13 

N  Chukchi  Winter  83  MY 

67 

-.033 

.62 

62 

.26 

26 

Beaufort  Summer  84  MY 

32 

-.127 

.34 

34 

.18 

18 

Known  Multi-Year 

266 

.001 

.53 

53 

.35 

35 

Heavy  Mixed  FY  and  MY 

1017 

.005 

.39 

39 

.26 

26 

Known  First-Year 

723 

-.145 

.31 

31 

.10 

10 

Summer  Beaufort  Sea 

504 

-.205 

.38 

38 

.19 

19 

Winter  Chukchi  Sea 

812 

-.145 

.43 

43 

.14 

14 

Winter  8ering  Sea 

398 

-.188 

.29 

29 

.11 

11 

0-4 


TABLE  D-5 


REGRESSION  COEFFICIENTS  FOR  THE  3  PARAMETER 
EXTREME  VALUE  DISTRIBUTION  FOR 
HIGHEST  TOTAL  FORCE  ON  A  STRINGER 


REGRESSION  COEFFICIENTS 
NUMBER  A1  A2 


TITLE 

OF  EVENTS 

C 

MN 

LT 

MN 

LT 

Beaufort  Summer  82 

167 

.001 

.89 

89 

.55 

55 

S  Bering  Winter  83 

173 

.005 

.41 

41 

.25 

25 

N  Bering  Winter  83 

241 

.234 

.75 

75 

.30 

30 

S  Chukchi  Winter  83 

299 

.001 

.73 

73 

.30 

30 

N  Chukchi  Winter  83 

513 

-.020 

.98 

98 

.36 

36 

Antarctic  Summer  84 

310 

.095 

.69 

69 

.23 

23 

Beaufort  Summer  84 

337 

.169 

.60 

60 

.36 

36 

N  Chukchi  Winter  83  MY 

67 

.001 

1.13 

113 

.55 

55 

Beaufort  Summer  84  MY 

32 

.210 

.52 

52 

.32 

32 

Known  Multi-Year 

266 

.001 

.89 

89 

.54 

54 

Heavy  Mixed  FY  and  MY 

1017 

.006 

.81 

81 

.41 

41 

Known  First-Year 

723 

.170 

.64 

64 

.29 

29 

Summer  Beaufort  Sea 

504 

-.066 

.66 

66 

.39 

39 

Winder  Chuk  chi  Sea 

812 

.001 

.87 

87 

.37 

37 

Winter  Bering  Sea 

398 

.165 

.60 

60 

.31 

31 

D-5/D-6 


APPENDIX  E 


DESIGN  OF  SHELL  PLATING  ON  ICEBREAKERS 


Failure  of  shell  plating  can  be  defined  in 
several  ways: 

°  initial  yielding 

°  unacceptable  deformation  (permanent  set) 

°  rupture  (after  large  deformation) 

°  fatigue 

Initial  yielding  is  not  normally  considered  as 
failure  in  icebreakers.  No  visible  damage  results  from 
yielding.  As  well,  design  against  initial  yielding  would 
be  unnecessarily  costly.  Fatigue  considerations  are 
normally  dealt  with  through  proper  selection  of  steel 
grades . 


For  modern  icebreakers,  constructed  of  high 
strength  and  highly  ductile  steels,  the  design  condition  is 
some  level  of  permanent  set  that  is  considered  acceptable. 
For  these  steels  the  rupture  strength  is  very  high  and  is 
of  less  interest. 

Post  yield  behavior  of  plating  depends  on  a 
number  of  parameters  and  mechanisms : 

°  plastic  hinge  formation  at  edges 
°  plastic  hinge  formation  at  center 
°  membrane  action  due  to  deflection 
°  lateral  deflection  of  plate  edges 
°  "pulling  out"  of  plastic  hinges  under 
high  membrane  loads 
°  plate  aspect  ratio 
°  plate  slenderness  ratio 
°  load  shape 

°  load  relaxation  due  to  deformation 

The  complexity  of  the  plastic  response  of  the 
plate,  together  with  scarcity  of  experimental  data,  has  so 
far  resulted  in  a  diverse  range  of  opinions  and  equations 
in  the  literature.  This  report  will  discuss  the  various 
approaches  and  propose  a  suitable  basis  for  design  of  shell 
plating. 

°  stress  strain  model 
°  yield  criteria 
°  ultimate  strain  at  rupture 


E-l 


Timoshenko  [1]  studied  long  flat  places  with 
built  in  and  simply  supported  edges.  All  of  his  solutions 
considered  elastic  action  only.  His  solutions  considered 
the  membrane  stresses  resulting  from  laterally  restrained 
supports.  This  work  provided  the  basis  for  later 
plasticity  analysis. 


Clarkson  [2]  considered  the  elastic  and  plastic 
resconse  of  rigid  clamped  plates  of  various  aspect  ratios, 
up  to  the  formation  of  a  single  plastic  hinge  at  the  edge 
(including  membrane  action).  Clarkson's  equation  (for  zero 
aspect  ratio)  is: 


P 


4 .56 


’  4/3  1/3  ' 

cr  /  E 

L  y 


4/3 

{  t/s  ) 


(  F  .  1  ) 


p  *  applied  load  (capacity) 


a  m  yield  strength 
E  =  young's  Modulus 


t  -  plate  thickness 

s  =  plate  span  (frame  spacing) 

The  resulting  permanent  deformation  was  equal  to  these 
caused  by  the  welding  process  during  fabrication. 

Johansson  [3]  considered  a  three  plastic  hinge 
mechanism  as  the  design  condition.  He  assumed  that  membr¬ 
ane  stresses  would  not  be  important  and  that  (practically) 
no  permanent  deformation  would  result.  His  strength 
formula  is: 


P 


4  <j 


ks 


(  F  .  2  ) 


Hooke  (5)  based  on  the  experimental  results  of 
Hooke  and  Rawlings  (4),  extended  Clarkson's  work  on  fixed 
edged  plates  of  non-zero  aspect  ratio.  In  (4)  it  is  noted 
that  the  tests  were  unable  to  prevent  lateral  movements  of 
the  plate  boundaries. 

Jones  (6]  examined  finite  plastic  deflection  of 
rectangular  plates  with  simply  supported  and  fully  clamped 
edges.  Using  a  rigid  perfectly  plastic  stress  model  and 


E-2 


assuming  no  lateral  deflections  at  the  boundaries,  he 
produced  the  following  formula  for  the  fully  clamped  case: 


P  =  P 


2 

— 

1 

2  fC  +  ( 3-2C)  3 

p  • 

1  +  - 

-  i 

c 

3 

it  j 

1(3-0  /  J 

(F.3) 


where  W  <  t 


!  =  12  a 

c  v 


*h  3  +  a 


(3  h i n c  e  load) 


C  =  a 

k  'I  3  +  a  -  a 
a  =  s/a  (aspect  ratio) 
a  =  height  of  plate 
W . =  permanent  deflection 

(Note  that  (F.3)  reduces  to  (F.2)  when  a  and  W  are  ecual  to 
0)  . 

For  the  simply  supported  case: 


P  -  P 


1  + 


rw 


t) 


2  rc  +  ( 3-2  c ) 

l  (3  -  C) 


(F.4) 


W  <  t/2 


o  r 


P  = 


p  4W 
c  t 
2 


1  +  C  (  C  -  2)  fl  - 
3  -  C 


2 

t  ^ 
2 

12W  , 


W  >  t/2 


(F.5) 


E-3 


for  plates  of  zero  aspect  ratio  this  reduces  to: 


p  * 8  (I)  (!) 


(F.6) 


Coburn  et  al  [7]  took  equation  (F.3)  for  W/t  -  1,  a  -  0  to 
obtain 


p-8,  g) 


(  F  .  7  ) 


and  suggested  a  plate  design  (inverting  ( F . 7 ) )  of; 


t  -  .353  s  /  P/<r, 


Coburn  et  al  also  reports  that  Clarkson's  equation  (F.l)  if 
extended  to  3  plastic  hinges  would  become  approximately 


c  4/3  /  1/3  rt'|  4/3  ( F .  8 ) 

P  -  9.12  L  y  /  E  IsJ 


Hughes  [8]  in  a  very  comprehensive  paper,  developed  a 
solution  for  plastic  plate  response  based  on  an 
elastic-perf ectly  plastic  stress  model  and  allowing  for  the 
freedom  of  the  edges  to  pull  in  (Figure  1).  Hughes  defines 
the  slenderness  ratio  0  as; 


<3  -  s/t  Toy/E 


(F.9) 
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Figure  1  ASSUMED  CONDITION  of  PLATING  for  ANALYSIS  by  HUG 


ps 


ps 


Figure  2  IDEAL  MEMBRANE  BEHAVIOUR  of  PLATE 


ps  ps 


Figure  3  MEMBRANE  BEHAVIOUR  in  the  PRESENCE  of  SPRINGS 


Hughes  explains  that  for  sturdy  plates  (0  <  2.4)  the 
deflections  under  load  will  be  quite  small  and,  because 
of  the  occurrence  of  edge  pull-in,  that  membrane  effects 
are  small  even  after  yielding.  It  is  noted  that  plates 
icebreakers  are  very  sturdy  (i.e.  s  ■  400  mm,  t  -  36  mm, 
-  360  will  result  in  0  -  .46!)  and  that  edge  pull-in  ecu 
be  very  important.  As  a  result  membrane  stresses  would 
small.  Using  a  non-dimensional  load  parameter  Q; 


on 

<3 


Q 


pE/ffy2 


(F.10) 


Hughes  developes  a  3  part  expression  for  the 
load/deformation  equation: 


Q  -  Qy  +  Kw  UQq  +  AQ1  •  RJ 

in  which; 


(F.ll) 


Qy 

*Q0 

AQ1 


Wp 

Wp 

0 


nondimensional  load 

increase  in  load  up 
the  edge  hinge 

slope  of  Qvs  Wp/WpQ 


at  initial  yield 
to  full  formation  of 

after  edge  hinge 


W 

l  >  1 

<1-(1-Rw)3)1/3  rw  <  1 

permanent  deformation 


permanent  deformation  at  completion  of  edge 
hinge  formation 


.  O702t 
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tr »- 


Equation  {F.ll)  is  based  on  snail  deflection  theory  and  is 
valid  for  Wpg  <  Wp  <  t.  Equation  (F.ll)  considers  the 
steel  to  follow  Hencky-VonMi ses  yield  criteria.  The 
equation  includes  the  effects  of  non  zero  aspect  ratios. 
The  complete  equations  are  as  follows: 


Qy 


1  +  0.5 


(F.12) 


(F.13) 
(  F  .  14  ) 


(F.15) 


or 

p  ■  3-“  %  (if 


(F.16) 
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For  a 


.2,  0 


.46,  and  Wp 


t,  (F.ll)  reduces  to 


■  5-53  %  (|)2  |F-13) 

To  examine  the  edge  pull-in,  consider  first  the  limiting 
condition  of  membrane  behaviour.  Figure  2  illustrates 
ideal  membrane  behaviour  with  pinned  ends.  Ratzlaff  and 
Kennedy  [9]  provide  a  clear  discussion  of  the  ideal 
membrane  situation  (as  well  as  a  very  informative 
discussion  of  many  aspects  of  the  problem) .  It  can  easily 
be  shown  that  for  the  case  shown  in  Figure  2 


P  -  8  W  T  ( F.19 ) 


T  »  Membrane  Tension 


This  is  true  for  any  stress  condition  (elastic  or  plastic) 
as  long  as  the  membrane  tension  T  is  constant  throughout 
the  plate.  The  increase  in  length  due  to  the  sag  is  given 
by  X; 


2 

X  -  8  W_  (  F.  20  ) 

3  s 

again,  regardless  of  the  state  of  stress. 
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which  foe  a  -  360  KPa  is; 


For  deflections  above  Wy,  the  tension  becomes  constant  at 
the  yield  value  (T  =■  cyt/(  1-vp^ )  )  which  results  in  a  load 
vs.  deflection  equation  of: 


8  <j 


where  vp  »  plastic  Poissons  ratio 
-  .5 


It  is  necessary  to  examine  the  assumption  of  edge  pull-in. 
To  do  this  the  neighbouring  plate  is  considered  as  a  spring 
(as  in  Figure  3).  To  get  the  value  of  the  spring  constant, 
assume  that  the  plates  to  either  side  are  semi-infinite  and 
edge  loaded  over  a  length  of  a.  Timoshenko  [10]  gives  the 
solution  for  this  condition.  The  deflection  of  the  center 
of  the  load  is: 


5  -  2  a  (a  log  a/2) 
re  E 


(  F  .  28  ) 


which  means  that  the  stiffness  (for  plates 
on  both  sides  combined)  is: 


K  -  re  E  t _  (  F  .  29  ) 

e  4  log  (a/2) 
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The  stiffness  of  the  plate  itself  is  determined  using 
equations  (F.20  and  F.21). 


K  =>  T/X  •  a 
P 

2 


(1-v  ) 


2 

a  w 

3  s 


K  =  E  t  •  a 
P  2 

s  (1-m  )  ( F . 3  0 ) 


It  is  possible  to  combine  Kq  and  K  and  simply  consider  the 
ideal  membrane  case  (Figure* 2)  but  with  a  more  flexible 
plate  of  the  combined  stiffness  Kc,  is 


K  =  1/ ( 1/K  +  1/K  ) 

c  p  e 


=  K  (1/(1  +  K  /K  ) ) 
P  P  e 


»  k  e  ( f . 31 ) 

p 


0  =  1/(1  +  1.4  a  log  a/2  (F.32) 

s 


E- 1 1 


by  determining  T  from  Kc 
deflection  ecuation: 


•  X,  a  modified  load  vs. 


P 


(1-v  ) 


6 


2 

+  (W) 


(F.33) 


Mote  that  (F.33)  is  only  (24)  times  9.  Essentially  E  can 
be  replaced  by  E  ©  to  simulate  the  edge  springs.  Again  the 
elastic  solution  is  only  valid  up  to  the  yield  stress,  at 
which  point  the  deflection  is; 


W 

y 


3 

y 

8  e  e 


s 


(F.34) 


9  -  1/p.  +  1.4  a  log  a/2j 

Given  -  360 

a  =  1 . 2  ra 
a/s  =  4 

The  onset  of  full  membrane  action  would  take  place  at  a 
deflection  of 


w  *  .036  •  s 


(F.35) 
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Any  deflection  greeter  than  3.6%  of  the  frame  scan  would 
result  in  complete  plasticity  and  therefore  complete 
membrane  action.  This  simplifies  the  problem  greatly 
because  T  becomes  constant  at  <r  t/(l-v  “).  Again  utiliz 
(F.19)  the  load/deflecticn  equation  for  a  full  plastic 
membrane  is: 


ing 


at  1 

p  *  8  W  •  y 


4  2 

(1-v  )  4  s  +  (4Ws) 

P 


cr 


p  =  10.67  \.sj  a 


4  1  + 


(F.36) 


(4  W/s) 


for  W/s  -  -0 . 1 


P  -  0.99  a  t 

y  s 


(F.37) 


It  is  worth  noting  that  when  the  full  plastic  membrane 
stress  is  activated,  the  plastic  hinge  at  the  edge  of  the 
plate  is  completely  pulled  out  resulting  in  no  edge  moment, 
and  the  ideal  membrane  condition  upon  which  equation  (F.37) 
depends.  Note  also  that  (F.36)  is  almost  identical  to 
(F.6),  which  would  be  expected,  since  both  involve  both 
membrane  action.  In  fact,  equation  (F.5)  probably  best 
describes  the  action  of  plates  subject  to  large  deflections 
and  has  the  advantage  of  considering  non  zero  aspect 
ratios.  As  a  result  it  is  suggested  that  a  suitable  design 
criteria  would  be  some  ratio  of  permanent  set  (W)  to  frame 


E-13 


spacing.  The  required  plate  thickness  would  be: 


t 


4.1 7x10" 2  s 


fP  'I 


V  YJ 


e 


w  =  allowable  permanent  set 
P  =  design  ice  load  (MPa) 
s  =  frame  spacing  (mn) 
o  y  a  yield  stress  (MPa) 

0  =  function  of  plate  aspect  ratio 


(  F .  38  ) 


(F.39) 


2  2 
H  3  +  0  -  (3) 


1  +  C  (C  -  2)  'i 
3  -  C  J 


C  -  3  H  3  +  0  -  0) 


0  is  Tabulated  in  Table  1. 


(  F .  40  ) 


0 


e 


o 

.1 

.2 

.25 

.33 

.5 

1 


3 

2.99 

2.96 

2.93 

2.86 

2.71 

2 
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No  tine  the  limited  ranee  of  8,  a  simpler  and  still  accurate 
formula  could  be: 


t  -  0.125  s  P  (±) 

a  lw  J  (  F .  4 1  ) 

y 


0.03  <  W  <  0.2 

s 
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